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DEFINITION 

Bubble  Parameter 

C' 

d 

Empirical  'drag  coefficient,'  based  on 
experimental  data,  used  to  calculate 
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Load  Factor 
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SUMMARY 


The  experimental  results  presented  herein  satisfy  one  major  goal  in 
demonstrating  the  technical  feasibility  of  two-phase  LMMHD  (liquid-metal 
MHD),  i.e.,  operating  an  MHO  generator  at  power  densities  equal  to  or  above 
that  anticipated  for  practical  power  systems.  Power  densities  of  up  to 
32  MWe/m3  and  efficiencies  higher  than  0.6  at  high  void  fractions  were 
attained  for  a small  20-kWe  generator.  Slip  ratio  data,  and  more-extensive 
pressure  distribution  and  voltage  profile  data  are  also  given. 

Barium  has  been  identified  as  an  attractive  additive  for  the  generation 
of  foams  with  NaK,  and  bismuth  is  a back-up  for  barium.  A guiding  hypothesis 
for  understanding  the  phenomena  contributing  to  foamability  in  liquid  metal 
systems  has  been  developed,  and  basic  data  on  surface  tension  obtained. 

A computer  code  has  been  developed  that  will  calculate  the  electrical 
end  loss  for  any  insulating  vane  geometry  and  magnetic  flux  density  distri- 
bution, a significant  accomplishment.  Initial  studies  with  this  code  have 
shown  that  generator  efficiencies  in  excess  of  0.8  are  attainable. 

Initial  measurements  of  local  flow  parameters  in  a NaK-nitrogen 
two-phase  liquid-metal  MHD  generator,  using  hot-film  and  resistivity 
probes,  are  reported.  Although  the  results  are  preliminary,  there  is 
evidence  that  large  electromagnetic  force  densities  have  a significant 
effect  on  the  size  and  distribution  of  the  gas  bubbles  or  voids. 

All  results  are  very  encouraging  for  the  development  of  large  high- 


efficiency  LMMHD  generators. 


I.  INTRODUCTION 


The  two-phase  LMMHD  (liquid-metal  MHD)  generator  combines  experimental 
and  analytical  studies  designed  to  increase  understanding  of  the  physical 
processes  occurring  in  two-phase  generators  and  to  establish  the  limits 
of  generator  performance.  For  eight  years  this  program  has  been  sponsored 
by  the  Power  Branch  of  the  Office  of  Naval  Research. 

The  generator  experiments  are  conducted  with  an  ambient-temperature 
NaK-N2  test  facility,  consisting  of  a two-phase  mixer,  MHD  generators 
magnet,  and  required  auxiliary  loop  equipment.  The  original  facility, 
used  in  all  tests  through  July  1975,  has  been  described  in  Pittenger  et  al 
(1972).  The  revised  facility,  completed  in  1977,  is  described  in  Petrick 
et  al  (1977).  Earlier  experiments  on  generators  different  geometries 
have  been  performed  and  the  results  presented  in  previous  annual  reports 
[Amend  et  al  (1973),  Petrick  et  al  (1975),  Petrick  et  al  (1976),  Petrick  et 
al  (1977),  and  Petrick  et  al  (1978)]. 

Analytical  models  of  the  generator  have  been  developed  to  aid  in 
the  basic  understanding  of  the  flow  characteristics  in  the  generator  and 
to  provide  a (one-dimensional)  computer  code  capable  of  predicting  approxima- 
tions of  generator  performance  from  such  independent  input  parameters  as 
flow  rate,  magnetic  field  strength,  and  geometry.  This  computer  code  and 
the  semiempirical  technique  by  which  the  bubble  parameter  (churn-turbulent 
drag  coefficient)  used  in  the  determination  of  the  interfacial  slip  ratio 
is  characterized  have  been  reported  in  detail  [Amend  et  al  (1973),  Petrick 
et  al  (1970)]. 

The  channel  contours  used  in  generator  experiments  are  compared  in 
Fig.  1.1.  The  first  channel  (LT-1)  was  designed  for  a slip  ratio  of  unity, 
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because  the  actual  slip  ratio  was  not  known,  but  higher  slip  ratios  were 
observed  experimentally,  and  it  was  noted  that  the  liquid  velocity  was  not 
constant;  it  decreased  along  the  generator  [Amend  et  al  (1973)].  A constant- 
area  channel  (LT-2)  was  then  tested,  and  it  was  observed  that  the  velocity 
increased  along  the  channel  [Amend  et  al  (1973),  Petrick  et  al  (1975)]. 
The  second  diverging  channel  (LT-3)  has  a smaller  angle  than  the  first 
channel,  as  shown  in  Fig.  1 . 1 , to  yield  close-to-constant  liquid  velocity 
at  the  design  point.  As  has  been  described  [Petrick  et  al  (1976),  Petrick 
et  al  (1977)]  improved  generator  performance  was  obtained  with  LT-3,  and 
the  liquid  velocity  was  almost  constant.  LT-3  was  also  tested  with  the 
revised  facility,  and  these  tests  demonstrated  the  expected  reduction  in 
slip  ratio  at  higher  liquid  velocities  (flow  rates)  [Petrick  et  al  (1978), 
Fabris  et  al  (1977a)].  The  new  high-pressure,  high-velocity  channel  (LT-4) 
was  designed  to  utilize  the  full  NaK-nitrogen  loop  pressure  and  flow-rate 
capabilities,  and  test  results  are  reported  herein.  The  first  channel 
for  the  high-temperature  sodium-nitrogen  test  facility  (HT-1)  has  been 
fabricated  and  tested  at  elevated  temperatures  under  open-circuit  conditions 
[Dunn  et  al  (1979)]. 

The  experimental  and  analytical  aspects  of  the  two-phase  LMMHD  generator 
program,  with  the  exception  of  the  high-temperature  generator  experiments, 
are  described  in  Section  II.  The  LT-4  test  results,  presented  in  detail 
in  Section  III,  include  the  effects  of  liquid  velocity  on  efficiency, 
power  density,  and  slip  ratio,  and  data  on  pressure  distributions 
and  voltage  profiles.  Further  positive  results  on  foams  and  surface-active 
agents  are  reported  in  Section  IV.  The  results  of  an  extensive  analytical 
study  of  end  effects  are  presented  in  Section  V,  and  some  initial  data 
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on  local  measurements  are  summarized  In  Section  VI.  The  conclusions 
drawn  from  the  work  presented  herein  are  given  in  Section  VII,  and  a 
brief  review  of  the  significant  results  from  the  eight-year  program  as 
a whole  is  in  Appendix  0. 
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II.  PROGRAM 


II.l  NaK-N2  Generator  Experiments 


The  NaK-Ng  ambient-temperature  flow  facility  used  for  these  tests 
has  been  described  previously  [see  Petrick  et  al  (1977)].  The  contour  of 
the  new  channel,  LT-4,  is  compared  with  the  other  channel  contours  tested 
in  Fig.  1.1.  This  channel  was  designed  for  the  higher  pressu,*°s  and  flow 
rates  of  the  revised  facility  and,  thus,  uses  certain  novel  constructor, 
features.  The  channel  contour  is  set  by  a fi lament -wound,  fiberglass  liner 
that  forms  the  insulating  walls  and  is  relatively  easy  to  replace.  The 
liner  was  wound  on  a mandrel  machined  to  the  desired  channel  contour,  and 
the  exterior  of  the  liner  was  then  machined  to  the  shape  required  to  accept 
the  copper  electrode  inserts  and  copper  compensating  bars.  The  copper 
structure  absorbs  all  of  the  mechanical  stress,  the  liner  serves  solely  as 
the  insulating  channel  walls.  Because  the  liner  is  one  piece,  i.e.,  no 
glued  joints,  it  is  easy  to  seal  the  electrodes  and  end  flanges  to  it  with 
0-rings;  thus,  no  leakage  occurs. 

The  internal  contour  of  the  channel  was  designed,  as  for  all  previous 
two-phase  MHO  generator  channels,  using  the  generator  portion  of  the  system 
code  described  in  Appendix  A of  Petrick  et  al  (1970).  In  that  code,  the 
mass,  momentum,  and  energy  balances  are  solved  in  a stepwise  fashion  along 
the  channel.  The  liquid  velocity  is  held  constant  (constant  ut  gives 
constant  induced  voltage  for  constant  B,  and  minimizes  the  internal  circula- 
ting currents  and  losses),  an  ideal  gas  and  perfect  heat  transfer*  between 


* See  Petrick  et  al  (1978)  or  Fabris  and  Pierson  (1978)  for  experimental 
confirmation  of  this  heat  transfer. 
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the  gaseous  and  liquid  phases  are  assumed,  and  the  liquid  and  gas  properties 
are  calculated  at  each  step.  The  channel  dimensions  and  powers  are  calcu- 
lated for  specified  inlet  conditions  and  outlet  void  fractions.  The  code 
can  handle  different  working  fluids  (a  version  exists  for  steam  in  place  of 
an  ideal  gas),  different  models  for  the  two-phase  electrical  conductivity, 
and  either  no  slip  between  the  phases  or  a slip  model  with  a bubble  parameter, 
Cq  , that  is  an  input  va.  ’able. 

The  contour  calculate4  for  the  generator  design  conditions  is  shown  in 
Fig,  1 1.1  (note  the  expanded  channel  height  scale).  The  bubble  coeff- 
icient Cj  = 2000  nr1  was  used  for  the  design,  based  on  the  analysis 
of  the  most-recent  LT-3  experiments.  The  extra  expense  of  machining  the 
calculated  contour  onto  the  liner  mandrel  was  rejected  as  unnecessary;  t>"js 
the  two  straight-line  segments  shown  in  Fig.  II. 1 were  used.  (Previous 
channels  used  only  one  straight-line  segment,  but  they  had  much  less 
divergence.)  The  straight-line  segments  were  extended  to  the  ends  of  the 
liner,  0.038  m (1.5  in.)  outside  of  the  electrodes. 

The  diagnostics  for  the  experiment  were  improved  for  LT-4.  The 
eight  pressure  transducers  used  previously  along  the  electrode  side  of 
the  channel  and  end  regions  were  replaced  by  seventeen  transducers,  and 
eight  transducers  were  added  to  the  insulating  side  in  the  end  regions 
just  outside  of  the  magnet  pole  face  (four  at  the  inlet,  four  at  the  exit), 
ihe  six  void  fraction  windows  spaced  along  the  channel  were  retained,  and 
two  each  added  in  the  entrance  and  exit  sections.  Three  sets  of  eight 
voltage  probes  spaced  between  the  electrodes,  i.e.,  across  the  insulating 
wall,  were  incorporated,  as  used  previously  in  LT-2.  The  sets  are  located 
0.0127,  0.1937,  and  0.3620  m downstream  of  the  beginning  of  the  electrodes. 


The  experimental  data  measured  with  LT-4  channel  are  discussed  in 
Section  III,  and  the  significant  parameters  for  all  of  the  experimental 
runs  from  October  1977  to  December  1978  are  summarized  in  Appendix  A. 
Significant  results  include  the  higher  efficiencies  obtained  at  higher 
liquid  flow  rates,  the  high  power  densities  measured  (comparable  to  or 
larger  than  those  anticipated  for  large,  practical  LMMHD  generators), 
ar.d  the  reduction  in  slip  into  the  liquid  mass  flow  rate  is  increased. 


« • 


Vi 


II. 2 Generator  as  Part  of  a Thermal  Engine 

In  practically  all  thermal  engines,  only  a gas  (vapor)  is  used  as  the 
working  fluid  th'it,  by  its  expansion,  converts  thermal  to  mechanical  (or 
electrical)  energy.  In  the  two-phase  LMMHD  generator,  the  mechanical 
energy  of  the  liquid  is  also  converted  to  electricity.  However,  this  does 
not  represent  the  conversion  of  thermal  energy,  and  the  liquid  mechanical 
energy  must  be  replaced  elsewhere  in  the  cycle.  Accordingly,  instead  of 
the  generator  isentropic  efficiency 

ni  * electrical  power  output/isentropic  power  input,  (11.1) 

a generator  "equivalent  turbine"  efficiency 


nt  * electrical  power  out-ideal  liquid  mech.  pump  work 
isentropic  power  in-ideal  liquid  rnech.  pump  work 


(11.2) 


can  be  defined  to  characterize  the  LMMHD  generator  as  a thermal  energy 
converter.  Note  that  is  useful  in  interpreting  experimental  genera- 
tor performance  because  the  liquid  work  is  eliminated,  a positive  value 
shows  that  the  electrical  power  out  exceeds  the  ideal  liquid  pump  work. 
However,  is  used  in  system  performance  studies  because  the  liquid 
pump  work  and  flow  losses  are  properly  accounted  for.  The  isentropic  power 
input  minus  the  ideal  liquid  mechanical  pump  work  represents  the  contribu- 
tion resulting  from  the  presence  of  the  gas.  Note  that  the  contribution  of 
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the  gas  is  much  larger,  because  of  the  heat  transferred  to  the  gas  from  the 
liquid  metal,  than  just  the  change  in  gas  enthalpy,  which  is  small  because 
of  the  small  temperature  change  across  the  generator.  The  specific  expres- 
sions for  each  of  these  quantities  are  developed  in  Appendix  B. 

The  equivalent  turbine  efficiency  can  be  more  easily  interpreted  in  the 
following  manner.  The  isentropic  power  input,  AHmix,  can  be  written  as 

A^mix  * ^ Vmixdp  = (^1  + Vg)AP»  (I I *3) 

where  ap  * pjn  - pout  is  the  pressure  differences  across  the  generator, 

• » t 

and  V1  is  the  effective  average  volumetric  flow  rate  such  that  V'ap  = / Vdp 

• I 

for  each  component.  In  Eq.  II. 3,  V^Ap  is  the  ideal  (isentropic)  liquid 

• I 

mechanical  pump  work,  and  VgAp  represents  the  mechanical  energy  contribution 
from  the  gas.  Note  that  the  latter  is  not  the  enthalpy  change  of  the  gas 
(most  of  the  energy  actually  comes  from  the  temperature  change  of  the  liquid), 
but  is  the  extra  isentropic  power  available  from  the  mixture  because  of  the 
presence  of  the  gas  and  the  heat  transfer  from  the  liquid  to  the  gas.  If  the 
weighted  average  mixture  void  fraction  is  defined  as 

« Vg/(Vg  + Vi,.  (II *4) 

it  follows  that 

• I 

ViAp/isentropic  power  input  -•  1 - a'  (II. 5) 

and 

• I 

VgAp/isentropic  power  input  c o'  (1 1.6) 

The  resultant  relations  between  "equivalent  turbine"  efficiency,  isentropic 
efficiency,  and  weighted  void  fraction  are  shown  in  Figs.  1 1.2  and  II. 3 for 
assumed  isentropic  efficiencies  of  0.90  and  0.50,  respectively.  One  notes 
that,  when  o'  exceeds  1 - ni  the  electrical  power  output  exceeds  the  pump 
work  input  and  the  "turbine  efficiency"  is  positive.  Further,  at  high  void 
fractions  approaches  n^. 
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Figure  II. 2 "Equivalent  turbine"  efficiency  vs  weighted  mixture  void 
fraction  for  a generator  isentropic  efficiency  of  0.9 
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Experimental  data  for  the  LMMHD  generator  experiments  can  be  recast 
in  terms  of  nt»  Table  II. 1,  shows  that  the  net  electrical  output  exceeds 
the  net  pump  work  input  at  higher  void  fractions,  with  only  a slight 
departure  from  isothermal  operation.  The  "equivalent  turbine"  efficiencies 
are  low  for  these  cases  because  of  the  low  generator  efficiencies  and 
low  average  void  fractions,  but  are  consistent  with  Fig.  1 1. 3.  With 
operation  at  higher  average  void  fractions  (such  as  foam  flow),  much  higher 
isentropic  and  turbine  efficiencies  are  expected,  Fig.  1 1.2.  Note  that 
the  temperature  difference  along  the  generator  is  not  measured  directly, 
but  is  calculated  from  the  experimental  data  using  Eq.  B.8,  Appendix  B. 


TABLE  II .1 

Experimental  Data  for  "Equivalent  Turbine"  Efficiencies 

Run  8/18/  8/18/  2/17/  2/18/  7/12/  11/29/ 

76-2  76-3  77-6  77-7  78-1  78-3 


Average  a 

Net  Pump 
work  (W) 

Electrical 
power  out  (W) 


2566  2544 


3205  3764 


aT  generator  0.4 

(K) 


4214  4554  21600 


0.464  0.430  0.524  0.530  0.606 


12500 


0.487 


0.147  0.197  0.399  0.476  0.457 


0.346 
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1 1 . 3 Foams 

The  most  significant  advance  of  the  year  in  the  surfactant  program  was 
the  discovery  that  barium  markedly  stabilizes  bubble  lifetime  and  promotes 
foaming  in  a NaK  alloy.  A number  of  other  systems  were  also  tested,  but 
none  yielded  beneficial  effects  approaching  those  of  barium.  Bubble  tests 
and  surface  tension  measurements  proceeded  in  parallel,  with  the  surface 
tension  data  yielding  the  scientifically  interesting  information  that 
barium  raises  the  surface  tension  at  concentrations  effective  for  bubble 
stability.  Exploration  of  this  finding  leads  to  the  development  of  a 
guiding  hypothesis  for  understanding  the  phenomena  contributing  to  foam- 
ability  in  liquid  metal  systems.  This  hypothesis  invok»s  a tendency  for 
attractive  interatomic  interactions  to  enhance  film  elasticity  and,  hence, 
bubble  stability.  The  extent  of  the  applicability  of  this  hypothesis  can 
only  be  determined  by  more  detailed  laboratory  studies.  Details  of  all 
aspects  of  the  work  are  reported  in  Section  IV. 

NaK  was  chosen  for  these  initial  experiments  because  it  allows 
room  temperature  tests,  eliminating  the  special  problems  associated  with 
handling  hot  liquid  metals.  Thus,  the  concept  of  using  surface-active 
agents  in  MHD  generators  can  be  evaluated  more  rapidly  and  inexpensively 
with  NaK,  the  existing  NaK  facility  is  available  for  future  flow  tests, 
and  a generator  of  known  performance  properties  with  pure  NaK  will  be 
available  to  observe  the  changes  with  a foam. 
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1 1. 4 local  Measurements* 

The  operation  of  two-phase  LMMHD  generators  depends  on  the  processes 
occurring  at  the  local  level,  i.e.,  interactions  between  the  gas  bubbles 
and  the  liquid  metal.  To  date,  however,  all  reported  experiments  have 
measured  only  gross  or  terminal  parameters  (voltage,  current,  pressure) 
or  averages  over  (part  of)  the  channel  (pressure,  void  fraction).  To 
obtain  a better  understanding  of  the  local  processes,  the  generator  experi- 
ments have  been  expanded  to  include  some  local  measurments.  Both  hot-film 
and  resistivity  probes  were  tested  in  the  previous  channel,  LT-3.  The 
initial  results,  reported  in  Section  IV,  give  evidence  that  large  electro- 
magnetic force  densities  have  a significant  effect  on  the  size  and  distribu- 
tion of  the  gas  bubbles  or  voids. 


1 1. 5 Analytical  Studies 
II. 5.1  End  Effects 

End  effects  in  LMMHD  generators  can  be  a significant  loss  mechanism, 
particularly  because  the  coupling  between  increased  terminal  voltage  and 
increased  end  loss  limits  the  attainable  voltage  in  high-efficiency  genera- 
tors. A computer  code  has  been  developed  that  will  calculate  the  electrical 
end  loss  (fringing  currents)  for  any  insulating  vane  geometry  and  magnetic 
flux  density  distribution,  a significant  accomplishment.  The  development 
of  the  code  and  some  initial  results  are  presented  in  Section  V,  and  the 
code  is  described  in  Appendix  C. 


* Supported  in  part  by  the  U.S.  Department'  of  Energy 


II. 5. 2 Imperfect  Compensation 

Unless  compensation  is  provided  for  the  induced  magnetic  field 
caused  by  the  currents  in  the  conducting  fluid,  the  total  magnetic  field 
varies  exponentially  along  the  channel  and,  for  a generator,  most  of 
the  energy  conversion  occurs  in  the  downstream  end  of  the  channel.  The 
result  is  a reduced  generator  efficiency.  (A  similar  phenomenon,  called 
armature  reaction,  occurs  in  dc  rotating  machines.)  To  eliminate  the 
induced  field  or  compensate  the  generator,  the  fluid  current  is  returned 
through  the  air  gap  of  the  electromagnet.  If  the  compensation  is  perfect, 
i.e.,  the  return  current  is  at  the  same  axial  location  as  the  fluid  current, 
the  induced  field  lies  along  the  flow  and  has  no  effect  on  the  generator. 
However,  because  of  both  the  end  currents  and  the  varying  fluid  properties 
and  velocity  along  the  flow,  the  compensation  will  be  imperfect,  as  has 
been  shown  experimentally  [see  Petrick  et  al  (1976)].  The  objective 
of  this  analytical  study  is  to  model  the  distortion  in  the  magnetic  field 
caused  by  imperfect  compensation  and  assess  the  consequent  effect  on 
generator  performance. 

The  study  has  been  divided  into  a one-dimensional  model  and  a two-dimen 
sional  model.  In  the  one-dimensional  model,  each  vector  quantity  has  only 
one  component,  which  is  a function  only  of  axial  position.  In  the  two-dimen 
sional  model  the  magnetic  flux  density  and  magnetic  intensity  have  two 
components  that  are  functions  of  two  variables,  whereas  all  other  vector 
quantities  have  one  component  that  is  a function  of  one  variable.  The  one- 
dimensional model  is  restricted  to  inf initely-permeable  iron  and  is  used 
to  indicate  general  trends,  whereas  the  two-dimensional  model  considers 
real-iron  effects. 
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The  solution  of  the  nonideal  compensation  problem  can  be  reached 
to  the  solution  of  two  subproblems,  1)  the  determination  of  the  current 
distribution  in  the  air  gap  caused  by  a given  magnetic  flux  density 
distribution,  and  2)  the  determination  of  the  magnetic  flux  density  produced 
by  a given  air  gap  current  distribution. 

The  whole  solution  is  achieved  by  iterating  between  two  subproblem 
solutions.  To  start,  the  applied  magnetic  field  is  assumed  to  be  the 
total  field  and  the  current  distribution  for  that  field  is  calculated 
(using  the  circuit  model  described  below).  The  magnetic  field  produced 
by  the  current  distribution  is  then  calculated  and  compared  to  the  field 
used  to  calculate  the  current  distribution.  If  the  difference  is  not 
within  acceptable  tolerances  for  a solution  (specified  input),  the  calculated 
magnetic  field  is  used  to  determine  a new  current  and  field  distribution. 
This  new  field  is  compared  to  the  field  calculated  in  the  previous  iteration, 
and  the  iteration  process  is  continued  until  the  field  converges  to  a 

solution. 

The  subproblem  of  the  current  produced  by  a given  magnetic  field 

is  best  considered  by  using  a circuit  model  to  include  the  effects  of 

the  external  circuit  connections.  Both  the  one-  and  two-dimensional 
models  use  the  same  circuit  model.  The  channel  is  divided  into  a number 
of  sections,  and  the  magnetic  field  and  current  density  are  assumed 
constant  over  each  section.  The  compensating  bar  is  also  divided  into 
a number  of  sections,  not  necessarily  the  same  number  as  in  the  channel. 
The  end  resistances  are  estimated  from  a theoretical  model,  assuming 
that  a specified  percentage  of  the  end  currents  from  the  model  are  eliminated 
by  insulating  vanes  or  field  overhang.  The  electrodes  can  be  modeled 
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as  real  conductors  by  considering  the  resistance  in  the  electrodes  between 
the  channel  sections.  An  example  of  the  circuit  model  for  five  channel 
sections  and  five  compensating  bar  sections  is  shown  in  Fig.  II. 4. 
The  known  magnetic  flux  density  determines  the  generated  voltages,  and 
the  circuit  model  is  solved  for  the  currents  in  the  channel  and  compensating 
bar  sections,  using  the  fundamental  tie-set  matrix,  a standard  method 
of  circuit  analysis. 

The  second  subproblem,  the  determination  of  the  magnetic  flux  density 
produced  by  a given  air-gap  current  distribution,  is  handled  differently 
in  the  two  models.  The  magnetic  field  is  determined  directly  from  Ampere's 
law  in  the  one-dimensional  model  but,  in  the  two-dimensional  model  the 
magnetic  vector  potential  is  calculated  by  using  the  magnet  design  code, 
TRIM,  which  also  considers  nonlinear  iron  relationships. 

Codes  for  the  one-dimensional  and  two-dimensional  models  have  been 
written  in  PL/1  and  have  been  debugged.  These  codes  are  now  being  used 
to  run  parameter  sweeps. 
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III.  GENERArOR  PERFORMANCE 

The  data  presented  here  includes  all  that  obtained  with  the  LT-4 
channel.  In  contrast  to  the  recent  data  reported  for  LT-3  [see  Petrick  et 
al  {1978)3,  the  LT-4  channel  allowed  operation  at  high  velocity  and  high 
pressure  (high  magnetic  flux  density)  and,  thus,  generator  operation  was 
extended  to  sign  if icantly-h igher  power  densities.  The  emphasis  in  the 
following  data  is  on  high-power-density  operation,  but  comparisons  with 
previous  IT -3  low-power  data  are  given  as  appropriate.  The  measured 
efficiencies  and  power  densities  are  discussed  in  Sections  III.l  and  III, 2, 
respectively.  The  slip  ratio  data,  pressure  distributions,  and  voltage 
distributions  are  presented  in  Sections  III. 3,  III. 4,  and  III. 5,  respectively. 
A summary  of  all  of  the  experimental  data  for  the  period  October,  1977  to 
December,  1978  is  given  in  Appendix  A. 

III.l  Efficiency 

The  generator  isentropic  efficiency  is  plotted  in  Fig.  III.l  for  all 
LT-4  data  with  Rl  = 0.4  mn , B = 1.2  T,  and  a separation  tank  pressure  of 
Pst  = 0.14  MPa  gauge  (20  psig).  Measurements  were  made  for  liquid  flow 
rates,  ill  , of  3,  6,  9,  and  12  kg/s.  The  emphasis  was  on  obtaining  high 

At 

average  void  fractons,  a,  and  the  data  were  very  encouraging,  because  the 
efficiency  remained  high  at  the  high  void  fractions.  (High  generator 
average  void  fractions  are  necessary  for  good  system  performance  or  equivalent 
turbine  efficiency,  see  Section  II, 2.)  The  efficiency  increases  sharply 
with  increasing  liquid  flow  rate  (velocity),  as  anticipated,  because  of  the 
reduced  slip  ratio  (ratio  of  gas  velocity  to  liquid  velocity).  There  is 
only  a single  point  at  12  kg/s  because  some  experimental  difficulties 
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occurred  with  the  resulting  high  gas  flow  rate  {above  the  facility  design 
value).  Some  of  the  lower-efficiency  data  points  at  3 and  6 kg/s  were 
obtained  with  very  high  gas  flow  rates,  and  this  probably  resulted  in 
transitions  to  . er.iiannular  flow  with  high  slip  ratios  (see  Section  III. 3 
and  Appendix  A) 

The  generator  isentropic  efficiencies  measured  with  LT-4  are  compared 
with  similar  data  for  LT-2  and  LT-3  in  Fig.  III. 2,  all  measured  for  Rl  = 

0.25  mo,  B = 1.2  T,  and  pst  = 0.14  MPa  gauge.  Note  that  there  is  no  data 
for  greater  than  6 kg/s  for  LT-2  and  LT-3  at  B = 1.2  T,  because  the 
older  channels  were  not  designed  for  the  higher  pressures  required.  The 
efficiencies  for  LT-3  and  LT-4  at  iti  = 6 kg/s  and  high  void  fractions  are 
in  very  good  agreement.  This  is  ir.  spite  of  the  different  taper  for  the  two 
channels  and  the  fact  that  LT-4  is  operating  far  off  design  conditions, 
and  indicates  that  efficiency  is  not  strongly  sensitive  to  off-design  opera- 
tion. (At  a - 0 the  efficiency  for  LT-4  is  significantly  lower,  because  the 
larger  taper  or  cross  sectional  area  change  results  in  increased  circulating 
currents  and  losses,  and,  for  the  LT-3  case,  there  was  flow  separation, 
which  significantly  decreased  the  circulating  currents.) 

The  higher  liquid  flow  rates,  iti^,  for  LT-4  resulted  in  substantial 
increases  in  efficiency  at  the  higher  void  fractions,  see  Fig.  III. 2,  and 
efficiencies  in  excess  of  0.60  were  measured  for  a greater  than  0.5. 

This  justifies  the  effort  in  extending  the  experiments  to  higher  liquid 
flow  rates  (velocities)  and  pressures.  To  place  these  efficiencies  in 
perspective,  they  were  obtained  with  a small  generator  (approximately  20 
kWe)  that  is  basically  wall  limited,  used  ordinary  NaK  without  any  surfactive 
agent  (see  Section  IV),  and  had  no  provision  for  insulating  vanes  or  a 
special  field  profile  to  minimize  end  loss  (see  Section  V).  Thus,  these 
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Figure  II 1.2  Generator  efficiency  for  LT-2, 
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results  present  very  encouraging  evidence  that  large  LMMHD  generators  can 
be  built  with  efficiencies  of  0.8  or  higher. 


III. 2 Power  Density 

Power  density,  i.e.,  electrical  power  output  per  unit  volume,  in  LMMHD 
generators  is  important  because  higher  power  densities  mean  smaller 
generators,  which  are  cheaper  to  build  and  easier  to  integrate  into  the 
power  plant,  smaller  and  less  expensive  magnets,  and  reduced  heat  losses. 

In  testing  small  LMMHD  generators,  it  is  important  to  reach  power  densities 
typical  of  those  anticipated  for  large,  practical  generators,  so  that  a 
comparable  degree  of  electromagnetic  interaction  is  obtained,  thus  minimizing 
the  chances  of  encountering  unanticipated  problems  in  scaling  to  larger 
generators. 

The  electrical  power  output  and  power  density  for  the  Rp  = 

0.25  mn  data  are  shown  in  Fig.  1 1 1. 3.  The  power  density  (and  power)  are 
proportional  to  the  velocity  squared  and,  thus,  increase  rapidly  with 

increasing  ill  . The  maximum  power  output,  21.6  kWe,  corresponds  to  a 

A 

power  density  of  31.8  MWe/m^  at  an  average  liquid  velocity  of  22.4  m/s. 

This  is  probably  the  highest  power  density  obtained  to  date  for  any 
MHO  generator  that  operated  continuously  for  more  than  one  minute.  For 
comparison,  the  LMMHD  generator  design  for  the  ECAS  study  had  a power 
density  of  11.1  MWe/m^  at  a = 0.75  and  a velocity  of  30.5  m/s. 


III. 3 SI *p  Ratio 

As  described  in  the  previous  report  [Petrick  et  al  (1978)],  the  two- 
phase  slip  ratio  is  expected  to  decrease  with  increasing  magnetic  Froude 
nombe',  defined  as  tne  ratio  of  the  inertial  to  the  Lorentz  forces,  or 
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Here,  fiL  is  the  liquid  density,  ut  the  liquid  velocity,  0 the  electrical 
conductivity  of  the  two-phase  mixture,  B the  magnetic  flux  density,  F the 
load  factor  (ratio  of  terminal  voltage  to  generated  voltage),  and  the 
typical  bubble-length  scale.  From  Eq.  III.l,  Fr^  increases  as  u increases, 
and  the  resulting  decreases  in  slip  ratio  and  slip  loss  are  the  primary 
motivations  for  operating  the  LMMHD  generator  at  higher  liquid  flow  rates 
(velocities)  and  pressures.  (The  use  of  surfactants  and  foams,  Section  IV, 
will  have  a similar  effect  by  decreasing  the  bubble  size,  D^). 

The  average  slip  ratios,  R,  for  all  LT-4  data  are  shown  as  functions  of 

~ and  it)  in  Figs.  III. 4 and  III. 5 for  R|_  = 0.25  itiq  and  0.4  mQ,  respectively. 

A 

As  ill  (or  u ) is  increased,  R decreases  and  approaches  unity  (within  the 

i 1 

the  accuracy  of  the  experimental  data). 

The  variation  of  R with  rh  is  shown  in  Fig.  1 1 1.6  for  " = 0.54  and 

i " 

0.71.  As  ill  increases,  R decreases  very  sharply,  contributing  signifi- 

£ t 

cantly  to  the  increase  in  generator  efficiency  with  higher  flow  rates. 

The  curve  for  1 ~ 0.71  terminates  at  a lower  4i  value,  because  of  the 

i 

facility  design  limitations.  Note  that  similar  results  were  reported 
previously  [Petrick  et  al  (1978)  or  Fabris  et  al  (1977a)],  but  only  at 
B = 0.6  T,  a quarter  of  the  present  power  density. 


III. 4 Pressure  Distributions 

The  new  generator  channel,  LT-4,  was  provided  with  23  pressure  trans- 
ducers, rather  then  eight  pressure  gauges,  as  previously  to  provide 
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additional  information  on  the  pressure  distribution  inside  the  generator 
channel.  Seventeen  transducers  were  spaced  along  the  short  (electrical)  side 
of  the  channel,  to  provide  a detailed  picture  along  the  generator  length, 
including  the  end  regions,  and  six  transducers  were  located  in  the  center 
of  the  wide  (insulating)  side  of  the  channel  just  outside  of  the  magnet,  to 
provide  information  on  the  behavior  in  the  end  regions.  The  pressure 
behavior  in  the  end  regions  is  of  interest  because  it  is  anticipated  that 
the  end  currents  circulating  within  the  channel  will  cause  the  pressure  to 
be  higher  in  the  center  of  the  wide  side  than  along  the  electrodes  at  the 
generator  inlet,  and  this  results  in  a higher  velocity  close  to  the  electrodes, 
commonly  referred  to  as  M-shaped  velocity  profiles  [see  Branover  (1974)]. 

The  pressures  for  two  pure-liquid  runs  are  shown  in  Fig.  III. 7.  Note 
the  rapid  pressure  drop  in  the  first  portion  of  the  channel  and  the  pressure 
recovery  near  the  downstream  end;  this  is  characteristic  of  pure-liquid  MHD 
flows  in  a diverging  channel.  Circulating  currents  cause  power  generation 
near  the  inlet  and  pumping  near  the  outlet  for  either  loaded  or  open-circuit 
conditions.  The  pressure  in  the  center  of  the  wide  side  is  higher  than  at 
the  walls  at  the  generator  inlet  by  approximately  140  kPa  (20  psi),  a 
significant  amount,  so  that  M-shaped  velocity  profiles  would  be  anticipated. 
Note  that  the  use  of  insulating  vanes  to  suppress  end  loss  (circulating 
currents)  will  also  suppress  M-shaped  profiles. 

The  pressures  for  two  typical  two-phase  runs  are  shown  in  Fig.  III. 8. 

Note  the  relatively-smooth  pressure  decreases  along  the  electrode  length, 
and  the  pressure  difference  of  up  to  approximately  100  kPa  (15  psi)  near 
the  generator  entrance.  For  two-phase  flow,  this  pressure  difference  could 
cause  some  gas-liquid  separation  in  addition  to  M-shaped  profiles;  again 
these  would  be  suppressed  by  insulating  vanes. 
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DISTANCE  ALONG  GENERATOR  (m) 

» III. 7 Pure  liquid  pressure  profiles  along  channel 


DISTANCE  ALONG  GENERATOR  (m) 

Figure  III. 3 Two-phase  pressure  profiles  along  channel 


III. 5 Voltage  Profiles 


Three  sets  of  eight  voltage  probes  were  equally  spaced  between  the 
electrodes  to  determine  the  voltage  profiles  across  the  insulating  walls. 

The  probes  were  mounted  flush  with  the  wall,  and  were  located  0.0127, 

0.1937,  and  0.3620  m downstream  from  the  beginning  of  the  electrodes.  To 
interpret  the  resulting  voltage  profiles,  note  that  Ohm's  Law  can  be 
written  as 

J 

- = u B - E . (III. 2) 

a 1 

Because  the  load  factor  F = E/ujB  is  typically  about  0.85,  both  utB  and 
E are  large  compared  with  J/0,  and  the  shape  of  the  voltage  curve  between 
the  electrodes  is  determined  primarily  by  the  variation  of  ut,  rather 
than  the  variation  of  0 between  the  electrodes.  Thus,  the  spatial  variation 
of  u^  can  be  inferred  from  the  voltage  profiles. 

One  set  of  voltage  profiles  for  pure  liquid  flow  (Fig.  III. 9)  shows 
convincing  evidence  for  M-shaped  velocity  profiles,  as  expected  from  the 
pressure  data  for  this  case  (Fig.  III. 7).  The  magnitude  of  the  M-shaped 
variation  seems  to  increase  slightly  along  the  channel.  The  reason  for  the 
asymmetry  between  the  two  sides  is  not  known. 

The  voUage  profiles  for  one  two-phase  case  are  shown  in  Fig.  III. 10. 

It  is  harder  to  interpret  two-phase  voltage  profiles,  because  0 may  vary. 
However,  it  appears  that  the  voltage  gradient  is  reduced  close  to  the 
electrodes,  perhaps  because  ut  is  reduced.  (A  decrease  in  a,  i.e.,  a, 
would  have  a similar  effect.)  Adjacent  to  the  two  regions  of  reduced 
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Figure  1 1 1. 10  Two-phase  voltage  profiles  between  electrodes 


^*2**MwfrS?aF!jia«!s^^  <**?&y  ^ 


voltage  gradient  are  regions  of  high  voltage  gradient  these  would  cor- 
respond to  an  M-shaped  velocity  profile,  with  the  peaks  of  the  M closer  to 
the  electrodes  than  in  the  pure-liquid  case. 

One  important  result  of  the  new  pressure  distribution  and  voltage  pro- 
file data  is  to  emphasize  the  importance  of  the  end  regions  on  influencing 
the  local  behavior  of  the  generator  even  far  from  the  ends.  The  analytical 
studies  aimed  at  understanding  and  minimizing  end  loss  are  described  in 
Section  V. 
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IV.  LIQUID-METAL  FOAMS  AND  SURFACE -ACTIVE  AGENTS 


The  program  has  concentrated  on  bubble  tests  (Section  IV. 1)  and 
surface  tension  measurements  (Secfion  IV. 2)  because  the  former  indicate 
foamability  directly,  and  surface  tension  is  one  quantitative  measure  of 
surface  activity.  As  previous’y,  the  program  is  focussed  on  NaK  alloy  at 
ambient  temperature,  as  it  would  be  used  in  the  initial  foam  LMMHD  generator 
experiments. 

IV. 1 Bubble  Tests  in  NaK  Alloy 
IV. 1.1  Barium  as  an  Additive 

Early  in  the  fiscal  year,  it  was  discovered  that  barium  was  effective 
in  promoting  bubble  formation  in  NaK  77  (76.7  wt  % K nominal).  The  initial 
mixture  contained  0.45  wt  % (0.11  mole  %)  barium.  In  the  standard  test,  a 
single  orifice  in  a Pyrex  tube  was  used  for  gas  introduction;  full  surface 
coverage  by  bubbles  was  achieved  immediately  on  starting  the  gas  flow 
(helium  at  5 cc/min).  At  375  cc/min,  a column  of  bubbles  several  centimeters 
high  was  observed.  It  is  appropriate  to  describe  this  aggregation  of 
bubbles  as  a foam.  When  the  Ba-NaK  solution  was  transferred,  helium 
was  blown  under  the  surface  of  the  liquid  metal  with  the  transfer  pipet  in 
a manner  analogous  to  blowing  bubbles  in  a beverage  with  a soda  straw. 

Surface  coverage  of  the  liquid  metal  with  bubbles  in  the  5-cm  diameter 
vessel  was  achieved  quickly. 

When  nitrogen,  instead  of  helium,  was  used  as  the  bubble-blowing  gas, 
it  was  observed  that  bubble  performance  did  not  degrade  over  a two  day  test 
period  at  ambient  temperature.  In  this  test,  the  system  contained  0.22  wt  % 
Ba.  There  was  immediate  surface  coverage  by  bubbles.  Two  minutes  after 
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gas  flow  was  stopped,  one  bubble  still  survived,  and  broke  only  when  the 
apparatus  was  shaken  to  determine  whether  it  actually  was  a bubble. 

Barium  was  also  observed  to  act  as  a wetting  agent  in  NaK.  Stainless 
steel  or  Pyrex  surfaces  that  normally  do  not  become  wetted  on  immersion  in 
NaK  at  ambient  temperature  were  wetted  immediately  on  immersion  in  NaK 
containing  small  amounts  of  barium. 

Some  tests  were  made  of  bubble  performance  at  higher  temperatures.  A 
sample  of  NaK  77  with  0.43  wt  % barium  was  tested  at  493  K.  The  tube  was 
wetted  immediately  on  immersion,  and  the  bubbles  were  larger  than  those 
generated  at  ambient  temperature  and  seemed  to  break  faster.  Approximately 
one-half  the  surface  would  be  covered  with  bubbles,  which  persisted  only 
approximately  5 seconds  after  the  gas  flow  stopped.  For  comparison,  at 
ambient  temperature,  the  same  sample  showed  nearly  the  same  or  slightly 
greater  coverage  with  bubbles  lasting  approximately  7 seconds. 

Pure  NaK  77  at  448  K did  not  wet  the  tube  on  immersion  or  after  a 
three-minute  soak.  Only  approximately  20  % of  the  surface  was  covered,  by 
a few  large  bubbles  that  lasted  only  1 to  2 seconds.  Clearly,  barium 
enhances  bubble  stability  even  at  higher  temperatures,  although  there 
appears  to  be  some  decline  in  this  stability  at  the  higher  temperatures. 

IV. 1.2  Calcium  and  Magnesium  as  Additives 

NaK  77  containing  0.71  wt  % calcium  did  not  support  bubble  formation, 
nor  did  it  wet  Pyrex  or  stainless  steel  as  barium-NaK  solutions  were 
observed  to  do.  However,  when  barium  (0.75  wt  %)  was  added  to  the  mixture, 
the  excellent  bubbling  and  wetting  characteristics  were  again  observed. 

The  results  indicate  that  barium  is  effective  as  a bubble  promoter  in  the 
presence  of  calcium,  despite  the  fact  that  calcium  itself  is  ineffective  in 
promoting  bubbles. 


With  0.11  wt  % magnesium  in  NaK  77,  there  were  favorable  signs  of 
bubble  enhancement.  This  and  subsequent  bubble  tests  were  performed  in  the 
tensionmeter  pot,  the  vessel  for  surface  tension  measurements,  using  the 
same  stainless  steel  orifice  tube  used  earlier  in  the  Pyrex  bubble  test 
apparatus.  This  change  was  made  to  eliminate  any  question  of  the  Pyrex 
contributing  a bubble-promoting  surfactant.  Adding  0.21  wt  % barium  to  the 
magnesium-containing  mixture  conferred  bubble  stabilization  and  wettability 
much  the  same  as  seen  for  the  calcium-containing  system. 

IV. 1.3  NaK  46  (Sodium  as  an  Additive  to  NaK  77) 

In  bubble  tests  with  pure  NaK  46  (46.2  wt  % K),  it  was  noticed  that 
the  geometry  of  the  test  setup  can  influence  bubble  formation,  a result 
that  indicates  the  importance  of  mechanical  design  in  exploiting  foaminess 
in  a LMMHD  machine.  It  was  found  that  bubbles  lasting  2 to  4 seconds  would 
appear  at  the  liquid  surface  if  the  L-shaped  bubbler  tube  was  held  so  that 
bubbles  rose  up  away  from  the  main  shaft  of  the  tube.  No  bubbles  appeared 
at  the  surface  if  the  tube  was  held  in  the  inverted  position.  In  that  case 
bubbles  presumably  coalesced  on  the  shaft  on  impact.  Adding  0.072  wt  % 
barium  to  the  NaK  46  caused  a dramatic  improvement  in  bubble  performance. 
Bubbles  could  easily  be  blown  to  cover  the  2-inch  (6-cm)  diameter  surface. 
Bubbles  could  be  formed  regardless  of  the  position  of  the  orifice  tube,  and 
it  took  35  to  45  seconds  for  all  of  them  to  break  after  the  gas  flow 
was  stopped.  A lower  concentration  of  barium,  0.020  wt  %,  was  not  effective 
in  improving  bubble  performance. 

IV. 1.4  NaK  85  (Potassium  as  an  Additive  to  NaK  77) 

With  pure  NaK  85  (85.15  wt  % K),  approximately  20  % of  the  surface  was 
covered  with  bubbles  that  lasted  3 to  5 seconds,  and  no  wetting  of  the  stain 
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less  steel  orifice  tube  was  noted  on  immersion  in  the  liquid  alloy.  With 
0.029  wt  % barium,  essentially  the  same  results  were  obtained.  With  0.16  wt  X 
barium,  the  orifice  tube  became  progressively  wetted  on  repeated  immersions. 
The  bubbles  covered  approximately  one-half  of  the  surface  and  lasted  3 to  5 
seconds.  Further,  bubbles  would  form  with  the  tube  in  any  position.  Overall, 
the  results  of  tests  of  Nak  46  and  NaK  85  show  better  performance  with  NaK 
46  and  an  apparent  threshold  of  effectiveness  for  the  barium. 

IV. 1.5  Silicon  as  an  Additive 

NaK  77  treated  with  high-purity  silicon  was  tested.  The  solubility  of 
silicon  appears  low,  so  the  system  was  probably  saturated.  On  immersion, 
the  orifice  tube  was  not  wetted.  The  performance  was  no  better  than  that 
of  NaK  77  alone;  approximately  25  % of  the  surface  would  be  covered  with 
bubbles  that  lasted  about  3 seconds.  When  0.080  wt  % barium  was  added  to 
the  Si/NaK  solution,  bubble  performance  and  wettability  tended  to  be 
improved. 

IV. 1.6  Lead  as  an  Additive 

NaK  77  with  0.083  wt  % lead  was  tested.  The  orifice  tube  was  not 
wetted,  even  after  repeated  immersion.  Approximately  one  third  of  the 
surface  would  be  covered  with  bubbles  lasting  3 to  5 seconds,  and  a single 
surviving  bubble  would  last  8 to  10  seconds.  The  performance  is  slightly 
better  than  that  with  NaK  77  alone.  Barium  (0.075  wt  %)  was  added  to  the 
Pb/NaK  solution.  The  orifice  tube  was  not  wetted,  and  approximately 
one-third  of  the  surface  became  covered  with  bubbles  lasting  some  3 seconds. 
With  an  increased  concentration  of  barium  (0.18  wt  X),  wetting  was  easily 
achieved,  and  approximately  one-half  of  the  surface  would  be  covered  with 
bubbles  lasting  about  5 seconds.  The  threshold  for  the  effectiveness  of 
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barium  in  this  case  is  influenced  by  the  fact  that  extensive  compound 
formation  is  possible  in  the  Ba-Pb  system. 

With  still  higher  concentrations  of  lead  and  at  higher  temperatures, 
the  following  results  were  obtained.  In  a sample  of  NaK  77-1.05  wt  % Pb, 
the  orifice  tube  did  not  become  wetted  on  immersion  or  on  swirling  at 
315  and  463  K.  At  both  temperatures,  bubbles  lasted  approximately  2 seconds; 
at  the  lower  temperature,  however,  somewhat  more  surface  coverage--as  much 
as  half--could  be  obtained.  At  the  higher  temperature,  25  to  30  % coverage 
was  observed.  This  performance  does  not  appear  to  be  significantly  better 
than  is  found  with  NaK  77. 


IV. 1.7  Gallium  as  an  Additive 

Gallium,  probably  at  the  saturation  level,  was  found  to  confer  slightly 
better  bubble  performance  than  NaK  77  alone,  but  no  wetting  occurred.  With 
0.17  wt  % barium  also  present,  the  orifice  was  wetted  after  three  immersions; 
30  to  50%  of  the  surface  would  be  covered  with  bubbles  lasting  approximately 
5 seconds. 


IV. 1.8  Cadmium  as  an  Additive 

Cadmium,  which  was  reported  on  as  an  additive  last  year,  was  reexamined. 

A mixture  with  0.20  wt  % Cd  did  not  wet  the  orifice  tube,  and  no  bubbles 
could  be  formed  with  the  orifice  tube  in  the  inverted  position.  In  the 
upright  position,  up  to  4 bubbles  could  exist  for  4 to  5 seconds,  a perfor- 
mance slightly  better  than  that  of  pure  NaK  77  under  these  conditions.  This 
situation  was  essentially  unchanged  when  barium  was  added  to  bring  the 
barium  concentration  to  0.20  wt  %,  and  then  to  0.43  wt  %.  However,  at 
0.70  wt  % Ba,  the  bubbler  tube  was  wetted  immediately  and  approximately  75  % 
surface  coverage  with  bubbles  lasting  about  6 seconds  was  obtained.  A single 
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bubble  would  survive  approximately  12  seconds.  The  reason  for  the  ineffec- 
tiveness of  initial  barium  additions  is  probably  the  reaction  between 
cadmium  and  barium;  compounds  ranging  from  BaCd^  to  BagCd  are  known.  The 
final  stoichiometry  in  this  test  mixture  was  equivalent  to  Ba2.gCd,  i.e., 
enough  barium  was  included  to  exceed  the  cadmium  capacity  for  formation  of 
Ba2Cd  and,  thus,  to  be  effective  in  bubble  stabilization. 

A test  with  cadmium  at  higher  temperature  was  also  included.  NaK  77 
with  0.58  wt  % cadmium  at  333  K did  not  wet  the  bubble  test  orifice  tube, 
and  provided  approximately  25  % surface  coverage  by  bubbles  lasting  2 to  3 
seconds.  This  performance  was  not  as  good  as  that  obtained  with  barium. 


IV. 1.9  Bismuth  as  an  Additive 

In  a bubble  test  with  NaK  77  containing  1.15  wt  % bismuth,  the  orifice 
tube  was  not  wetted  on  immersion  or  swirling.  Bubbles  could  be  blown  to 
produce  approximately  50  % surface  coverage  lasting  5 to  7 seconds.  At 
higher  temperatures,  the  performance  decreased,  as  with  barium,  so  that 
only  25  % surface  coverage  could  be  achieved  at  453  K with  the  bubbles 
lasting  1 to  2 seconds.  Overall,  the  performance  was  about  as  good  as  with 
some  barium-containing  systems,  but  not  as  good  as  that  found  in  the  most 
favorable  cases  with  barium.  Further  work  should  be  done  with  bismuth  to 
confirm  and  extend  these  data,  particularly  if  nitrogen  reacting  with 
barium  at  ambient  temperature  becomes  a problem. 


IV. 1.10  Selenium  as  an  Additive 

Selenium  was  retested  under  the  new  conditions.  A saturated  sample 
did  not  wet  the  orifice  tube.  Approximately  30  to  50  % surface  coverage 
could  be  achieved,  with  bubbles  lasting  2 to  3 seconds.  This  is,  perhaps, 
slightly  better  performance  than  that  achieved  with  pure  NaK  77,  and  might 
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be  the  result  of  bubble  stabilization  by  undissolved  selenium  particulate. 

The  same  performance  prevailed  after  addition  of  0.17  wt  % barium.  With 
0.34  wt  % barium,  the  bubbler  tube  was  nearly  completely  wetted  on  immersion. 
Approximately  half  the  surface  was  covered  with  bubbles  lasting  5 to  7 
seconds,  and  the  bubbler  tube  was  effective  either  upright  or  inverted. 

IV. 1.11  Pure  Potassium  and  Potassium-Barium  Systems 

Studies  of  pure  potassium  and  potassium-barium  mixtures  were  included, 
to  provide  reference  data  for  the  NaK  77  system,  which  is  a potassium-rich 
system.  At  348  K (potassium  melts  at  337  K),  the  potassium  did  not  wet  the 
stainless  steel  bubbler  tube,  and  bubble  stability  was  not  observed.  This 
behavior  is  to  be  expected  from  a pure  material. 

With  mixtures  of  potassium  containing  0.11  wt  % barium  at  358  K,  and 
0.33  and  0.63  wt  % barium  at  353  K,  the  same  lack  of  wetting  and  absence 
of  bubble  stability  was  observed.  Barium  is  not  effective,  at  least  at 
these  concentrations,  in  promoting  bubble  stability  in  potassium;  this  is 
consistent  with  the  observation  of  the  elevation  of  surface  tension  by 
barium  in  potassium. 

The  sodium  concentration,  then,  is  an  important  variable  in  relating 
the  nonfoaming  K-Ba  systems  to  the  foaming  NaK  77-Ba  system.  A sample  of 
potassium  containing  1.19  wt  % sodium  and  0.62  wt  % barium,  held  at  353  K, 
was  found  to  support  an  occasional,  fairly  large  bubble  that  lasted  for 
a fraction  of  a second.  Although  there  was  no  wetting  of  the  orifice  tube, 
this  bubble  test  did  indicate  a tendency  for  increased  film  stability. 
However,  using  the  same  material  with  more  sodium  added,  to  increase  the 
composition  to  15.8  wt  % sodium  (and  0.53  wt  % barium),  resulted  in  a 
marked  change  in  solution  characteristics.  The  orifice  tube  was  wetted 
almost  completely  on  the  first  immersion.  As  regards  bubble  stabilization, 


approximately  one-third  of  the  surface  became  covered,  with  the  last 
surviving  bubble  remaining  for  approximately  one  minute.  At  a higher 
temperature,  498  K,  and  when  compared  with  other  experiments  at  this 
temperature,  some  improvement  in  performance  was  found.  The  orifice  tube 
was  wetted  after  a single  brief  immersion,  but  the  bubble  lifetime  was 
shorter.  Approximately  25  % of  the  surface  was  covered  with  bubbles,  with 
the  last  surviving  bubble  breaking  after  5 to  7 seconds. 

IV. 1.12  NaK  83  with  Barium  as  an  Additive 

The  composition  of  NaK  83  (82.96  wt  % potassium)  is  about  the  most 
potassium-rich  NaK  alloy  that  has  potential  use  at  ambient  temperature  in  a 
LMMHD  machine.  Bubble  performance  in  pure  NaK  83  indicated  only  slight 
bubble  persistence.  However,  when  the  system  contained  0.27  to  0.99  wt  % 
barium,  the  wetting  characteristics  of  the  solution  improved  markedly,  and 
bubbles  lasting  as  long  as  5 seconds  covered  20  to  35  % of  the  surface. 

For  a mixture  containing  1.28  wt  % barium,  a test  at  478  K showed  bubbles 
covering  50  % of  the  surface  and  lasting  for  approximately  5 seconds. 

IV. 1.13  Effect  of  Nitrogen  on  Bubble  Performance  with  Barium  as  an  Additive 
The  1.28  wt  % barium  system  mentioned  above  was  used  for  a test  of 
the  reactivity  of  the  dissolved  barium  with  nitrogen.  After  three  days  of 
bubbling  N2  through  the  sample  at  temperatures  that  were  increased 
stepwise  (323,  373,  423,  and  468  K),  a bubble  test  at  468  K showed  that  the 
orifice  tube  was  not  wetted  initially,  but  did  wet  after  swirling  during 
immersion.  Bubbles  could  be  blown  to  give  approximately  25  % surface 
coverage  and  lasted  for  about  1 second.  The  effect  of  the  barium  was 
definitely  diminished,  quite  possibly  because  of  reaction  with  nitrogen  at 
the  higher  temperature.  However,  when  the  same  sample  was  cooled  to 


-44- 


ambient  temperature  and  retested,  35  to  50  % of  the  surface  was  covered 
with  bubbles  that  lasted  5 to  7 seconds,  with  lone  survivors  lasting 
for  10  to  15  seconds.  It  appears  that,  even  if  considerable  barium  was 
removed  from  solution  by  reaction  with  nitrogen  at  high  temperatures, 
sufficient  barium  remained  for  bubble  stabilization  to  be  appreciable  at 
ambient  temperature. 

IV. 2 Surface  Tension  Measurements 

During  the  year,  apparatus  was  installed  to  measure  surface  tension  as 
a function  of  temperature.  Data  obtained  with  this  equipment  for  a particu- 
lar system  are  more  extensive  and  are  regarded  as  more  reliable  than  the 
earlier  measurements,  in  part  because  of  their  greater  range  and  also 
because  of  the  accumulated  experience  with  the  technique.  Linear  regression 
equations  expressing  the  surface  tension-temperature  data  are  tabulated  in 
Table  IV. 1,  and  curves  for  these  data  appear,  where  appropriate,  in  the 
following  discussion.  Table  IV. 2 contains  a complete  tabulation  of  surface 
tension  data  for  temperatures  at  and  above  ambient. 

IV. 2.1  The  Na-K-Ba  System  and  the  Rationale  for  Bubble  Stabilization 
by  Surfactants 

The  most  noteworthy  observation  gleaned  from  our  studies  on  surface 
tension  in  the  NaK  77-Ba  systems  was  that,  for  systems  showing  enhanced 
bubble  stability,  the  surface  tension  was  not  depressed  but  was,  in  fact, 
higher  than  that  of  pure  NaK  77.  Usually,  foamability  and  film  stability 
are  correlated  with  lower  values  of  surface  tension.  Actually,  it  is  more 
accurate  to  state  that  the  phenomenological  correlation  is  between  film 
stability  and  surface  activity  of  a solute.  The  question,  then,  is  one  of 
reconciling  surface  activity  with  increased  surface  tension. 
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Linear  Regression  Equations  for  Surface  Tension  Data 


y (inN/m)  = a - bT(C°)  with  95%  Confidence  Intervals  Shown 


System 

K 

K + 0.33  wt  % Ba 
K + 0.63  wt  % Ba 
K + 0.62  wt  % Ba  + 1.19  wt  % Na 
K + 0.53  wt  % Ba  + 15.8  wt  % Na 
NaK77.17  (Nav77)** 

NaK77  + 0.15  wt  % Ba 
NaK77  + 0.43  wt.%  Ba 
NaK77  + 0.58  wt  % Cd 
NaK77  + 1.15  wt  % Bi 
NaK77  + 1.05  wt  % Pb 
NaK  82.96  (NaK83) 

NaK83  + 0.27  wt  % Ba 
NaK83  + 0.52  wt  % Ba 
NaK83  + 0.75  wt  % Ba 
NaK83  + 0.99  wt  % Ba 

NaK83  + 1.28  wt  % Ba 


0.0717 

0.0777 

0.0781 

0.0724 

0.0679 

0.0690 

0.0518 

0.1773 

0.0618 

0.1276 

0.0583 

0.0601 

0.0569 

0.0396 

0.0667 

0.0599 

0.0195 

0.0623 


t 0.0069 
t 0.0093 
+ 0.0045 
t 0.0234 
1 0.0134 
t 0.0069 
t 0.0043 
t 0.0258 
t 0.0042 
+ 0.0099 
+ 0.0049 

t 0.0062 

+ 0.0091 
+ 0.0084 
+ 0.0061 

t 0.0070  (high  temp.) 
t 0.0169  (low  temp.) 
t 0.0099 


111.9  t 0.3 

113.2  t 0.5 

121.7  + 0.3 
123.4  + 1.2 

125.7  + 0.7 

120.3  + 0.4 

110.4  + 0.2 

169.4  + 1.6 

117.2  t 0.2 

148.7  t 1.4 

102.2  + 0.3 

116.4  1 0.4 

125.6  + 0.6 

125.4  + 0.5 

140.7  + 0.3 

140.4  + 0.3 

136.7  t 0.4 

138.5  t 0.5 


Temperature  range  25-200°C  (298-473  K) 

*fhe  number  in  a designation  such  as  NaK77  refers  to  the  mass  percentage  of 
potassium. 


Table  IV. 2 

Surface  Tension  Measured  by  Maximum  Bubble  Pressure  Method  at  Ambient  T emperature 


Identification 

Number 


Description 


Surface  Tension 
niN/m 


1 

NaK77  previously  treated  with  butanol.  Orifices 
freshly  cleaned. 

113.7 

2 

Fresh  supply  of  NaK77.  No  additives. 

114.8 

3 

NaK77  treated  with  dicyclopentadiene  at  323-333  K. 

117.9 

4 

Same  system  after  additional  treatment  with 
dicyclopentadiene  at  -v388  K and  after  standing 
overnight. 

114.3 

5 

Fresh  NaK  sample. 

114.0 

6 

Same  sample  with  0.22  wt  % Hg. 

113.0 

7 

Same  system  one  day  later  and  after  developing 
hot-wetting  technique  for  orifices. 

114.4 

8 

Same  system  with  1.18  wt  % Hg. 

114.3 

9 

Same  system  with  6.38  wt  % Hg. 

116.6 

10 

Fresh  NaK.  Helium  is  bubble-blowing  gas. 

114.7 

11 

Same  system,  but  using  nitrogen  to  blow  bubbles. 

113.6 

12 

Same  system  with  0.070  wt  % barium. 

He  is  bubble  gas. 

114.4 

13 

Same  system,  two  days  later. 

115.2 

14 

Same  system  with  0.45  wt  % Ba. 

118.3 

15 

Fresh  sample  of  NaK77  with  0.17  wt  % calcium. 

Helium  bubbles. 

116.0 

16 

Same  sample  seventeen  hours  later. 

115.8 
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Identification 

Number 

Table  IV. 2 (continued) 

Description 

Surface  Tension 
mN/rn 

17 

Same  sample  after  heating. 

120.4 

18 

Same  sample  after  additional  heating. 

116.7 

19 

Same  sample  three  days  later. 

120.1 

20 

Additional  calcium  added  to  preceeding  sample 
to  make  concentration  0.45  wt  %. 

117.9 

21 

Same  sample  after  standing  overnight. 

118.4 

22 

Same  sample  after  exploratory  bubbling  test 
with  Pyrex  pi  pet. 

118.5 

23 

Same  sample  after  heating. 

118.3 

24 

Same  sample  after  standing  overnight. 

119.1 

25 

Additional  calcium  added  to  same  sample  to 
bring  calcium  content  to  0.71  wt  %. 

Helium  bubbles. 

118.2 

26 

Same  sample  after  standing  five  days. 

118.9 

27 

Same  sample  with  0.18  wt  % barium. 

118.8 

28 

Same  sample  with  additional  Ba  to  give 

0.51  wt  % barium. 

124.6 

29 

Same  sample  after  heating  for  47  minutes. 

124.1 

30 

Same  sample  with  additional  Ba  to  give 

0.75  wt  % barium-. 

122.1 

31 

Same  sample  after  bubbling  through  for  3 days. 

122.4 

32 

NaK77  sample  with  0.005  wt  % magnesium. 

121.7 

33 

Same  sample  with  additional  Mg: 0.029  wt  % Mg. 

137.4 

34 

Same  sample  after  reheating. 

130.3 

35 

Same  sample  after  bubbling  helium  through 
overnight. 

130.1 
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Identification 

Number 

Table  IV. 2 (continued) 

Description 

Surface  Tension 
mN/m 

36 

Same  sample  with  additional  Mg : 0 . 11  wt  % Mg. 

126.7 

37 

Same  sample  after  reheating. 

124.0 

38 

Same  sample  after  bubbling  helium  through 
for  three  days. 

124.1 

39 

Same  sample  after  bubble  tests  and  reheating. 

123.4 

40 

Same  sample  with  0.21  wt  % barium  (no  heating). 

123.2 

41 

Same  sample  after  heating. 

123.8 

42 

Same  sample  three  hours  later. 

124.0 

43 

Same  sample  after  bubbling  helium  through 
overnight. 

124.4 

44 

NaK  with  46.43  wt  % K,  designated  NaK46. 

126.8-127.5 

45 

Same  sample  with  0.020  wt  % Ba. 

125.0 

46 

Same  sample  two  weeks  later. 

127.9 

47 

Same  sample  with  0.072  wt  % Ba. 

127.0 

48 

NaK  with  85.15  wt  % K,  designated  NaK85. 

109.9 

49 

Same  sample  16  hours  later. 

109.5 

50 

Same  sample  with  0.029  wt  % Ba. 

120.1 

51 

Same  sample  after  bubble  test  and  cleaning 
tensionmeter. 

121.8,  118.0 

52 

Same  sample  with  0.16  wt  % Ba. 

118.0,  119.3 

53 

Pure  NaK77. 

116.0 

54 

Same  sample  with  0.068  wt  % selenium. 

116.5 

55 

Same  sample  with  0.017  wt  % barium. 

116.4 

. 56 

Same  sample  with  0.34  wt  % barium,  total. 

114.3 

-49- 


In  developing  an  understanding  of  the  behavior  of  surfaces  in  NaK 
77-Ba  alloy  systems,  it  is  important  to  recognize  that  NaK  77  is  a two-com- 
ponent system.  Bradhurst  and  Buchanan  (1961)  have  reported  a few  surface 
tension  values  for  the  sodium-potassium  system,  which  indicated  that 
potassium  is  surface  active  in  sodium.  According  to  their  data,  the 
general  nature  of  which  is  consistent  with  the  results  reported  here,  the 
surface  tension  decreases  slightly  as  potassium  concentration  increases  on 
the  potassium-rich  end  of  the  system.  NaK  77  is  potassium-rich  with  a 
potassium  atom  fraction  of  0.66.  The  slight  bubbliness  achievable  with 
pure  NaK  77  might  derive  from  the  surface  tension  gradient  still  present  at 
this  concentration. 

Adding  barium  to  the  NaK  alloy  results  in  a three-component  system, 
rather  than  the  simpler  two-component  system.  It  is  a potassium-rich 
mixture,  and  it  seems  more  reasonable  to  examine  the  phenomena  from  the 
viewpoint  of  a potassim  system,  rather  than  a sodium  system. 

The  surface  tension  data  in  Figure  IV. 1,  all  for  measurements  coming 
from  this  study,  show  that  the  dissolution  of  sodium  in  potassium  increases 
the  surface  tension  of  the  sodium-potassium  mixture.  The  NaK  83  system  was 
chosen  because  it  is  as  close  as  is  possible  to  the  potassium-rich  end  of 
the  binary,  while  still  maintaining  a liquid  at  ambient  temperature.  To 
give  a feel  for  the  precision,  data  points  are  shown  in  Fig.  IV. 1 . This 
precision  is  typical  of  most  cases  examined.  For  simplicity,  only  the 
linear  regression  lines  are  shown  in  other  figures. 

Figure  IV. 2 shows  that,  for  concentrations  in  the  range  of  interest, 
barium  raises  the  surface  tension  of  pure  potassium,  in  other  words,  that 
barium  is  surface  inactive.  From  the  data  given  in  Figs.  IV. 1 and  IV. 2, 
the  conclusion  may  be  drawn  that  both  sodium  and  barium  are  surface  inactive 
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SURFACE  TENSION  froN/ml 


In,  or  raise  the  surface  tension  of  potassium.  When  both  additives  are 
present,  the  situation  is  like  that  shown  in  Tig.  IV. 3,  where  all  the 
relevant  data  are  plotted.  Here  substantial  increases  in  surface  tension 
are  seen  with  increasing  additions  of  sodium  to  a "base"  system  containing 
0.63  wt  X Ba  in  K.  Furthermore,  the  system  with  15.8  wt  X sodium  is 
virtually  the  same  composition  as  NaK  83  in  respect  to  the  sodium  and 
potassium  components.  Therefore,  the  uppermost  curve  implies  that  adding 
barium  to  NaK  83  will  also  raise  the  surface  tension. 

This  expectation  is  supported,  at  least  over  the  limited  concentration 
range  shown  in  Fig.  IV. 4.  The  surface  tension  of  NaK  83  (the  lowest  curve) 
is  taken  as  the  reference  case.  Generally,  the  surface  tension  increases 
on  adding  barium  to  NaK  83,  at  least  up  to  the  0.99  wt  X Ba  composition. 
However,  surface  tension  data  for  1.28  wt  X barium  in  NaK  83  parallel,  and 
are  slightly  below  the  curve  for  0.99  wt  X Ba,  implying  that  a maximum 
surface  tension  was  reached  at  0.99  wt  X Ba.  The  individual  data  are  not 
shown  in  the  figure,  but  are  tabulated  in  Table  IV. 1.  Although  appearing 
to  indicate  a maximum  in  surface  tension  in  this  concentration  region, 
these  two  curves  are  close  enough  together  that  the  experimental  error 
associated  with  both  of  them  does  not  yet  allow  an  unequivocal  conclusion 
about  a maximum.  At  this  stage  of  understanding,  only  a tentative  statement 
about  a maximum  is  offered,  for  use  in  the  following  discussion. 

Figure  IV. 4 also  shows  values  for  the  slopes  of  these  curves.  For 
0.99  wt  X barium,  there  is  a break  in  the  curve  at  approximately  90°C  (363  K). 
Examining  all  the  data  above  90°C  (363  K),  one  can  see  that  the  slope  goes 
through  a minimum  for  the  solution  containing  0.52  wt  X barium.  Further, 
for  the  data  above  90°C  (363  K)  the  values  of  the  slopes  are  the  same 
within  a 95  X confidence  range.  On  the  low  temperature  side  [below  90°C 
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Figure  IV. 3 Surface  tension  vs  temperature  for  potassium- 
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Figure  IV.4  Surface  tension  vs  temperature  for  NaK  83-Barium  combinations 
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(363  K) 3 » the  1.28  wt  % barium  system  did  not  show  a break,  so  that  there 
are  two  minima  in  the  values  of  the  slopes  in  the  low  temperature  region--one 
for  the  0.52  wt  % barium  case,  and  one  for  the  0.99  wt  % barium  system. 
Assuming  that  the  slopes  of  the  surface  tension-temperature  curves  reflect 

surface  entropy  (S  = - ^ t where  S is  entropy,  y is  surface  tension,  and  T 

is  temperature),  it  appears  that  surface  entropy  goes  through  two  minima 
below  90°C  (363  K)  and  through  one  minimum  above  90°C  (363  K). 

A rationalization  for  this  occurrence  can  be  developed  on  the  assumption 
that  barium  is  surface  active  in  NaK  in  the  sense  of  positive  adsorption, 
rather  than  in  the  sense  of  depressing  the  surface  tension  of  the  mixture. 

How  surface  activity  coupled  with  increasing  surface  tension  can  come  about 
can  be  understood  in  terms  of  the  general  form  of  the  Gibbs  absorption 
isotherm 

dr  = - ^ r.  du.  , (IV. 1) 

where  ir  is  the  surface  tension,  is  the  surface  excess  of  component  i, 
and  y is  the  chemical  potential  of  component  i.  For  the  three  component 
Na-K-Ba  system  it  is 

dY  = -rNa  dMNa  - rK  - rBa  ^Ba*  (IV*2) 

On  adding  barium  to  a NaK  sample,  one  can  be  sure  of  a positive  sign  for 
dyBa*  To  the  extent  that  dilution  is  influential,  dy^a  and  dyj(  will 
be  negative.  Conceivably,  therefore,  a positive  sign  for  rga  (positive 
adsorption  of  barium)  could  be  coupled  with  positive  signs  for  rNa  or 
r|(.  However,  chemical  interactions  could  also  lead  to  positive  signs  for 
duNa  or  dM«,  instead  of  the  negative  signs  assumed  on  the  basis  of 
dilution.  In  the  absence  of  detailed  thermodynamic  information,  the  point 
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to  be  made  from  these  considerations  is  that  one  may  consider  the  possibi- 
lity of  positive  adsorption,  coupled  with  a rise  in  surface  tension,  when 
dealing  with  a three-component  system,  whereas  this  coupling  is  not  allowed 
thermodynamically  in  a two-component  system.  If  barium  is  positively 
adsorbed,  the  barium-enriched  surface  can  be  viewed  as  having  a composition 
markedly  different  from  the  composition  of  the  bulk,  with  behavior  different 
from  the  bulk  in  terms  of  the  phase  diagram  relationships. 

Unfortunately , a ternary  phase  diagram  for  the  Na-K-Ba  system  is  not 
available;  thus,  the  changes  in  surface  tension  must  be  correlated  from 
knowledge  of  the  compounds  identified  in  the  Na-Ba  and  K-Ba  binary  systems. 
Of  the  two,  only  the  Na-Ba  diagram  (Fig.  IV. 5)  has  been  worked  out  in 
sufficient  detail,  the  K-Ba  diagram  (not  shown)  is  rather  featureless  and 
without  strong  interactions.  The  significant  feature  of  the  Na-Ba  diagram 
is  that  Na  and  Ba  form  two  compounds,  NaBa  and  Na4Ba.  Furthermore,  the 
temperature  range  of  NaBa  stability  is  greater  than  that  of  Na4Ba. 

Although  it  is  not  claimed  that  compound  formation  occurs  in  the  surface 
film,  it  seems  reasonable  to  expect  that  the  attractive  interactions  that 
lead  to  such  compounds  in  the  bulk  could  also  persist,  and  could  lead  to 
ordering  in  the  surface  film.  In  addition,  these  attractive  interactions 
are  expected  to  have  an  elasticity-enhancing  effect,  favoring  film  and 
bubble  stability. 

It  is  of  interest  to  speculate  further  on  the  correlation  between  the 
information  given  in  the  Na-Ba  phase  diagram  and  the  slopes  of  the  surface 
tension  curves  for  the  NaK  83-Ba  mixtures.  At  low  temperature,  the  phase 
diagram  shows  that  two  compounds  or  two  kinds  of  interactions  are  possible, 
whereas  only  one  such  interaction  persists  at  higher  temperature.  This 
situation  might  be  regarded  as  being  reflected  in  the  two  minima  in  surface 
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entropy  below  90°C  (363  K),  and  in  the  single  minimum  above  90°C  (363  K). 

With  the  conceptual  framework  for  understanding  the  potassium-rich 
systems  developed  on  the  basis  of  NaK  83,  similar  reasoning  can  be  applied 
to  the  NaK  77  alloy  system  that  is  of  present  interest  in  the  LMMHD  program. 
Figure  IV. 6 shows  a series  of  surface  tension-temperature  curves  for  pure 
NaK77,  and,  for  reference,  pure  potassium  is  included.  For  the  NaK77 
alloy,  an  increase  in  surface  tension  is  noted  for  the  0.43  wt  % barium 
system.  This  increase  is  very  substantial,  as  is  the  steepness  of  the 
slope.  Both  features  need  further  checking,  but  the  increase  in  surface 
tension  is  qualitatively  consistent  with  data  collected  in  earlier  measure- 
ments of  several  NaK  77-Ba  alloy  mixtures  at  ambient  temperature.  However, 
an  additional  feature  of  interest  is  the  decrease  in  surface  tension 
observed  with  the  0.15  wt  % barium  mixture.  The  decrease  in  surface 
tension  values  at  low  barium  concentration  might  be  viewed  as  a vestige  of 
the  behavior  of  barium  in  pure  sodium,  as  reported  by  Addison  et  al  (1962). 
Some  of  the  Addison  data  are  shown  in  Fig.  IV. 7.  The  curve  for  150°C  (423  K) 
has  been  displaced  upward  by  50  mN/m  for  clarity.  The  minimum  in  the 
surface  tension  at  about  4 at  % is  at  a rather  high  concentration  on  a 
weight  basis,  approximately  20  wt  %.  It  seems  indicated  that,  in  comparing 
the  present  data  with  the  Addison  data,  the  minimum  in  surface  tension 
should  shift  to  lower  barium  concentrations  with  increase  in  potassium 
content. 

The  hypothesis  that  attractive  interactions  of  the  type  tending  to 
form  compounds  occur  in  the  surface  receives  additional  support  from  two 
types  of  calculations.  The  first  relates  to  the  variation  in  the  surface 
entropy  invoked  in  connection  with  the  minimum  values  of  the  slopes  of  the 
surface  tension-temperature  curves  of  the  NaK  83-Ba  systems.  On  the  assump- 
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Figure  IV. 7 Surface  tension  vs  atom  % barium  in  sodium  at  150  and  400°C 


tion  that  there  is  an  analogy  between  the  entropy  of  mixing  and  the  surface 
entropy  of  a mixture,  an  investigation  was  made  of  the  effect  that  the 
addition  of  a third  component  makes  on  the  entropy  of  mixing  in  the  binary 
system,  with  the  added  condition  that  one  or  two  compounds  can  form.  The 
calculation  uses  the  expression  for  the  entropy  of  mixing, 

S = - R Z.x.  In  x.  , (IV. 3) 

where  is  the  atom  fraction  of  species  i,  and  R is  the  gas  constant. 

It  is  further  assumed  that  one  of  the  binary  components  is  decreasing  in 
concentration  with  formation  of  the  compound  in  which  it  is  involved,  and 
concurrently  the  concentration  of  that  compound  is  increasing.  If  further 
addition  of  the  reactive  third  component  results  in  the  formation  of 
a second  compound  at  the  expense  of  the  first,  then  the  concentrations  of 
all  species  are  adjusted  accordingly.  The  result  of  such  an  approach  is 
that  the  entropy  of  mixing  could  be  made  to  go  through  either  one  or  two 
minima,  depending  on  whether  one  or  two  compounds  are  invoked.  Recalling 
that,  in  a given  temperature  regime,  one  or  two  minima  were  found  in  the 
slopes  of  the  surface  tension-temperature  curves  for  a series  of  compositions 
for  the  NaK  83-3a  mixture,  this  calculation  can  be  regarded  as  supporting 
the  idea  that  attractive  interactions  could  be  operating  in  the  surface  to 
stabilize  a preferred  atom  aggregate  composition. 

Again  using  the  data  for  the  NaK  83-Ba  system,  the  second  calculation 
follows  a rigorous  application  of  the  Gibbs  definition  of  a surface,  which 
sets  the  contribution  of  the  solvent  equal  to  zero  in  the  overall  calcula- 
tion. In  the  NaK  83  series,  the  barium  additions  were  in  small  increments, 
and  simultaneous  equations  for  the  Gibbs  adsorption  isotherm  were  set  up 
as 

dy  = -rNa  din  aua  - rBa  din  aga  . (IV. 4) 
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where  a is  the  activity,  and  the  other  terms  are  as  defined  previously. 

(The  term  for  potassium  does  not  appear  because  the  conventional  Gibbs 
definition  of  the  surface  is  being  followed.)  When  activities  were  taken 
equal  to  atom  fractions,  the  values  for  rNa  and  rBa  were  physically 
nonsensical,  values  of  the  order  of  -10"^  mole/cn^  were  obtained, 
whereas  10“^  mole/cm^  would  be  reasonable.  Further,  using  values  of 
df  that  were  at  the  extremes  of  the  95  % confidence  ranges  resulted  again 
in  nonsensical  values.  However,  more  physically  acceptable  values  for  r 
could  be  obtained  if  the  values  of  a^  were  reduced,  as  they  would  be  for 
nonideal  solutions  with  an  activity  coefficient  of  less  than  one.  Such 
nonideality  is  typical  of  systems  exhibiting,  or  showing  tendency  toward, 
compound  formation.  The  calculations  with  reduced  activity  are,  therefore, 
taken  to  be  supportive  of  the  idea  of  attractive  interactions  existing  in 
the  surface. 

As  for  the  more  practical  problem  of  stabilizing  a liquid  metal  foam, 
it  is  suggested  that  liquid  metal  films  cannot  have  available  one  of  the 
mechanisms  that  operates  to  stabilize  aqueous  and  organic  films,  namely, 
the  repulsive  effect  of  electrical  double  layers  on  opposite  sides  of  the 
film.  The  conductive,  metallic  medium,  in  effect,  short  circuits  such  a 
condition.  In  dealing  with  liquid  metals,  greater  reliance  must,  therefore, 
be  placed  on  the  film-stiffening  effects  of  attractive  interatomic  forces. 

On  the  basis  of  these  ideas,  systems  with  the  potential  for  such  film 
enhancing  effects  can  be  chosen  from  among  those  in  whch  a component  is 
positively  absorbed  and  for  which  there  is  evidence  for  compound  formation 
in  the  bulk. 


IV. 2. 2 Surface  Tensions  of  Systems  of  NaK  77  with  Additives  Other  Than 
Barium 

IV. 2. 2.1  Cadmium 

Figure  IV. 8 shows  the  results  from  surface  tension  measurements  on  a 
NaK  77-0.58  wt  % cadmium  alloy.  Cadmium  is  seen  to  lower  the  surface 
tension  of  NaK  77. 

IV. 2. 2. 2 Bismuth  as  an  Additive 

Data  for  the  NaK  77-1.15  wt  % Bi  system  show.  Fig.  IV. 9,  that  bismuth 
raises  the  surface  tension  of  NaK  77.  A linear  regression  line  for  the 
surface  tension  of  pure  NaK  77  is  shown  for  reference.  The  scatter  for  the 
points  as  shown  is  considerably  greater  than  in  prior  work,  and  might  be 
the  result  of  increased  tendencies  to  precipitate  alkali  metal  bismuth 
compounds  on  the  orifice  lips. 

A significant  feature  of  the  regresssion  line  is  that  the  slope  (a 
measure  of  surface  entropy)  is  approximately  double  that  of  the  curve  for 
pure  NaK  77,  i.e. , 0.123  mN/m  deg  vs  0.068  mN/m  deg.  It  may  prove  significant 
that  instances  of  increased  surface  entropy  have  been  observed  in  the 
systems  that  favor  film  stability,  i.e.,  Ba  and  Bi  in  NaK  77.  It  is 
noteworthy,  too,  that  bubble  stability  has  been  associated  in  these  studies 
with  increases,  rather  than  decreases,  in  surface  tension. 

IV. 2. 2. 3 Lead  as  an  Additive 

Surface  tension  data  for  the  NaK  77-1.05  wt  % Pb  system.  Fig.  IV, 10, 
indicate  significant  reductions  in  surface  tension.  This  fact,  together 
with  the  negligible  effect  of  lead  on  bubble  stability,  further  weakens  the 
conventional  correlation  between  foamability  and  decreased  surface  tension. 

If  the  lowered  surface  tension  is  an  indication  of  positive  absorption  of 
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Figure  IV. 9 Surface  tension  vs  temperature  for  NaK  77-Bismuth  combinations 
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Figure  IV. 10  Surface  tension  vs  temperature  for  NaK  77-lead  combinations 


lead,  it  might  be  possible  to  take  advantage  of  this  phenomenon  to  promote 
compound-like  surface  interactions  with  an  additional  solute,  to  stabilize 
films  and  bubbles. 

IV. 2. 3 Effect  of  Nitrogen/Barium  Interaction  on  Bubble  Stabilization 
Because  nitrogen  reacts  with  bulk  barium  to  form  Ba3N£  and  is 
known  to  react  with  barium  in  molten  sodium  at  high  temperatures  to  form 
Ba3N,  the  potential  for  a Ba-N2  reaction  in  NaK  should  be  evaluated  for 
cases  in  which  barium  is  used  as  a foam  promoter.  In  the  present  experi- 
mental program,  no  change  in  surface  tension  was  observed  for  a NaK  77-0.22 
wt  % Ba  solution  after  nitrogen  had  been  bubbled  through  it  for  two  days  at 
ambient  temperature.  A further  examination  of  the  reactivity  question  was 
done  in  a series  of  surface  tension  measurements  at  higher  temperatures 
using  a NaK  83-1.28  wt  % Ba  alloy  (as  in  Section  IV. 1.13)  Nitrogen  flowed 
continuously  through  tl  ' orifice  tubes  of  the  apparatus  and  measurements 
were  made  at  nominal  temperatures  of  323,  373,  423,  and  468  K.  The  results 
are  plotted  in  Fig.  IV. 11,  with  reference  curves  for  pure  NaK  83  and  NaK  83 
containing  0.27  and  1.28  wt  % barium.  The  data  indicate  that,  over  a 
period  of  about  3 days,  during  which  the  temperature  rose  stepwise  from  323 
to  423  K,  the  surface  tension  values  were  essentially  parallel  to  the 
surface  tension  curve  for  the  initial  liquid  metal  system.  The  near 
parallelism  tends  to  indicte  a relatively  slow  reaction  rate  below  423  K. 
However,  at  468  K,  the  surface  tension  departs  from  the  curve  noticeably, 
even  in  2 to  4 hours,  and  in  3 days  has  unmistakably  fallen  to  a value 
below  that  characteristic  of  a system  with  0.27  wt  % barium.  This  is  taken 
to  indicate  that  barium  has  been  removed  from  the  solution  by  combining  in 
a compound  with  nitrogen. 
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Figure  IV. 11  Surface  tension  of  NaK  83-1.28  wt%  Ba  after  exposui 
bubbling  nitrogen.  Indicated  times  are  for  length 
exposure  at  the  respective  temperatures,  times  are 
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For  the  perspective  of  utilization  of  a NaK-Ba  alloy  in  a LMMHD 
machine,  the  results  indicate  that  the  depletion  of  barium  by  reaction  with 
nitrogen  at  ambient  temperature  is  a slow  process.  The  fact  that,  even  at 
373  to  423  K,  the  rate  is  low,  together  with  the  fact  that  the  system,  which 
became  depleted  in  barium  at  468  K,  still  showed  good  bubble  characteristics 
when  retested  at  ambient  temperature,  suggests  that  good  bubble  performance 
at  ambient  temperature  could  be  expected  for  long  periods  of  operation. 
Exactly  how  long  remains  to  be  determined  by  further  laboratory /engineering 
evaluation  of  the  kinetics  of  the  Ba-N£  reaction  in  NaK. 


-70- 


V.  END  EFFECTS 


V. 1 . Introduction 

The  overall  efficiency  of  LMMHD  generators  can  be  significantly 
decreased  by  the  effects  associated  with  the  passage  of  the  working  fluid 
through  the  regions  at  the  upstream  and  downstream  edges  of  the  electrodes, 
where  a strong  gradient  of  the  magnetic  flux  density  exists,  unless  care  is 
taken  in  the  design  of  the  generator.  These  effects,  usually  called  "end 
losses",  include: 

a)  short-circuiting  electric  current  loops  closing  through  the  fluid 
regions  beyond  the  electrodes  (where  the  magnetic  flux  density  is 
lower),  i.e.,  the  "end  current  loss";  and 

b)  an  additional  pressure  drop  due  to  the  net  ponderomotive  force 
determined  by  the  above-mentioned  end  currents,  further  referred 
to  as  "the  end  pressure  loss."  This  is  in  excess  of  the  pressure 
drop  caused  by  the  load  current. 

The  effect  of  the  end  loss  on  the  generator  efficiency  can  be  better 
understood  by  the  aid  of  the  equivalent  electric  circuit  shown  in  Fig.  V.l. 
The  current  density  inside  the  channel  is  given  by  Ohm’s  law.  The  component 
perpendicular  to  the  electrodes  is 

0y=  U8);  (V.l) 

where  a is  the  electric  conductivity  of  the  working  fluid,  $ the  electric 
field  potential  created  by  the  voltage  difference  between  the  electrodes, 
and  u the  (constant  and  uniform)  fluid  velocity.  The  distribution  of  the 
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z component  of  the  flux  density  B is  represented  in  Fig.  V.l;  the  x and  y 
components  of  B are  taken  as  zero. 

In  the  region  between  the  electrodes,  the  current  density  has  the  sign 
of  the  "velocity- induced"  term  of  Eq.  V.l,  uB.  Beyond  the  electrodes, 
where  the  magnetic  field  strength  is  reduced,  the  "electrostatic"  term 
- 34 

jy  becomes  predominant  and  Jy  changes  sign,  giving  rise  to  a counter- 

current  or  end  current  Ie.  Thus,  beyond  a certain  curve  where  the 
current  changes  sign  the  fluid  acts  as  an  electric  shunt,  whose  resistance 
(the  end  resistance)  is  represented  by  Re.  Note  that  Re  includes  the 

effect  of  both  ends.  The  generator  internal  resistance  is  R ^ ; 

where  h is  the  electrode  half  spacing,  A the  aspect  ratio  (ratio  of  electrode 
length  2a  to  electrode  spacing  2h),  and  w the  channel  height  parallel  to 
B.  The  expression  for  R,  is  approximate  because,  depending  on  the  mag- 
netic flux  density  distribution,  some  of  the  generated  current  may  flow 
outside  of  the  region  between  the  electrodes. 

Various  methods  are  available  to  cont^o1  end  losses,  including: 

1.  Increasing  the  aspect  ratio  A.  This  leads  to  a decrease  in  the 
value  of  Rj  without  affecting  Re,  and  thus  proportionally 

less  end  loss.  However,  attainable  values  of  A are  limited  by 
the  available  pressure  difference  and  the  required  terminal 
voltage. 

2.  Shaping  the  decay  of  the  magnetic  flux  density  beyond  the  elec- 
trodes in  such  manner  as  to  increase  Re  (and  decrease  Rj. 

This  method  will  be  evaluated  in  Sections  V.4  through  V.6. 
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lev: 
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3.  Inserting  electrically  insulating  vane(s)  parallel  to  the  flow  in 
the  regions  outside  of  the  electrodes,  thereby  increasing  the  end 
current  path  length  and  Re  without  affecting  Ri  [Rossow  et  al 
(1961)].  This  method  implies  a certain  increase  in  the  viscous 
losses  in  the  end  regions;  the  increased  viscous  loss  has  to  be 
weighted  against  the  reduced  end  loss.  Vanes  are  considered 
throughout  the  remainder  of  Section  V. 

Several  aspects  of  the  above  methods  have  been  treated  in  the  literature 
by  analytic  two-dimensional  [Fishman  (1959),  Sutton  et  al  (1961),  Lavrentiev 
(1967,  1968),  Gubarev  et  al  (1970)]  and  three-dimensional  [Vasil iev 
and  Lavrentiev  (1970)]  approaches,  as  well  as  by  analog  modelling  [Moszynski 
(1967)].  A few  experimental  investigations  have  been  reported  [Moszynski 
(1967),  Petrick  et  al  (1968)].  The  effect  of  the  fringing  magnetic  field 
on  the  pressure  losses  has  also  been  estimated  [Hoffman  and  Carlson  (1971)]; 
however,  no  investigation  has  been  reported  on  the  effect  of  the  fringing 
field  upon  the  viscous  losses  associated  with  the  insertion  of  insulating 
vanes  in  the  flow.  Methods  for  the  reduction  of  the  end  losses  in  a 
two-phase  flow  have  also  been  proposed  [Petrick  et  al  (1968),  Moszynski 
and  Agrawal  (1968),  Branover  et  al  (1977)]  and  recently  tested  experi- 
mentally [Branover  et  al  (19/8)]. 

The  purpose  of  the  present  study  is  to  give  a quantitative  evaluation 
of  the  effectiveness  of  different  methods  for  end  loss  reduction  and  to 
supply  useful  data  for  the  design  of  the  LMMHO  generators.  The  model  is 
described  in  Section  V.2,  and  the  reference  case--no  magnetic  field  exten- 
sion past  the  electrode  ends--covered  in  Section  V.3.  Previous  analytical 
solutions  are  summarized  in  Section  V.4,  and  the  numerical  solution  is 
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described  in  Section  V.5.  The  numerical  results  are  compared  with  the 
analytical  solutions  and  analog  data  in  Section  V.6  with  zero  or  one 
insulating  vane.  The  results  with  insulating  vanes  are  presented  in 
Section  V.7  for  an  exponential  magnetic  field  extension  and  Section  V.8 
for  a natural  extension.  The  recommendations  for  the  reduction  of  end 
losses  are  in  Section  V.9. 

V.2  Model 

Like  most  of  the  previous  studies,  the  present  investigation  is 
based  on  a two-dimensional  approach  (ignoring  the  direction  parallel  to  B). 
The  channel  configuration,  magnetic  field  distributions  considered,  and  the 
notations  used  are  presented  in  Fig.  V.2. 

The  assumptions  made  are: 

a)  The  fluid  velocity  u is  uniform  and  constant  everywhere  in  the 
channel,  and  is  directed  along  the  x-axis.  Viscous  effects  are 
neglected. 

b)  The  magnetic  field  is  directed  along  the  z-axis  and  depends  only 
on  x.  The  magnetic  Reynolds  number  is  negligible  (compensated 
channel ). 

c)  The  electrodes  are  perfect  electrical  conductors;  the  walls  and 
vanes  are  perfect  insulators. 

d)  The  vanes  are  of  infinitesimal  thickness  and  do  not  disturb  the 
flow. 

V.3  The  Reference  Case:  No  Magnetic  Field  Extension 

A magnetic  field  of  constant  flux  density  along  the  electrodes  and  of 
zero  flux  density  beyond  them  was  chosen  as  the  "reference  case"  for 
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Figure  V.2  The  channel  configuration  and  magnetic  field  distributions 
considered:  l.-The  reference  case;  2. -Exponentially  de- 
caying overhang;  and  3. -Natural  overhang. 
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estimating  the  effectiveness  of  the  different  methods  chosen  for  the 
reduction  of  the  end  losses  because  of  its  simplicity  and  the  avail- 
ability of  established  analytical  solutions.  Note  that  it  is  neither 
realistic  nor  desirable  as  the  abrupt  termination  leads  to  large  end 
currents.  The  above  case  has  been  studied  both  analytically  [Sutton 
et  al  (1961),  Fishman  (1959)]  and  by  analog  modelling  [Moszynski  (1967)]. 
Also,  a preliminary  attempt  at  optimizing  the  number  of  insulating  vanes  to 
give  the  maximum  efficiency  is  available  for  this  case  [Moszynski  (1967)]. 

It  was  found  by  Sutton  et  al  (1961)  that  for  an  aspect  ratio  A > 0.7, 
the  resistance  ratio  Ri/Re  can  be  approximated  by 


Rj_  a 1 2 In  2 

Rc 


(V.2) 


\e  A IT 

In  this  case  the  maximum  attainable  efficiency  can  easily  be  related  to 
the  aspect  ratio,  as  plotted  in  Fig.  V.3.  They  also  determined  efficiency 
curves  as  a function  of  aspect  ratio  and  load  factor  F,  as  shown  in  Fig. 

V.4  for  A = 1,  2,  and  3.  (The  other  curves  in  Fig.  V.4  will  be  explained 
later,  and  all  curves  are  used  for  comparison  with  the  numerical  solutions.) 


V.3.1  The  Effect  of  Insulating  Vanes  on  the  End  Resistance 

If  it  is  assumed  that  an  "inf initely-long"  insulating  vane  at  the 
center  line  has  the  effect  of  "doubling"  the  apparent  aspect  ratio  of  the 
channel,  and  that  n inf initely-long  equidistant  insulating  vanes  multiply 
the  aspect  ratio  by  a factor  n+1,  then  the  effect  of  the  vanes  on  the 
maximum  attainable  efficiency  is  easy  to  infer  from  Eq.  V.2  and  Fig.  V.3. 
However,  there  is  no  analytical  or  experimental  confirmation  for  this. 

An  electric  analog  model  [Moszynski  (1967)]  seemed  to  indicate 
that  a vane  of  a length  7h  (2h  being  the  distance  between  the  electrodes) 
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Figure  V.3  The  maximum  attainable  efficiency  as  a function 
of  aspect  ratio  with  no  magnetic  field  extension. 
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Figure  V.4  Analytical  efficiency  curves  from  Sutton  et  al 
(1961)  for  aspect  ratios  1,  2,  and  3. 
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has  practically  the  same  effect  on  Re  as  an  infinitely-long  vane.  A 
later  analytical  solution  by  Lavrentiev  (1968)  found  that  to  attain  a 
shunt  resistance  Re  within  5%  of  that  for  an  infinitely-long  vane, 
the  finite  vane  should  have  a length  of  at  least  22h.  Also,  a vane 
7h  long  would  yield  a value  for  Re  15%  lower  than  for  an  infinite  vane. 

The  data  from  Moszynski's  (1967)  analog  experiment  may  also  be  used  to 
find  the  influence  of  the  number  and  length  of  the  insulating  vanes  on 
the  ratio  Ri/Re  and,  hence,  on  the  efficiency.  The  analog  data  seem 
to  be  well  correlated  by 

/Ri>  /Ri ' exp  [-Vn~  /t\0. 7 ] for  CKi<4h,  (V.3) 

^Rjn  " VRjnv  6-5  W 

e e 

as  shewn  in  Fig.  V.5,  where  n is  the  number  of  vanes,  i is  the  vane  length, 
and  the  index  nv  refers  to  the  case  with  no  vanes.  The  analog  data  do  not 
confirm  the  hypothesis  that  n infinitely-long  equidistant  vanes  multiply 
the  apparent  aspect  ratio  by  a factor  of  n+1,  as  shown  in  Fig.  V.5.  However, 
if  the  Re  values  are  corrected  by  15%,  as  suggested  by  the  study  of 
Lavrentiev  (1968),  an  infinitely-long  vane  would  indeed  have  the  effect  of 
doubling  the  apparent  aspect  ratio. 

V.3. 2 The  Effect  of  the  Insulating  Vane(s)  on  the  Viscous  End  Loss 

To  calculate  the  viscous  loss  caused  by  the  insertion  of  the  vanes, 
the  assumptions  needed  are  [Moszynski  (1967)]: 

a)  The  friction  factors  for  fully-developed  pipe  flow  apply, 

b)  The  vanes  are  thin  and  shaped  to  eliminate  contraction  and  enlarge- 
ment losses,  and 

c)  The  electromagnetic  interactions  have  a negligible  effect  on  the 
pressure  drop. 
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dependence  of  R^/Rg  on  the  number  and  length  of  the  insulating  vanes 


Experimental  measurements  of  the  pressure  drop  in  a varied  rectangular 
channel  [Moszynski  (1967)]  have  proven  that  the  first  two  assumptions  are 
valid  within  10*  accuracy.  The  third  assumption  is  clearly  true  for  the 
reference  case--no  magnetic  field  exists  beyond  the  electrodes  where  the 
vanes  are  situated.  In  an  actual  generator  the  velocity  profile  would  be 
modified  by  the  magnetic  field  and  this  assumption  might  not  be  valid. 

V.4  Analytical  Solutions  for  Magnetic  Field  Extensions 

The  cases  of  magnetic  field  overhang  examined  in  the  literature 
include:  a "constant"  overhang  (magnetic  field  decreasing  abruptly  to 
zero  at  a certain  point  beyond  the  electrodes),  a linearly  decaying  overhang, 
and  an  exponentially  decaying  overhang.  The  most  significant  reduction 
in  the  end  losses  was  found  by  Sutton  et  al  (1961)  to  be  associated  with 
the  exponential  overhang,  as  shown  in  Fig.  V.4.  The  characteristic  length 
of  the  exponential  field  overhang,  the  distance  from  the  point  where  the 
decay  begins  to  the  point  where  the  flux  density  has  decreased  to  1/e  of 
its  initial  value,  is  referred  to  as  "the  e-fold  length"  Le  (see  Fig.  V.2). 

The  results  presented  in  Fig,  V.4  suggest  that  an  increase  in  the 
e-fold  length  Le  brings  about  an  increase  in  efficiency-  This  dependency 
will  be  quantitatively  analyzed  in  Section  V.7.2.1.  Note  that  for  Le  = 

2h  the  maximum  efficiency  is  almost  twice  that  of  the  reference  case  (see 
also  Table  V.l). 

V.5  The  Numerical  Technique 

The  analytical  studies  referenced  above  determine  the  electric  field 
distribution  in  a longitudinal  section  of  the  MHD  channel  by  solving  the 
Laplace  equation  in  a domain  obtained  by  conformal  mapping  of  that  section 
or  of  a part  of  it.  The  insertion  of  an  insulating  vane  at  the  center- line 
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of  the  channel  makes  the  conformal  mapping  more  difficult,  but  the  case 
of  a semi-infinite  insulating  vane  at  center-line  was  solved  analytically 
by  Lavrentiev  (1967)  for  the  case  of  a magnetic  field  with  arbitrary  longi- 
tudinal distribution.  However,  no  analytic  solution  is  available  for  the 
cases  of  exponential  or  arbitrary  field  overhangs  and  one  or  more  vanes  of 
finite  length.  To  evaluate  the  reduction  in  end  losses  in  this  latter 
case,  a numerical  solution  was  sought  in  the  present  study,  as  described 
below. 

A computer  program  (ELLPACK)  for  solving  elliptical  partial  equations, 
presently  being  developed  in  a cooperative  effort  by  Harvard,  Texas  and 
Purdue  Universities,  was  adapted  for  finding  the  electric  field  distribution 
and  the  efficiency  of  the  LMMHD  generator.  In  an  early  stage  of  this  investi- 
gation, the  ELLPACK  program  [Rice  (1977)]  was  used  to  solve  the  Laplace  equa- 
tion in  a rectangular  domain  with  Dirichlet,  Neumann,  and  mixed-type  boundary 
conditions.  However,  the  solution  with  more  than  one  insulating  vane  present 
required  significant  modifications  in  the  program.  The  final  form  of  the 
program  is  described  in  Appendix  C. 

V.5.1  The  Mathematical  Model 

A schematic  drawing  of  the  channel  and  the  magnetic  flux  densities 
considered  is  given  in  Fig.  V.2.  Because  of  symmetry,  a complete  solution 
for  the  electric  field  can  be  obtained  by  solving  the  Laplace  equation  in 
one  quadrant  of  the  x-y  plane. 

The  mathematical  model  used  for  the  upper  right  quadrant  is: 
a2«>(x,y)/3x2  + 32*(x.y)/3y2  = 0,  0<x<xmax,  o<y<h  (V.4a) 

♦(x,h)  = FhuB0>  0<x<a  (V.4b) 

3*{x,h)/3y  = uB0b(x),  a<x<xmax  (V.4c) 
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♦ (xmax,j')  = 0, 

♦ ( x ,0 ) = 0, 

3*(0,y)/3x  = 0, 


0<y<h 

0<x<xmax 

0<y<h 


(V.4d) 

(V.4e) 

(V.4f) 


The  presence  of  an  insulating  vane  on  the  axis  between  the  points  (v,0)  and 
(ymax,0)  was  modelled  by  requiring  that  the  electric  current  density  normal 
to  the  vane  is  zero,  i.e.,  by  replacing  Eq.  V.4e  with: 

3*(x»0)/3y  = uB0b(x),  v_<x<vmax.  (V.4g) 

A convenient  feature  is  to  normalize  the  distance  between  the  electrodes 
as  2h  = 2 or  h = 1.  In  this  way,  the  x coordinate  of  the  electrode  end 
becomes  numerically  equal  to  the  aspect  ratio,  i.e.,  a = A.  The  infinite- 
length  channel  is  approximated  by  a finite-length  one,  with  end  points 
±xmax,  where  xmax  is  much  larger  than  either  a or  h. 

The  magnetic  flux  density  is  written  as  B = B0b(x),  where  B0  is  the 
flux  density  at  x = 0.  For  the  exponential  overhang  case  the  distribution 


function  is 

hfxl  J1'  °ixia 

^ ' lexp[-(x-a)/Le],  a<x  . 


(V.5) 


The  boundary  condition  Eq.  V.4d  requires  some  prior  knowledge  of 
the  .-lectric  potential  distribution,  namely  the  distance  from  the  electrodes 
where  the  electric  potential  can  be  considered  to  be  zero.  In  the  analog 
experiments  of  Moszynski  (1967)  it  was  found  that  the  region  between  x = lOh 
and  x = 20h  contributed  no  more  than  2%  to  the  end  current  shunt  resistance. 
This  was  used  in  choosing  xmax  = 16h. 
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V.5.2  The  Efficiency  Computation 

To  determine  the  net  current  which  can  be  drawn  from  the  system  (per 
unit  channel  width  in  the  2 direction),  Jy  has  to  be  integrated  along  a 
line  parallel  to  the  x-axis;  by  symmetry  this  can  be  2Jy  from  x » 0 
to  x = “ (to  x * xmax  in  the  program).  The  value  obtained  should  be 
independent  of  y for  y between  -h  and  +h.  The  total  electric  power  Pe  to 
the  load  per  unit  channel  width  is  the  product  of  this  current  and  the 
load  voltage 

VL  = 2F*0=  2FuB0h.  (V.6) 

The  result  is 

Pe=  4FuB0ho  j[  [-3$(x,y)/3y  + uB0b(x)]dx.  (V.7) 

The  generator  efficiency  is 


ni  = Pe/W’ 


(V  .8) 


where  W is  the  rate  of  flow  work  done  against  the  pressure  difference 
Ap  per  unit  width, 

W = u* Ap*2h* 1 • . (V.9) 

The  pressure  difference  Ap  is  found  by  integrating  the  ponderomotive  force 
«►  ♦ 

J x B along  the  channel  and  averaging  over  y,  i.e.. 


+h  +« 

Ap  = 1^  J*  / [Jy(x,y)B(x)]dx  = 
-h 


h • 

/ dy  Jdx[-34(x,y)/3y  + uB0b(x)]B0b(x) 
I 5 0 


(V  10) 
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It  can  be  shown  (see  Appendix  C)  that  the  group  ou^Bq?  is  a common 
factor  in  the  expressions  for  both  Pe  and  W,  and  that  the  terms  do  not  appear 
separately.  Therefore,  in  this  model  the  efficiency  n does  not  depend  on 
any  of  these  parameters. 

To  check  the  validity  of  the  current  calculations,  the  electric  current 
integration  was  performed  at  three  different  parallel  locations--y=h/4,  h/2, 
and  3h/4.  This  resulted  in  differences  of  up  to  8%  between  the  first  and 
last  level  locations.  However,  when  the  integration  was  performed  along 
the  path  MNP  in  Fig.  V. 2 this  difference  fell  below  1%.  An  alternative  way 
to  reduce  to  about  1%  the  differences  between  the  integration  results  is  to 
set  the  right  boundary  of  the  domain  at  xmax  = 32h.  Therefore,  for  all 
the  cases  with  no  more  than  one  insulating  vane  the  value  xmax  = 32h  is 
used. 

V.5.3  Plots  of  the  Electric  Potential  Inside  the  Channel 

The  solution  of  Eq.  V.4  was  provided  by  the  computer  both  as  tabulated 
values  of  the  electric  potential  4>(x,y)  for  every  mesh  point,  and/or  as 
plots  of  the  equipotential  lines  of  ♦.  The  latter  form  of  output  is  useful 
for  getting  a better  understanding  of  the  physical  effect  of  different 
magnetic  field  parameters,  insulating  vanes,  etc.  on  the  electric 
field  (and  implicitly  the  current  density)  distribution  inside  the  channel. 
Some  typical  plots  are  presented  in  Figs.  V.6  and  V.7. 

V . 6 C omparison  of  the  Nu me rical  Results  with  the  Existing  Analytica 1 

Solutions  and  Ana)og~TTafa 

A useful  test  of  the  computer  program's  accuracy  is  a comparison  of 
its  results  with  those  of  a few  of  the  existing  exact  solutions  resulting 
from  an  analytical  approach. 
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Figure  V.6  Electric  potential  distributions  for  different  load 
factors,  A = 1,  L /h  = 1/2,  and  no  insulating  vanes 
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V.6.1  The  Reference  Case 

The  maximum  efficiencies  reported  by  Sutton  et  al  (1961),  about 
0.29  and  0.40  for  A » 1 and  A = 2,  respectively,  as  well  as  the  efficiency  vs. 
load  factor  curves,  are  reproduced  b*  the  computer  code  with  an  accuracy 
within  IX.  This  is  demonstrated  by  comparing  Figs.  V.4.a  and  V.4.b  with 
Fig.  V.8.  Note  that  on  Fig.  V.8  and  the  following  figures  all  distances  are 
normalized  by  h,  the  half  distance  between  the  electrodes,  V and  Vmax  are 
the  positions  of  the  left  and  right  edges  of  the  insulating  vanes,  MAXEFF 
is  the  maximum  value  of  the  efficiency,  and  MAXLOAD  is  the  load  factor 
at  MAXEFF. 

V. 6.2  An  Exponential  Magnetic  Field  Overhang  with  No  Vanes 

The  dependence  of  the  efficiency  on  the  load  factor  has  been  computed 
by  Sutton  et  al  (1961)  on  the  basis  of  an  analytical  solution  for  the 
electric  potential,  see  Fig.  V.4.  The  computer  code  reproduced  the 
data  for  Le  = h and  2h  wu  less  than  3%  error  for  both  A = 1 and  A * 2 (see 
Figs.  V.9  and  V.10).  In  addition,  the  computer  results  for  Le  = h/2 
(one-quarter  of  the  interelectrode  distance)  are  shown. 

V.6.3  No  Magnetic  Field  Overhang  with  an  Insulating  .ane  on 
the  Center  Line 

The  electric  analog  of  this  configuration  reported  by  Moszynski  (1967)  was 
intended  to  determine  the  ratio  Rj/Re  for  a variety  of  vane  lengths  and 
positions  in  order  to  compute  the  generator  efficiency  according  to  the 
equation 

. j-Lp  0 - F(1  * Rj/Re)].  (V-U) 
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Moszynski  found  that  increasing  the  vane  length  i beyond  a value  of  7h  (3.5 
times  the  interelectrode  distance)  does  not  have  any  significant  effect  on 
the  Rj/Re  ratio;  and,  therefore,  a vane  of  this  length  would  be  practically 
equivalent  to  an  infinite-length  vane.  Using  the  data  for  R^/Re  reported 
by  Moszynski  (1967),  the  efficiency  vs.  load  factor  curves  for  the  cases  of 
no  vanes  and  center-line  vanes  of  lengths  4h  and  6h  are  plotted  in  Fig.  V . 11 . 

The  computer  code  accurately  (within  IX)  reproduced  the  curve  for  the 

"no  vane"  case  of  Moszynski  (1967),  Fig.  V.12.  However,  a significant 

difference  was  found  for  the  last  two  cases(*  = 4h  and  *•  = 6h).  This 

difference  is  apparently  due  to  the  fact  that  the  analog  experiment  was 

based  on  measuring  the  currents  set  up  by  the  electric  potential  between 

the  electrodes  and  could  not  simulate  the  "velocity  induced"  current 
♦ > 

due  to  u x B.  Moreover,  increasing  the  vane  length  beyond  * = 7h  in  the 
computer  code  brought  about  a further  increase  in  efficiency.  This  means 
that  the  asymptotic  limit  found  by  Moszynski  (1967)  was  not  confirmed  by 
the  numerical  technique.  This  agrees  with  the  theoretical  work  reported  by 
Lavrentiev  (1968). 

V . 7 An  Exponential  Magnetic  Field  Extension  With  Insulating  Vanes 
V.7.1  One  Inf initely-Long  Insulating  Vane  at  the  Center  Line 

When  the  magnetic  field  extends  exponentially  beyond  the  electrodes, 
the  insertion  of  insulating  vanes  on  the  center  line  from  the  points  x = ± a 
where  the  electrodes  end  up  to  the  boundaries  x = ±xmax  increases  the 
maximum  attainable  efficiency,  as  shown  in  Table  V..1  for  Le/h  = 1 and  2.  The 
efficiency  curves  for  aspect  ratios  of  1 and  2 and  for  vanes  with  *■  = 15h 
and  * * 14h  are  shown  in  Figs.  V . 13  and  V.14,  respectively. 


Efficiency 


Load  Factor 


Figure  V.ll  Efficiency  curves  computed  with  the  analog  ^ ta 
for  R.j/Re  from  Moszynski  (1976)  for  A = 2 
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TABLE  V.l 

Maximum  Attainable  Efficiencies  for  Exponential  Field  Overhang  and  One 

"Inf initely-Long"  Vane 

(The  percentages  in  the  brackets  indicate  the  increase  with  respect  to 

the  reference  case) 


As'  ‘ct 
ra^io 

Maximum  Efficiency 

Reference 

Case 

Exponential  overhang, 
no  vanes 

Exponential  overhang, 
one  "infinitely  long"  vane 

__ 

Le/h=I 

Le/h=2 

Le7n=I 

1 

.28 

.57  (203*)  " 

.65  (232*1) 

766  C235tf)~ 

.69  (246*)  “' 

2 

.40 

.65  (162%) 

717  (T80f) 

i 

.75  imrr 

.77  (192%) 

V.7.2  The  Influence  of  the  Magnetic  Field  Distribution  and  the 
Insulating  Vane  Length 
V.7.2.1  Magnetic  Field  Overhang  Length 

For  the  case  of  an  exponentially  decaying  magnetic  field  overhang 
it  is  obvious  that  an  increase  in  Le  leads  to  an  efficiency  increase,  see 
Figs  V.13  and  V.14.  While  this  increase  is  not  large  for  an  inf initely-long 
vane  on  the  center  line  (see  the  last  two  columns  of  Table  V.l),  better 
results  may  be  attained  by  optimizing  the  vane  length  and  position  (especially 
when  viscous  loss  is  included). 

V.7.2. 2 Vane  Length 

Since  it  is  obvious  that  an  infinitely-long  insulating  vane  would 
be  technologically  inconvenient,  it  seems  reasonable  to  choose  a length 
such  that  only  a small  efficiency  improvement  would  be  obtained  by  a 
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further  increase  in  length.  The  maximum  attainable  efficiency  for  A * 2 
and  Le  = h is  plotted  in  Fig.  V.15  against  the  insulating  vane  length. 

V. 7.2.3  Vane  position 

The  axial  position  of  a short  vane  has  been  found  to  influence 
significantly  the  efficiency.  The  efficiency  curves  of  Fig.  V.16  for 
different  axial  positions  of  a short  vane  of  constant  length  t = 8h  indicate 
that  the  best  efficiency  is  obtained  when  the  vane  begins  at  a distance 
equal  to  Le/2  outside  of  the  electrode  end.  This  was  confirmed  by 
computations  made  for  both  Le/h  = 1 and  2. 

V.7.3  The  Influence  of  Additional  Insulating  Vanes 

As  mentioned  in  Section  V.5.2,  the  computer  program  was  modified 
to  provide  the  solution  for  the  electric  field  with  insulating  vanes  at  the 
center  line  and/or  at  intermediate  positions  within  the  channel.  Since  the 
numerical  integration  technique  was  modified  for  the  new  version  of  the 
program,  tests  were  run  for  cases  with  an  available  analytical  solution  and 
for  the  case  of  one  insulating  vane  on  the  center  line.  The  results 
obtained  previously  (see  Sections  V.6.2  and  V.7.2)  were  reproduced  within 
the  accuracy  allowed  by  the  discretization  mesh  employed. 

To  run  the  new  program,  which  required  an  increased  number  of  horizontal 
mesh  lines,  within  a memory  field  length  of  150  K,  the  right  boundary  of 
the  domain  had  to  be  reduced  from  xmax  = 32h  to  xmax  = 16h,  then  to 
xmax  = This  change  did  not  affect  the  results  by  more  than  abou-t 

3%  (compare  the  maximum  efficiency  from  the  first  case  of  Fig.  V.16  with  that 
of  line  4 of  Table  5.2). 

V.7.3.1  Two  Insulating  Vanes 

For  a certain  vane  length,  the  best  efficiency  was  obtained  when 
the  two  insulatiing  vanes  were  approximately  equidistant  (within  the  accuracy 
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of  the  mesh  spacing),  Table  V.2  and  Fig.  V.17.  (With  two  vanes,  there  is 
no  vane  on  the  center  line.)  The  optimum  axial  position  for  the  edge  of 
the  vanes  was  found  to  be  between  x = a and  x = a + Le/2  (probably  around 
a + Le/4)  for  A = 2,  Fig.  V . 18 . With  two  insulating  vanes,  the  influence 
of  Le  is  relatively  weak,  as  shown  in  Fig.  V.19. 

The  efficiency  increase  aused  by  an  increasing  vane  length  does  not 
seem  to  reach  an  asymptotic  v._1-je  even  for  = lOh,  Figs.  V.20  and  V.21. 

Quantitatively,  the  improvement  in  efficiency  obtained  by  the  inser- 
tion of  two  vanes  of  a certain  length  with  respect  to  a single  vane  of  the 
same  length  was  found  to  be  about  4%,  as  shown  by  lines  4 and  12  of  Table 
V.2. 

V.7.3.2  Three  or  More  Insulating  Vanes 

The  increase  in  efficiency  attained  by  the  insertion  of  three  vanes, 
with  respect  to  the  two-vane  case,  is  only  about  1.5%,  Fig.  V.22.  The 
influence  of  Le  is  again  very  weak.  Fig.  V.23. 

With  five  vanes,  only  about  1.5%  increase  in  efficiency  is  obtained 
with  respect  to  the  case  of  three  vanes  of  the  same  length,  lines  15  and  16 
of  Table  V.2. 

V.8  The  Natural  Magnetic  Field  Overhang  with  Insulating  Vanes 

While  the  exponential  extension  of  the  magnetic  field  was  proven  to 
be  more  efficient  than  a "constant"  or  "linearly  decaying"  field  overhang, 
it  was  not  shown  how  it  compares  in  shape  or  efficiency  with  a real  magnetic 
field  distribution.  This  should  determine  the  magnetic  field  shape  that 
yields  the  best  efficiency. 
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Table  V.2 


Maximal  Efficiency  Attainable  with  Various  Insulating  Vanes 


7 2 

8 2 


10  2 


cc 

ui 

K-  Q£ 
ar  •-> 


*-i  o 
-J  UJ 

a. 


>•  < 
oaa. 


«/>  *-*  ui 
KQZ 

M LU  < 


O <-*  UI 

o o z 

UJ  UJ  < 
lOI> 


Mil  v/h  U-,  / h y-j/h  Vyh  |Jt2/  hjy2/h 


MAXIMUM 

ATTAINABLE 

EFFICIENCY 


Vh  J«"  >2 


8 2.5  j 

8 2.25 

8 0.5  2,25 

8 0.5  2.0 

8 0.5  2.25 

8 0.5  2.5 

8 0.37  2.25 

8 0.5  2.25 

8 0.63  2.0 

8 0.37  2.25 

' 8 0.37  2.5 

8 0.37  2.75 

8 2.25  8 0.5  2.0 


16  0.58 
16  0.687 
14  0.69 


14  0.73 


14  0.77 


8 2.25  8 0.33  2.25  8 0.66 


14  0.78 


EFFICIENCY 


Figure  V. 1 7 Efficiency  vs  load  factor  as  a function  of  v 
for  two  Insulating  vanes,  L /h  = 1 and  A = 2 
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Figure  V.18  Efficiency  vs  load  factor  as  a function  of  vane 

position  for  two  insulating  vanes,  L /h  = 1,  and  A 
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Figure  V.19  Efficiency  vs  load  factor  as  a function  of  L 

for  two  insulating  vanes  of  length  &/h  = 8,  A = 2 
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Figure  V.21  Efficiency  vs  the  length  of  two  insulating  equidi- 
stant vanes,  A ■ 2,  and  L /h  * 1 
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Figure  V. 22  Efficiency  curves  for  vanes  of  length  t/h  * 8 
A = 2,  and  LVh  = 1 
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The  longitudinal  distribution  of  the  magnetic  field  created  by  two 


plane-parallel  pole  pieces  can  be  described  with  a good  approximation 
by  [see  Laverentiev  (1968)] 


b(x)  « 


1 - Ai  exp[  vj ( Jx  I - x)/h]  for  lx  I < x 


A2exp[-^(  ix  1 - x)/h] 


for  lx  I > x , 


(V . 12 ) 


where  Aj  + A2  * 1, 

Am  = A2V2, 

and  2X  the  length  of  the  pole  piece.  This  magnetic  field  distribution 
is  shown  in  Fig.  V.24  for  values  of  **  * tyh  (the  total  air  gap  « between 
the  magnetic  poles  normalized  by  the  interelectrode  distance  h)  of  0.2... 1.4. 
The  slope  of  the  field  decay  is  determined  by  the  gap  *;  a smaller  gap 
leads  to  a more  abrupt  decay.  The  actual  values  used  in  Eg.  V.12  for 
Fig.  V.24  were  A^  = 0.17,  A2  = 0.83,  ^ a 5.29/**,  and  V2  c 1.07/**. 

It  should  also  be  noted  that,  if  the  length  of  the  pole  piece  is 
equal  to  the  electrode  length,  x « then  at  the  edge  of  the  electrode 
the  magnetic  field  has  decreased  to  83%  of  its  value  at  x * 0.  Since  the 
magnetic  field  is  not  constant  (and  equal  to  B0)  between  x « 0 and  x * a, 
the  voltage  between  the  electrodes  will  be  taken  as 


VL  = Fu2h 
a 


B0b(x)dx. 


(V. 13) 


In  Fig.  V.25  the  real  distribution  of  the  magnetic  field  for  **  * 0.5, 
1.0,  and  1.5  is  compared  with  the  exponential  decay  fur  Le  3 0.5,  1.0, 
and  1.5,  respectively.  This  comparison  will  be  useful  for  the  interpreta- 
tion of  the  following  results.  For  instance,  it  seems  evident  that  if  the 
magnet  edge  is  beyond  the  electrode  edge,  x>a,  the  voltage  between  the 
electrodes  given  by  Eq.  V.13  will  be  higher  for  the  same  B0,  improving 
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the  efficiency.  Quantitatively,  this  can  be  seen  in  Fig.  ¥.26,  where  the 
electrode  edge  has  the  position  a « 2,  i . e . , A * 2.  The  less-predictable 
result  is  that,  even  for  *>a,  the  efficiency  is  about  20%  lower  than  it  was 
for  the  case  of  an  exponential  field  overhang.  For  instance,  with  ■ 0.5 
and  x « 2.5  * 1.25a,  the  efficiency  is  about  0.49,  as  compared  with  0.58  for 
an  exponential  overhang  with  Le  « 0.5  (see  Table  V.2,  line  2).  The 
insertion  of  an  insulating  vane  at  the  center  line  improves  the  efficiency, 
Fig.  V.27,  but  the  latter  still  remains  below  the  value  attainable  with  an 
exponential  field  overhang,  line  3 of  Table  V.2. 

Again,  the  axial  position  of  the  vane  influences  the  maximum  attainable 
efficiency.  Fig.  V.27. 

V . 9 Recommendations  for  the  Reduction  of  End  Losses 

In  the  preceding  sections  three  ways  to  reduce  the  end  current  losses 
were  considered: 

a)  increasing  the  aspect  ratio  of  the  channel, 

b)  providing  a magnetic  field  distribution  approximating  an 
exponential  decay  beyond  the  electrodes,  and 

c)  inserting  insulating  vane(s)  in  the  regions  where  the  fluid 
enters  and  leaves  the  magnetic  field. 

The  choice  of  the  channel  aspect  ratio  is  primarily  determined  by 
the  rated  electric  power  and  the  allowable  pressure  drop.  However,  from 
the  point  of  view  of  the  end  losses,  a larger  aspect  ratio  leads  to  a 
better  efficiency.  The  data  obtained  in  the  present  study  indicate, 
curve  1 of  Fig.  V.28,  that  for  aspect  ratios  above  5 a maximum  efficiency 
higher  than  0.80  can  be  obtained  without  any  insulating  vanes,  provided  a 
proper  shape  of  the  magnetic  flux  density  distribution  is  achieved. 
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Figure  V.27  The  optimum  position  of  a center-line  vane 

edge  for  the  real  decay  of  the  magnetic  field 
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Once  a specific  aspect  ratio  is  chosen  for  the  channel,  the  most 
important  parameter  affecting  the  efficiency  was  found  to  be  the  magnetic 
field  distribution  along  the  longitudinal  axis.  If  the  field  is  created 
by  plane-parallel  poles  of  a length  equal  to  the  electrode  length,  its 
strength  decays  toward  the  ends  of  the  electrodes  leading  to  a lower 
average  voltage  between  them  and,  therefore,  to  an  efficiency  lower  than 
that  obtained  with  a magnetic  field  of  uniform  amplitude  over  the  entire 
electrode  length  (due  to  added  circulating  currents).  If  this  situation  is 
compensated  for  by  making  the  magnetic  poles  longer  than  the  electrodes,  to 
obtain  a uniform  field  along  the  latter,  then  the  field  intensity  beyond 
the  electrodes  will  not  decay  sharply  enough,  giving  rise  to  currents  which 
increase  the  ponderomotive  force  without  contributing  to  the  electric  power 
drawn  from  the  generator  so  that  the  efficiency  is  not  improved.  The  best 
solution,  according  to  the  results  of  this  study,  is  to  shape  the  magnetic 
field  to  achieve  a field  of  uniform  strength  over  the  entire  length  of  the 
electrodes  and  decaying  (quickly)  beyond  the  electrode  edge  according  to  an 
exponential  curve. 

The  efficiency  was  found  to  be  higher  as  Le  is  increased.  For  a 
channel  without  insulating  vanes,  t.he  increase  in  efficiency  is  not  linearly 
dependent  on  Le.  For  instance  (see  Fig.  V. 29),  for  A = 2,  Le  = 2h  increases 
the  efficiency  by  92%  with  respect  to  the  reference  case  (for  which  le  = 0); 
and  doubling  Le  (to  4h)  brings  about  a further  efficiency  increase  of 
only  2.5%.  The  dependence  shown  in  Fig.  V.29  remains  qualitatively  correct 
for  channels  with  a higher  aspect  ratio;  however,  a relatively  larger 
value  of  Le  would  be  required  to  obtain  the  same  percentages  of  efficiency 
increase. 
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Figure  V.29  The  Influence  of  lf  on  the  maximum  attainable  efficiency 
for  an  "infinitely^long"  insulating  vane,  A * 2 
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It  should  be  noted  that,  If  insulating  vanes  are  present  in  the 
region  beyond  the  electrodes,  the  dependence  of  the  efficiency  on  Le  is 
different.  For  aspect  ratios  of  A * 3 and  5,  for  instance,  an  asymptotic 
value  for  the  maximum  attainable  efficiency  seems  to  be  reached  for  Le  of 
about  3h,  Fig.  V.30. 

The  emphasis  in  the  present  study  was  on  the  influence  of  the  insulating 
vanes  on  the  efficiency.  For  Le  * 2h,  a further  efficiency  increase  of 
10%  can  be  achieved  by  the  insertion  of  a center-line  vane  of  a length  i « 

8h.  Additional  vanes  of  the  same  length  lead  to  a further  increase  In  the 
maximum  attainable  efficiency  of  only  a few  percent,  Fig,  V.28.  This 
fact,  together  with  the  technological  difficulties  implied  by  an  increased 
number  of  vanes,  seem  to  Indicate  that  the  most  efficient  solution  for  the 
end  loss  reduction  would  be  the  insertion  of  a single  insulating  vane 
on  the  center  line.  It  has  to  be  stressed  that  not  only  the  length,  but 
also  the  axial  position  of  the  vane,  is  of  critical  Importance.  With  the 
magnetic  field  "shaped"  to  approximate  an  exponentially-decaying  longitudiual 
distribution,  the  best  efficiency  Is  obtained  when  the  vane  edge  Is  positioned 
at  a distance  of  about  LP/2  (measured  axially)  from  the  electrode  end. 

Since  for  different  load  factors  the  electrostatic  field  created  by  the 
electrodes  has  different  intensities,  the  optimal  position  of  the  vane  Is 
dependent  on  the  load  factor.  However,  the  distance  from  the  electrode  end 
indicated  above  may  be  considered  valid  for  the  range  of  load  factors  which 
correspond  to  maximum  generator  efficiency  (F  s 0.8  ....0.95). 

Obviously,  the  same  increase  in  efficiency  may  be  attained  in  two  or 
more  different  ways.  For  Instance  (see  Fig.  V.32),  for  a channel  with  A a 2 
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and  an  exponential  field  overhang  characterized  by  Le  - h,  the  efficiency 

can  be  increased  by  15%  by  either  a)  a center-line  insulating  vane  of 

length  14h  beginning  at  the  end  of  the  electrodes,  b)  a center-line  vane  of 

length  8h  beginning  at  a better  axial  position,  or  c)  a center-line  vane  of 

length  4h  and  an  increase  in  Le  to  2h. 

The  choice  between  the  last  two  alternatives  would  depend  upon 
technological  considerations.  In  its  last  version,  the  computer  program 
used  in  this  study  can  provide  data  describing  the  dependence  of  the 
maximum  attainable  efficiency  with  respect  to  any  one  of  the  parameters 
characterizing  the  vanes  and  the  magnetic  field.  A few  examples  are 
given  in  Figs.  V.29-V.32.  In  Figs.  V.28  and  V.30,  the  values  of  the  aspect 
ratio  are  comparable  with  the  LMMHD  generator  models  presently  being 
tested  (see  Section  III).  For  the  specific  parameters  chosen  the  maximum 
attainable  efficiencies  are  in  excess  of  80%. 


VI.  LOCAL  MEASUREMENTS 


VI.l  Introduction 


The  two-phase  LMMHD  generator  is  complex  fluid  dynamically  because  it 
employs  two  working  fluids--a  gas  as  the  thermodynamic  fluid  and  a liquid 
metal  to  provide  the  high  electrical  conductivity.  [For  a description  of 
the  cycle,  see  Brunsvold  and  Pierson  (1976),  Pierson  (1978),  or  Pierson  et 
al  (1979).]  The  gas  expands  and  performs  work  on  the  liquid  metal  in  the 
generator,  but  the  very  large  electromagnetic  force  acts  solely  on  the 
liquid.  In  addition,  because  of  the  high  ratio  of  liquid  to  gas  mass  flow 
rates  and  liquid  to  gas  heat  transfer,  the  gas  expansion  occurrs  at  almost 
constant  temperature  (hence  two-phase  LMMHD  cycles  have  high  thermodynamic 
efficiencies,  a significant  attractive  feature)  and  the  enthalpy  change  in 
the  two-phase  mixture  comes  almost  entirely  from  the  liquid.  To  maintain 
high  generator  and  cycle  efficiency,  good  interfacial  heat  and  mechanical 
energy  transfers  are  required,  as  has  been  demonstrated,  see  Fabris  and 
Pierson  (1978)  or  Petrick  et  al  (1978).  These  gas-liquid  interactions  are 
very  dependent  on  the  two-phase  flow  pattern. 

System  studies  have  shown  that  cycle  efficiency  will  improve  greatly 
as  the  average  void  fraction  in  the  generator  is  increased,  because  this 
reduces  the  amount  of  recirculating  liquid  and  associated  loss.  However, 
as  the  void  fraction  is  increased,  two-phase  flow-pattern  transitions  may 
occur,  limiting  generator  (and  cycle)  performance.  Good  generator  perform- 
ance has  been  measured  at  average  void  fractions  up  to  approximately  0.7, 
with  efficiencies  in  excess  of  0.60  and  power  densities  up  to  approximately 
32  MWe/m3  (see  Section  III).  The  use  of  foams  has  been  proposed  by 
Fabris  et  al  (1975)  to  prevent  flow  transitions  and  allow  higher  average 
void  fractions,  but  this  has  yet  to  be  tested  in  a generator. 
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All  reported  experimental  measurements  have  been  terminal  or  average 
parameters.  No  experimental  data  is  available  on  local  gas  and  liquid 
velocities,  bubble  sizes,  local  temperatures,  etc.  Local  measurements  in 
any  two-phase  flow,  with  or  without,  a magnetic  field,  are  difficult,  and  no 
measurements  inside  a LMMHD  generator  with  a magnetic  field  have  been 
reported.  There  is  not  even  published  data  on  hot- film  measurements 
in  pure-liquid  alkali-metal  flows.  Thus,  it  was  necessary  to  develop  the 
experimental  techniques,  and  this  should  be  read  as  an  initial  report,  not 
a presentation  of  data  obtained  by  established  measurement  methods.  The 
experimental  facility  is  described  briefly  in  Section  VI. 2,  and  the  hot- 
film  probe  and  resistivity  probe  results  are  in  Sections  VI. 3 and  VI. 4, 
respectively.  Hot-film  probes  yield  both  velocity  and  void  fraction  data, 
while  resistivity  probes  yield  only  void  fraction  directly  but  are  smaller, 
i.e.,  respond  to  smaller  bubbles,  and  easier  to  use. 


V 1 . 2 Experimental  Facility 

The  ambient-temperature  NaK-nitrogen  LMMHD  facility  used  for  the 
two-phase  LMMHD  measurements  has  been  described  in  Petrick  et  al  (1977). 

The  generator  channel  used  is  387  mm  long,  with  102-mm  electrode  spacing 
and  height  that  varies  linearly,  from  1.41  mm  at  the  inlet  to  1.91  mm 
at  the  outlet,  parallel  to  the  applied  magnetic  field.  The  load  resistance 


Recently  some  LMMHD  resistivity-probe  pipe-flow  data  have  been  published, 
but  they  are  not  in  the  correct  parameter  range.  For  example,  Michiyoshi 
et  al  (1977)  studied  bubbles  in  mercury  at  low  velocities,  less  than  1 m/s, 
where  the  electromagnetic  force  density  is  not  dominant.  Also,  it  is 
important  to  be  aware  that  mercury  behaves  differently  than  NaK  in  two-phase 
flows  because  of  the  substantial  difference  in  their  Weber  numbers  (ratio 
of  inertial  to  surface  tension  forces)  that  govern  the  disintegration  of 
liquid  sheets  [see  experimental  result  of  Dombrowski  and  Fraser  (1954)]. 

The  ratio  of  mercury  to  NaK  Weber  numbers  is  four. 
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is  continuously  variable  and  is  connected  to  the  electrodes  by  means  of 
compensating  bars  located  on  both  sides  of  the  channel.  The  electro- 
magnet has  a 406  mm  diameter  pole  face  and  is  capable  of  flux  densities 
in  excess  of  1.2  T at  the  required  air-gap  length.  Instrumentation  is 
provided  for  measurement  of  terminal  voltages  and  currents,  gas  and  liquid 
flow  rates,  pressures,  temperatures,  and  average  void  fractions. 

For  the  NaK-nitrogen  hot-film  experiments  a TSI  quartz-coated  rugged- 
ized  side-flow  probe  3.18  mm  in  diameter.  Model  1269  W,  was  used.  The 
probe  was  connected  to  a TSI  Series  1050  constant-temperature  anemometry 
system,  and  the  output  was  displayed  on  an  oscilloscope  screen  and  recorded 
on  a Sangamo/Tandberg  Model  TIR  115  four-channel  FM  tape  recorder.  The 
analog  signals  were  subsequently  digitized  at  the  rate  of  40UU  samples  per 
second,  and  the  digital  signals  used  for  some  of  the  figures  in  Section  VI. 3. 
The  minimum  detectable  bubble  size  is  limited  to  approximately  10  mm  by  the 
probe  size,  sample  rate,  and  two-phase  velocity. 

Appreciable  dc  drift  was  observed  in  the  output  signal  from  the 
hot-film  probe,  but  is  not  considered  important  here  because  only  qualita- 
tive behavior  and  gas- liquid  changes  are  of  interest,  rather  than  absolute 
velocities.  (After  the  experiments  reported  herein,  the  probe  was  gold- 
plated  by  vapor  deposition,  resulting  in  reduced  dc  drift  because  of 
improved  NaK  wetting  of  the  probe.)  A sample  calibration  curve,  taken  from 
the  digitized  data,  is  shown  in  Fig.  VI. 1 only  to  give  an  indication  of  the 
voltage-velocity  magnitudes  for  the  data  in  Section  VI. 3.  (The  voltage  is 
the  anemometery  system  output  with  the  probe  at  the  center  of  the  channel 
for  B * 0 and  a pure-liquid  flow;  the  velocity  is  calculated  as  the  average 
velocity  from  the  measured  volume  flow  rate  multiplied  by  1.225,  the  1.225 
being  the  correction  factor  based  in  the  velocity  profile  shape  in  a 
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curve  of  the  hot-film  probe  and  anomometry  system 


rectangular  channel  at  high  Reynolds  numbers.)  The  voltages  shown  in 
Fig.  VI. 1 must  be  multiplied  by  four  for  use  with  the  oscilloscope  signals. 

The  resistivity  probes,  basically  thin  wires  insulated  except  at  the 
very  tip,  were  made  by  one  of  the  authors.  To  improve  the  void  fraction 
and  bubble  frequency  measurements,  a discriminator  circuit  was  developed  by 
Dunn  (1978)  that  produces  a square-wave  signal  with  the  value  1 V when  the 
probe  is  in  gas  and  0 V when  the  probe  is  in  liquid;  thus,  a dc  voltmeter 
connected  to  this  signal  reads  local  void  fraction  directly.  Note  that  the 
time  response  of  the  probe  is  typically  one  order  of  magnitude  faster 
with  NaK-nitrogen  than  with  air-water. 

A resistivity  probe  with  the  discriminator  circuit  was  calibrated 
against  the  more  established  y-ray  attenuation  measurements  by  Gawor 
(1978).  Agreement  to  within  three  percentage  points  was  obtained  where 
bubble  diameters  were  greater  than  approximately  0.5  mm. 


VI. 3 Hot-Film  Probe  Results 

Oscilloscope  traces  of  the  hot-film  probe- anemometry  system  output 
signals  are  shown  in  Fig.  VI. 2 for  pure-NaK  non-MHD  flows  (B  = 0)  with  the 
probe  at  the  center  of  the  channel  cross  section,  i.e.,  halfway  between 
the  electrodes  and  centered  between  the  insulating  walls,  and  250  mm  from 
the  beginning  of  the  electrodes.  As  expected,  a substantial  increase  in 
the  frequency  and  the  intensity  of  the  fluctuations  is  apparent  as  the  mass 
flow  rate,  rtit,  is  increased.  The  influence  of  B is  shown  in  Figs.  VI. 3 
and  VI. 4,  as  measured  25  and  38  mm,  respectively,  from  the  center  of  the 
channel,  i.e.,  toward  one  electrode.  (The  distance  from  the  center  to 
either  electrode  is  51  mm.)  The  magnetic  field  clearly  damps  the  fluctua- 
tions, although  some  weak  large-scale  "eddies"  remain.  Comparison  of  the 


-129- 


Figure  VI. 3 Hot-film  probe  signals  in  pure  NaK  flow:  measured  38  mm  from  the  channel  center- 
line,  5 ms/DIV  on  abscisa,  10  mV/DIV  on  ordinate 


curves  is  limited  because  of  signal  drift,  but,  for  the  same  conditions, 
there  is  more  fluctuation  closer  to  the  electrodes.  Similar  measurements 
in  mercury  have  been  reported  previously  [see,  for  example,  Gardner  and 
Lykoudis  (1971)]. 

Plots  of  the  digitized  probe  signals  for  two-phase  NaK-nitrogen  flows 
are  shown  in  Figs.  VI. 5 to  VI. 10,  all  taken  at  the  channel  center.  The 
square  of  the  bridge  output  voltage  from  the  anemometer  is  proportional  to 
the  heat  transfer  rate;  thus,  the  high  voltages  occur  when  the  probe  is  in 
contact  with  NaK  and  the  low  voltages  (except,  possibly,  for  a liquid  film) 
when  the  probe  is  in  contact  with  nitrogen.  The  rate  of  voltage  increase 
is  greater  than  the  rate  of  voltage  decrease,  because,  when  the  probe 
starts  to  encounter  liquid,  the  increase  is  governed  by  the  time  constant 
of  the  probe  and  the  electronics,  whereas,  when  the  probe  starts  to  en- 
counter gas,  the  probe  remains  wetted  by  a liquid  film  (at  least  for  a 
short  time).  The  frequency  response  of  the  recording-digitizing  process 
was  not  determined,  hence,  the  plots  should  be  regarded  only  as  qualitative 
data. 

To  interpret  the  hot-film  sigr.als,  it  is  helpful  to  refer  to  previous- 
ly obtained  results  with  a conical  hot-film  orobe  0.5  mm  in  diameter  in  an 
air-water  flow  with  and  without  the  addition  of  small  amounts  of  liquid 
soap.  Without  the  soap,  the  bubbles  were  approximately  3-4  mm  in  diameter 
and  were  clearly  visible  in  the  output  signal,  whereas  with  the  soap  the 
bubbles  were  much  smaller  in  diameter  and  very  few  bubbles  were  indicated 
by  the  signal. 

Figures  VI. 5,  VI. 6,  and  VI. 7 were  all  obtained  for  B = 0.  The  first 
two  are  for  the  same  flow  conditions,  but  Fig.  VI. 6 shows  the  signal  at  a 
time  when  larger  numbers  of  smaller  bubbles  were  observed.  In  Fig.  VI. 7, 
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Figure  VI. 9 Hot-film  probe  signals  in  NaK-N?  flow:  measured  at  the  channel 
centerline,  B = 0.6  T,  m_  = 0.065  kg/s,  m = 5.0  kg/s 


the  gas  flow  rate  is  four  times  as  high  and  corresponds  to  an  average  void 
fraction  of  approximately  0.65.  In  all  cases,  bubbles  (voids  may  be  a 
more-appropriate  term,  because  they  are  not  necessarily  spheres)  appear  to 
differ  in  size  and  distribution  with  time.  Especially  for  the  last  figure, 
there  may  be  many  pieces  of  liquid  (drops,  threads,  sheets,  etc.)  inter- 
mingled with  the  voids;  this  is  characteristic  of  churn-turbulent  flow. 

Figures  VI. 8,  VI. 9,  and  VI. 10  were  measured  for  two-phase  MHD  flows 
with  a load  resistance  of  0.1  mn,  a relatively-small  resistance  value  that 
results  in  heavy  loading  of  the  channel,  and  magnetic  flux  densities  of  B = 
0.3,  0.6,  and  1.2  T,  in  that  order.  The  mass  flow  rates  of  Fig.  VI. 8 are 
the  same  as  those  of  Figs.  VI. 5 and  VI. 6,  whereas  Figs.  VI. 9 and  VI. 10  are 
close  to  the  mass  flow  rates  of  Fig.  VI. 7.  For  all  cases  with  an  applied 
magnetic  field,  especially  at  the  higher  flux  densities,  there  is  an 
apparent  decrease  in  the  number  of  detectable  bubbles.  Also,  the  fluctua- 
tions in  the  liquid  flow  appear  to  increase  in  magnitude  as  B increases, 

and  the  higher-frequency  components  appear  to  be  damped.  The  cause  of  the 

apparent  decrease  in  the  number  of  bubbles  is  unknown  at  present.  The 

probe  was  moved  to  different  locations  (but  still  halfway  between  the 
insulating  walls)  with  similar  results.  Note  that,  for  two-phase  flows, 
the  fluctuations  appear  to  increase  in  magnitude  with  increasing  B,  whereas 
for  pure-liquid  flows,  the  fluctuations  are  damped. 

VI .4  Resistivity  Probe  Results 

The  resistivity  probe  used  to  obtain  the  data  shown  in  Figs.  VI. 11  to 
VI. 13  was  inserted  through  a pressure  tap  located  250  mm  downstream  of  the 
beginning  of  the  electrodes.  The  probe  used  a stainless  steel  wire  with  an 
uninsulated  tip  approximately  0.1  mm  in  diameter  and  0.25  mm  in  length,  on 
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a stem  approximately  5 mm  long,  which  was  pointed  into  the  flow.  (Note 
that  the  length  of  the  uninsulated  tip  is  also  the  dimension  of  the  small- 
est detectable  bubble  or  void  parallel  to  the  flow.)  It  is  convenient  to 
use  bubble  frequency  (bubbles/second)  to  interpret  the  resistivity  orobe 
data,  because  the  discriminator  circuit  converts  the  signal  to  square  waves 
that  are  easy  to  count. 

Bubble  frequencies  for  m4  = 5 kg/s  and  tog  = 0.05  kg/s  are  shown  in 
Fig.  VI.  11  as  functions  of  B and  the  distance  of  the  probe  from  the  center 
of  the  channel  along  the  line  to  one  electrode  (total  distance  of  51  mm). 
For  B = 0,  the  bubble  frequencies  show  a triple  maximum,  similar  to  the 
void-fraction  profiles  obtained  with  a resistivity  probe  by  Fabris  et  al 
(1977b)  for  air-water  flow  with  an  identical  mixer  geometry.  As  B is 
increased,  the  bubble  frequency  drops  sharply  to  very  low  values.  Note 
that  the  decrease  in  bubble  frequency  with  increasing  B is  less  closer  to 
the  electrodes;  this  might  possibly  indicate  bubble  migration  toward  the 
electrodes,  as  would  be  expected  from  the  end-region  currents  that  produce 
M-shaped  MHD  velocity  profiles.  (Recent  measurements  for  a similar  LMMHD 
generator  confirm  the  existence  of  the  pressure  distribution  in  the  end 
regions  caused  by  these  currents,  see  Section  III. 4.)  The  addition  of  one 
or  more  insulating  vanes  in  the  end  regions  to  supress  end  currents  would 
also  suppress  this  bubble  migration;  thus  migration  is  not  a matter  of 
concern  in  regard  to  practical  LMMHD  generators. 

The  most-significant  feature  shown  in  Fig.  VI.  11  is  the  low  bubble 
count  at  all  locations  for  high  B;  this  confirms  the  hot-film  data.  Bubble 
migration  toward  the  electrodes  cannot  be  the  sole  cause,  because  so  many 
bubbles  there  would  break  the  current  paths.  The  most-likely  explanation 
is  that  the  large  electromagnetic  force  density  and  the  resulting  high 
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pressure  gradient  "shear"  the  gas  bubbles  or  voids  into  smaller  pieces, 
which  may  be  too  small  to  be  detected  clearly  by  either  probe.  (Such  small 
voids  could  cause  the  high-frequency  fluctuation  shown  in  Figs.  VI. 8 
through  VI. 10  for  the  hot-film  probe.)  Evidence  of  void  break-up  is 
discussed  below.  Another  factor  could  be  the  influence  of  the  electro- 
magnetic force  density  on  the  void  shape  and  motion  around  the  probe. 

Also,  the  voids  are  smaller  because  of  the  higher  static  pressure  at  the 
probe  location  with  B t 0;  however,  this  is  not  enough  of  a factor  in 
itself. 

Additional  experiments  were  performed  to  clarify  some  of  the  effects. 
At  low  flow  rates  or  velocities,  the  electromagnetic  force  density  is 
reduced;  thus,  changes  with  B becomes  more  gradual  and  easier  to  detect. 

In  addition,  voids  are  easier  to  detect  at  lower  velocities.  Data  are 
shown  in  Fig.  VI. 12  for  flow  rates  one-fifth  that  of  Fig.  VI. 11  for  the 
same  quality  (same  void  fraction  if  pressure  and  slip  ratio  are  the  same). 
At  B = 0,  the  bubble  frequency  is  very  low,  indicating  a slug  flow  pattern 
with  large  voids.  As  B is  increased,  the  bubble  frequency  increases  by 
more  than  one  order  of  magnitude  before  decreasing.  This  indicates  that 
the  large  voids  are  broken  up  into  smaller  voids  and  many  are  not  detected 
at  large  B values,  as  discussed  above.  Significantly,  more  bubbles  are 
observed  close  to  the  electrode,  as  occurred  with  the  higher  flow  rate. 

Data  are  shown  in  Fig.  VI. 13  for  a high  flow  rate  but  with  the  gene- 
rator open-circuited,  so  that  the  electromagnetic  force  density  is  greatly 
reduced.  Again,  the  bubble  frequency  increases  as  B is  increased;  this 
time  without  a decrease  at  larger  B values. 
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VII.  CONCLUSIONS 

The  new  generator  channel,  LT-4,  has  been  tested  at  high  velocities 
and  high  pressures  (large  magnetic  interactions).  Power  densities  of  up  to 
32  MWe/m^  were  obtained  in  this  small,  20-kWe  generator,  equal  to  or 
better  than  those  anticipated  for  LMMHD  generators  in  practical  power 
systems.  Efficiencies  higher  than  0.6  at  average  void  fractions  well  above 
0.5  were  measured,  which  is  very  good  performance  for  a small,  wall -dominated 
generator.  The  rapid  decrease  in  slip  ratio  with  increasing  flow  rate, 
demonstrated  previously  at  low  power  densities  [Petrick  et  al  (1978)],  was 
confirmed  for  high  power  densities.  In  addition,  more-detailed  pressure 
and  voltage  data  were  obtained.  These  postive  results  meet  the  goals  of 
the  revised  facility  and  channel  LT-4,  and  justify  the  existence  of  this 
experimental  program. 

The  use  of  foams  in  the  two-phase  LMMHD  generators  allows  the  generator 
to  operate  at  higher  average  void  fractions,  and  leads  to  higher  "equivalent 
turbine"  efficiencies,  reduced  losses,  because  there  is  less  liquid  to 
recirculate,  and  higher  system  performance.  The  foam  experiments  have 
concentrated  on  bubble  tests  and  surface  tension  measurements,  and  much 
basic  data  has  been  obtained.  Barium  has  been  identified  as  an  attractive 
additive  for  NaK  because  it  1)  markedly  stabilizes  bubble  lifetime,  2)  pro- 
motes foaming,  and  3)  is  a wetting  agent.  Although  barium  is  known  to 
react  with  nitrogen  (but  not  helium),  no  degradation  of  bubble  performance 
was  observed  at  temperatures  close  to  ambient;  some  slow  degradation  was 
observed  at  temperatures  above  approximately  400  K,  with  the  rate  of 
degradation  increasing  with  temperature.  Bismuth  was  identified  as  a 
back-up  for  barium.  A guiding  hypothesis  for  understanding  the  phenomena 
contributing  to  foamability  in  liquid  metal  systems  has  been  developed. 


End  current  losses  in  LMMHD  generators  can  be  controlled  by  the 
proper  choice  of  1)  channel  aspect  ratio,  2)  magnetic  flux  density  distribu- 
tion outside  of  the  electrodes,  and  3)  number  and  dimensions  of  insulating 
vanes.  A computer  code  has  been  developd  for  the  calculation  of  the  end 
current  loss  as  the  above  parameters  are  varied.  Initial  studies  with  this 
code  have  shown  that  generator  efficiencies  in  excess  of  0.8  are  attainable. 

Initial  measurements  of  local  flow  parameters  in  two-phase  LMMHD  flows 
have  been  presented.  Both  hot-film  and  resistivity  probes  were  used,  and 
similar  results  were  obtained.  Clearly,  the  behavior  of  the  voids  is 
influenced  when  the  electromagnetic  interaction  parameter  is  high,  although 
definitive  statements  are  limited  at  present  by  the  inability  of  the  probes 
to  detect  the  smaller  voids  that  result  for  higher  electromagnetic  force 
densities.  There  is,  however,  definite  evidence  of  void  break-up  and 
migration. 

Significant  advances  have  been  made  in  understanding  the  operation  of 
two-phase  LMMHD  generators.  All  results  are  very  encouraging,  and  indicate 
that  large  high-efficiency  two-phase  LMMHD  generators  can  be  developed. 
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R assumed  1 for  these  runs,  not  measured 


9.038  .254  1.672  .541  1.338  .955  11.015  .537  .306 


APPENDIX  B 


GENERATOR  POWER  EQUATIONS 

The  electrical  power  output  is  determined  from  the  measured  generator 
terminal  voltage  V|_  and  the  load  resistance,  Rl,  i.e., 

2 

Electrical  power  output  = V[_  /Rj_ . (B.l) 

The  ideal  liquid  mechanical  pump  work,  the  power  input  from  the  liquid  to 
the  generator,  is  the  pump  work  plus  the  kinetic  energy  change  of  the 
liquid  (assuming  there  are  no  losses  and  that  the  pumping  process  is 
isentropic), 

MP  + rti .(u?  ,„-u?  t) 

-^,(8.2) 


Ideal  liquid  mechanical  pump  work 


T 


where  mt  is  the  liquid  mass  flow  rate,  pt  the  liquid  density,  vp  = pin 
“ Pout  the  pressure  difference  across  the  generator  for  this  idealized 
case,  and  ut  the  liquid  velocity.  The  total  power  input  is  the  enthalpy 
and  kinetic  energy  changes  across  the  generator  for  both  the  liquid  and  gas 
components, 

2 2 
in"Ut  out ^ 

Total  power  input  = rfit  [ah^  + ! ] + 


2 2 

n.g  uhg  + lu3»1f  g-."*’  ] 


(B.3) 


where  vh  = hjn  - hout  denotes  the  enthalpy  change.  The  enthalpy  change 
for  the  liquid  is 


lh,  ' S.  (Tin  ' T«,t>  * <’  - *T  T>f  - 


(B.4) 
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where  T is  the  absolute  temperature,  cpt  the  specific  heat,  the 
superscript  denotes  an  average,  and 


* pa  C 


For  an  ideal  gas. 


4h9  * Cpg  (T-n  - Tout>  ' 


(B.5) 


(B.6) 


The  isentropic  power  input  is  found  by  assuming  that  the  expansion  of 
the  two-phase  mixture  through  the  generator  is  isentropic, 

as  = o = a [C  in Jai <■ [In j + 

•*  Kn  in  y 


iht[C_  In  ^ - !l  (p  p.  n 
1 pt  Tin  pt  out 


(B.7) 


This  leads  to 


T f « T-  exp  4 aT<P°ut  'Pi"’  + *9R9MpcuAn) 
out  ‘in 


1Vp9  ♦ *.  v 


(B.8) 


for  the  exit  temperature  of  the  two-phase  mixture.  The  isentropic  power 
input  is  calculated  from  Eq.  B.3  using  Tout  from  Eq.  B.8. 
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THE  USE  OF  ELLPACK77  FOR  SOLVING  THE  LAPLACE  EQUATION  ON  A REGION  WITH 


INTERIOR  SLITS,  APPLICATION  TO  A PROBLEM  IN  MAGNETOHYDRODYNAMICS 

2 

j1  Gherson  , an 
July  18,  1978 


12  2 
Robert  E.  Lynch  , Paul  Gherson  , and  Paul  S.  Lykoudis 


1.  Introduction.  We  describe  modifications  of  the  Purdue  University 
version  of  ELLPACK77  (see  Rice  [1977])  to  solve  the  Laplace  equation 
subject  to  mixed  boundary  conditions  on  the  boundary  of  a rectangule 
and  Neumann  conditions  on  straight  line  segments  in  the  interior  of 
the  rectangle. 

In  Section  2 a brief  description  is  given  of  the  physical  origin 
of  the  problem,  namely  the  analysis  of  the  performance  of  a 
magnetohydrodynamic  electric  generator.  Section  3 contains  a 
mathematical  model  of  the  generator.  Section  4 contains  a new  result: 
we  show  that  one  can  determine  the  efficiency  of  the  model  for  alj,  values 
of  one  of  its  parameters  by  solving  a single  pair  of  boundary  value 
problems.  Section  5 very  briefly  describes  how  the  efficiency  of 
the  generator  can  be  increased  by  inserting  insulating  vanes  in 
the  fluid;  in  the  mathematical  model,  these  are  modeled  by  slits 
in  the  interior  of  the  rectangular  domain  on  which  Neumann  boundary 
conditions  are  specified.  Section  6 contains  a brief  outline  of  the 
ELLPACK  system  and  its  application  to  this  problem.  Section  7 gives 
the  complete  system  of  finite  difference  equations  whose  solution  is 
taken  as  an  approximation  to  the  electric  potential  in  the  generator. 

The  ELLPACK  module  5-POINT  STAR,  which  was  written  by  Ronald 
F.  Boisvert  of  Purdue  University,  was  modified  to  solve  this  problem. 

We  thank  him  for  his  help  during  our  modification  of  his  routines. 
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Departments  of  Computer  Science  and  Mathematics,  Purdue  University. 
Department  of  Nuclear  Engineering,  Purdue  University. 
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2.  Physical  phenomenon ♦ In  this  section  we  describe  the  idealized 
physical  situation  which  we  consider. 

An  inviscid  incompressible  fluid  flows  with  constant  speed  v in 
the  x-direction  in  a channel  between  a pair  of  planes  at  y * -h 
and  at.  y = h.  All  variation  in  the  z-direction  is  neglected.  The 
fluid  is  electrically  conductive  with  uniform  conductivity  o. 

Except  between  x = -xg  and  x = xg,  the  channel  walls  are 
electric  insulators.  Between  -xg  and  xg  they  are  electrodes  with 
infinite  electric  conductivity. 

A magnetic  field  of  intensity  B is  imposed  in  the  z-direction. 

It  is  uniform  across  the  channel  and  its  magnitude  depends  only  on 
x.  Its  magnitude  is  symmetric  with  respect  to  x,  it  is  constant 
for  x between  -xa  and  x . and  it  tends  to  zero  as  x tends 
to  infinity;  a special  case  is  B(x)  « 0 for  jx|  > xg. 

The  Lorentz  force  which  acts  on  the  electric  charges  in  the 
fluid  causes  charge  separation  and  charges  of  different  sign 
collect  near  the  walls  of  the  channel.  This  creates  an  electric 
field  in  the  fluid  with  potential  0.  We  use  20q  to  denote 
the  resulting  potential  difference  of  the  electrodes  and  take  0Q 
to  be  the  potential  of  the  electrode  on  the  wall  at  y = h.  The 
channel  height  is  2h  and  it  follows  from  Maxwell's  equations  that 
20q  is  given  by 

20o  « 2hvB0 

The  device  becomes  an  electric  generator  when  wires  are  attached 
to  the  electrodes  and  ends  of  an  electric  resistance;  this  produces  an 
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electric  current  in  the  wires  and  resistance.  The  potential  difference 
between  the  electrodes  is  reduced  to  a fraction  n of  its  open 
curcuit  value  20Q;  the  fraction  depends  on  the  load  and  n is 
called  the  load  factor. 

The  current  density  J in  tl. . fluid  is 

J « o(F  + vxE)  = a(-grad  t>  + vxtf) . 

We  write  the  magnitude  of  the  magnetic  field  as  BQb(x)  with  b(x) 
equal  to  unity  for  x between  -xg  and  xg.  The  z-component 
jz  of  the  current  density  is  zero  and  its  other  two  components  are 

Jx(x,y)  * cj(-3(Kx .y)/3x) , 

jy(x,y)  = o[-30(x,y)/3y  + vBQb(x)]. 

We  call  the  term  in  j which  involves  (5  the  electrostatic  contribution 

•J 

and  the  term  which  involves  the  magnetic  field  the  magnetic  contribution. 

When  thei*a  is  an  external  current,  then  j is  positive  in 

the  vicinity  of  the  electodes  because  there  the  magnetic  contribution 

dominates  the  electrostatic  one.  The  magnetic  field  decreases  rapidly 

as  x increases  from  x and  for  large  x,  j is  negative  because 

“ y 

there  the  electrostatic  contribution  dominates  the  magnetic  one  and 
this  give  rise  to  a counter-current.  The  net  amount  of  current  which 
can  be  drawn  from  the  system  is  diminished  by  this  counter-current. 

It  is  part  of  cur  research  to  determine  the  effect  on  the  performance 
predicted  by  the  mathematical  model  of  thin  insulators  (vanes)  placed 
in  the  fluid;  as  expected,  these  vanes  decrease  the  counter-current. 
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To  determine  the  net  current  which  can  be  drawn  from  the  system 
per  unit  channel  width  (in  the  z-direction) , one  integrates  j 
along  a line  parallel  to  the  x-axis  and  by  symmetry,  one  can  integrate 
2jy  from  x = 0 to  x = The  value  obtained  is  independent  of 
y for  y between  y = -h  and  y = h.  The  total  electric  power  P 
per  unit  channel  width  is  the  product  of  the  current  and  2o0q  = 2nvBgh. 
One  obtains 

P = 4nvBQho  [-3*Kx,y)/3y  + vBQb(x)]  dx,  P independent  of  y. 

The  total  power  W per  unit  channel  width  needed  to  maintain  the 
flow  is  the  integral  of  vj^BQb  over  the  whole  channel.  This  is 
given  by 

W = 4 vo  /p  dy  /“  dx  [-3tf(x,y)/3y  + vBQb(x)]  BQb(x). 

Considering  only  these  two  powers,  the  efficiency  E of  the  device 
is  taken  to  be 


E = P/W. 

[We  remark  that  one  of  the  things  which  is  neglected  in 

this  model  is  the  effect  of  viscosity  which  gives  rise  to  friction 

effects  and  makes  the  fluid  speed  nonuniform  across  the  channel. 

A second  thing  which  is  neglected  is  the  variation  of  the  magnetic 

field  across  the  channel  and  between  -x„  and  x^.  A third 

e e 

thing  which  is  neglected  is  the  finite  width  of  the  channel.] 
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3.  Mathematical  model.  We  approximate  the  infinite  channel  with  a 

finite  one  with  end  points  ±xm,„  where  xmiU  is  much  larger  than 

r max  max  3 

either  x or  h.  The  mathematical  model  of  the  Idealized  situation 
e 

described  in  Section  2 is  taken  to  be 


(3-la) 
(3-lb) 
(3-lc) 
(3-ld) 
(3-le) 
(3-1 f) 


d2t{x>y)/dxZ  + 92<>(x,y)/3y2  = 0, 


0(x,h)  = nhvBQ, 
30(x,h)/3y  = vBQb(x), 


*(xmax’y)  = °’ 


0(x,O)  = 0, 
30(O,y)/3x  = 0, 


0<x<xmax‘  0<V<h» 


0 < x < xe, 


xe  < x < xmax* 


0 < y < h, 

R3  3T1  * 

0 < x < x , 

”■  “ max  * 

0 < y < h. 


The  function  b which  specifies  the  magnetic  field  is  equal  to 
unity  for-  x between  0 and  x and  in  some  of  our  experiments  it 

C 

was  equal  to  zero  for  x larger  than  xg;  in  others  it  was 


(3-2) 


b(x) 


0 < x < x. 


(exp(-[x-xe]/efold) , *e<x 


where  efold  is  a constant  which  specifies  the  rate  of  decay  of 
the  magnetic  field  beyond  the  end  of  the  electrode. 

The  channel  half-width  h is  taken  as  the  unit  of  length  and 
we  set 

x = hX,  y = hY,  *(x,y)  = rf(hX.hY)  = ♦'(X.Y), 


x e " ^ e’  xmax  ~ ^max  ’ 


b(x)  = b(hX)  = B(X),  efold  - hEFOLD 
30(x,y)/3x  = (l/h)D4>'  (X,Y)/3X , i»i(x,y)/3y  * (l/n)3$' (X,Y)/3Y. 
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With  the  exception  of  (3-lc),  the  equations  (3-1)  and  (3-2)  are 
transformed  by  making  the  changes: 


x -*•  X,  y -»•  Y,  b-^B,  h -♦■  1,  efold  -*■  EFOLD. 


Equation  (3-lc)  becomes 


3*'(X,1)/3Y  = hvBQB(X) , 
Now  set  <J> ' = hvBQ«l»  to  get  the  system 
(3-3a) 


32<f(X,Y)/3X2  4-  32*>(X,Y)/3Y2  = 0, 


(3-3b) 

(3-3c) 

(3-3d) 

(3-3e) 

(3-3f ) 
(3-4) 


4>(X,1)  = n. 
34>{X,1)/3Y  = B(X), 

*<Xmax‘Y)  * °- 
*(X.O)  « 0, 


34>  ( 0 , Y )/3X 
B(X) 


= 0, 


= < 


1. 


X C X c Xmav. 
e max 


0 < X 4 W 0 < V < ’• 


0 < X < X , 

»r-.  ii-  O 


Xe<X<Xmax’ 


0 < Y < 1, 


0 < X < Xmav, 
*“  — max 


0 < Y < 1, 


0 < X < X , 
«*  *•-  0 


(exp(-[X-X  1/EF0LD) , X <X. 

6 C 


The  power  P and  the  rate  of  work  W then  become 


(3-5a)  P = d4'V2BqO  fQmX  [-W(X,Y)/3Y  + B(X)]  dX 


.2n2 


(3-5b)  W = 4vtBjo  /0  dY  /0max  dX  [-34>(X,Y)/3Y  + B(X)]  B(X) 

Thus,  in  this  model,  the  efficiency  E = P/W  does  not  depend 
2 2 

on  the  value  of  v BQo. 

An  analysis  of  this  mathematical  model  by  conformal  mapping  is 
given  by  Sutton,  Hurwitz,  and  Poritsky  [1961]. 
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4.  Efficiency  as  a function  of  load  factor.  In  this  section 
we  present  a result  which  does  not  seem  to  be  in  the  literature  which 
significantly  reduces  the  amount  of  computation  to  determine  the 
efficiency  E - P/W  as  a function  of  load  factor  n.  Namely,  for 
a fixed  geometry  and  magnetic  field,  the  efficiency  can  be  determined 
by  solving  a single  pair  of  boundary  value  problems. 

Suppose  that  the  geometry  and  the  function  B are  fixed.  Then 

the  solution  of  (3-3)  depends  only  on  the  value  of  n.  Let  $1  denote 

the  solution  when  n is  equal  to  some  given  value  n-j . Consider 

the  solution  for  an  arbitrary  value  of  n.  We  can  write 


* = *1  + (n  - V 


where  T 

satisfies 

(4-la) 

02f(X,Y)/3X2  + 924'(X,Y)/3Y2  = 0, 

0<X<W  0 < Y < 1. 

(4-lb) 

Y(X,1)  = 1, 

°<  x < xe. 

(4-lc) 

Of (X ,1 )/0Y  = 0 

Xe<  X<  W 

( 4- Id) 

*<WY>  ■ °- 

0 < Y < 1 , 

(4-1  e) 

* 

O 

II 

O 

X 

0 < X < X , 

— - max 

(4-1 f) 

of(n,Y)/ox  = 0, 

0 < Y < 1, 

The 

efficiency  E(n)  for  a given  value 

of  load  factor  n is  then 

given  by 

(4- 2a) 

E(n)  - n[a  + (n  - n-j )6]/Cy 

+ (n  - 1^)6] 

where  the  constants  a,8,y,6  are  given  by 
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(4- 2b) 
(4- 2c) 
(4-2d) 
(4-2e) 


a = /0max  C-34»1  (X ,Y)/DY  + B(X)]  dX 
8 = /0maX  *(X,Y)  dX 

Y = /0  dY  C-3*,  (X,Y)/aY  + B(X)]  B(X)  dX 
6 * f0  dY  /omaX  C-3Y(X,Y)/3Y]  B(X)  dX 
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5.  Reduction  of  the  counter-current.  The  efficiency  of  the  device 
is  limited  by  the  amount  of  counter-current  mentioned  in  Section  2. 

This  counter-current  can  be  reduced  by  inserting  insulating  vanes 
parallel  to  the  X-axis  in  the  fluid.  These  vanes  are  symmetric 
with  respect  to  X = 0 and  to  Y = 0 and  a pair  might  be  on  the  X-axis. 

We  denote  the  left  (L)  and  right  (R)  endpoints  of  these  vanes  in 
the  first  quadrant  by 

(X{L,k),y{k))  and  {X(R*k\y^),  k = 1,2,... 

The  mathematical  models  (3-3)  and  (4-1)  are  augmented  by  the  boundary 
conditions 


(5-la) 


3«f(X,Y(k))/3Y  = B(X) , X(L,k)  < X < X(R,k\  k = 1,2,... 


(5-lb) 


3H'(X,Y(k))/3Y  = 0, 


X^L,k)  < X < X(R’k>,  k = 1,2,... 


If  one  of  these  vanes,  say  for  k - 1,  is  on  the  X-axis,  then 
(3-3e)  and  (4-le)  are  replaced  with 


(5-lc) 


<i>(X,0)  = H'(X ,0)  =0,  0 < X < x^L,k),  X^R,k)  < X < X, 


“ max 


together  with  (5-la)  and  (5-lb)  with  k = 1. 

The  thickness  of  the  vanes  is  neglected  in  the  mathematical  model 

and  thus  the  rectangle  [0,Xmax]x[0, 1]  contains  slits.  Although 
the  Y-derivative  of  the  potential  is  continuous  across  these 
slits  (because  of  (5-1  ),  the  potential  itself  is  discontinuous 
at  the  slits. 
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tool  for  testing  the  performance  of  the  various  components  which  make 
up  a program  designed  to  obtain  approximations  to  solutions  of  elliptic 
partial  differential  equations.  ELLPACK  is  made  up  a a number  of 
modules , for  example,  the  DISCRETIZATION  module  constructs  a set  of 
linear  algebraic  equations  whose  solution  gives  an  approximation  to 
the  solution  of  the  user  specified  partial  differential  equation. 

Each  module  has  several  different  versions,  for  example  there  are 
several  different  discretizations  one  can  use  which  vary  from  the 
standard  divided  central  difference  approximation  to  the  most  recently 
developed  methods  of  approximating  elliptic  partial  differential  equations. 
There  have  been  numerous  contributors  to  the  components  which  make 
up  ELLPACK  and  these  include  people  at  several  different  universities. 

The  current  version  of  ELLPACK,  called  ELLPACK77,  treats  problems 
with  the  domain  of  the  partial  differential  equation  the  interior 
of  a rectangle  or  a cube.  A new  version  of  ELLPACK  is  being  prepared 
which  will  treat  general  domains. 

An  ELLPACK  progrm  consists  of  (a)  a series  of  statements  written 
in  the  ELLPACK  user  oriented  language  grouped  to  form  ELLPACK  segments 
and  (b)  a set  of  user  supplied  Fortran  FUNCTION  and  SUBROUTINE  subprograms. 

The  ELLPACK77  segments  are: 

(i)  EQUATION,  which  specifies  the  elliptic  partial  differential 
equation  to  be  solved; 

(ii)  BOUNDARY,  which  specifies  the  rectangular  domain  and  the 
conditions  on  the  solution  of  the  partial  differential 
equation  at  the  boundaries. 
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(iii)  GRID,  which  specifies  the  number  of  vertical  and  horizontal 
mesh  lines,  a rectangular  mesh  is  generated; 

(iv)  OPTIONS,  which  specifies,  for  example,  the  amount  of 

ELLPACK  generated  information  about  the  execution  of  the 
program; 

(v)  DISCRETIZATION,  which  specifies  which  of  the  discretization 
method  is  to  be  used; 

(vi)  INDEX,  which  specifies  the  indexing  of  the  algebraic 
equations  and  the  unknowns; 

(vii)  SOLUTION,  which  specifies  which  of  several  linear 
algebraic  equation  solvers  is  to  be  used; 

(viii)  OUTPUT,  which  specifies  the  ELLPACK  generated  output, 
for  example:  print  a table  of  the  va^es  of  the 
approximation  at  mesh  points,  construct  a contour 
plot  of  the  approximation; 

(ix)  SEQUENCE,  which  specifies  the  order  of  the  exectuion  of 
the  ELLPACK  segments  and  the  number  of  times  the 
sequence  is  to  be  executed; 

(x)  FORTRAN,  which  specified  that  the  lines  of  the  program 
which  follow  are  user  supplied  Fortran  FUNCTION  and 
SUBROUTINE  subprograms,  these  are  for  the  boundary 
conditions  and  other  needed  routines  as  well  as 
a SUBROUTINE  called  TEST  which  is  described  below. 

The  current  version,  ELLPACK77,  of  ELLPACK  can  be  used  to  generate  an 
approximation  U to  the  potential  or  f,  the  solutions  of  (3-1)  and 
(4-1).  Paul  Gherson  of  Purdue  used  ELLPACK77  to  find  estimates  of  <t> 
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for  various  values  of  the  electrode  endpoint  Xg,  the  endpoint  of  the 
rectangle  Xmax,  different  magnetic  fields  B,  different  load  factors 
and  for  various  conditions  on  Y = 0,  namely  the  boundary  condition 
(3-3e)  and  the  condition  (5-1)  for  k = 1 and  Y^  = 0.  His 
results  are  given  in  the  report  by  Gherson  and  Lykoudis  [1978]. 

One  of  the  discretization  modules  5-POINT  STAR  which  was 
written  by  Ronald  F.  Boisvert  was  modified  by  Robert  E.  Lynch  to 
handle  the  case  of  insulator  vanes  in  the  interior  of  the  rectangular 
domain.  In  addition  to  modifying  the  routines  ased  by  5-POINT  STAR, 
Gherson’ s BOUNDARY  segment  was  changed  slightly  and  an  additional 
DISCRETIZATION  module  was  specified  in  order  to  increase  the  sizes  of 
certain  of  the  ELLPACK  generated  arrays  so  that  the  modifed  5-POINT 
STAR  would  operator  properly.  Also,  Gherson ’s  FORTRAN  seqmenc 
was  completely  rewritten  tn  order  to  handle  all  the  various  cases 
to  be  treated;  the  new  FORTRAN  segment  has  about  three  times  the 
number  of  lines  of  Gherson’ s.  The  new  segment  has  about  900  lines 
with  numerous  comments. 

In  some  of  our  recent  test  cases,  the  SEQUENCE,  segment  was: 

SEQUENCE.  LOOP  = 6 

DISCRETIZATION 

INDEXING 

SOLUTION 

OUTPUT 

TEST 

The  statement  LOOP  = 6 instructs  the  ELLPACK  system  to  execute  the 
SEQUENCE  six  times.  As  with  Gherson’s  SUBROUTINE  TEST,  ours  set 
initial  values  of  parameters  with  DATA  statements;  some  of  these 


specify  values  such  as  the  location  of  the  endpoints  of  the  interior 
insulating  strips,  the  value  of  EFOLO,  and  so  on;  others  control 
the  operation  and  they  are  "switches",  for  example  the  LOGICAL  variable 
MAGFON  (MAgnetic  Field  ON)  is  set  to  .TRUE,  when  B is  to  be  nonzero 
and  it  is  set  to  .FALSE,  when  B is  to  be  zero.  During  the  execution 
of  the  SEQUENCE  the  odd  numbers  of  times  (1 ,3,5),  MAGFON  is  set  to 
.TRUE,  and  the  approximation  U gives  estimates  of  the  potential 
for  a fixed  value  of  load  factor  »ij ; a different  value  of  EFOLD 
is  used  so  that  <J>  for  is  obtained  for  three  different  values 
of  EFOLD.  During  the  even  number  of  times  (2,4,6),  MAGFON  is  set  to 
•FALSE.,  the  load  factor  is  set  to  unity,  and  the  approximation  U 
gives  estimates  of  the  potential  ¥.  The  values  of  the  integrals 
in  (4-2)  are  estimated  by  the  Trapezoid  Rule,  values 
of  efficiency  for  n = 0,1/40,2/40,...  is  computed,  and  a graph 
of  efficiency  versus  load  factor  is  generated.  Sample  output  is  shown  below. 


7.  Discrete  model.  The  ELLPACK  module  5-POINT  STAR  which  we 
modified  uses  finite  difference  approximations  to  derivatives.  There 
are  three  basic  ones  which  we  display  below  in  terms  of  a given  function 
u of  z : 

d2u(0)/dz2  = [u(-Az)  - 2u(0)  + u(az)]/Az2  + (Az2/12)d^u(0)/dz4  •*■... 
du(0)/dz  = O3u(0)  + 4u(Az)  - u(2az)]/2az  - (Az2/3)d3u(0)/dz3  + ... 
du(0)/dz  = [u(-2Az)  - 4u(Az)  + 3u(0)]/2Az  - (Az2/3)d3u(0)/dz3  + ... 

Use  of  the  values  of  the  divided  differences  as  approximations  to  the 

values  of  the  derivatives  leads  to  local  error  which  is  the  order  of 

2 

Az  . Provided  that  the  u is  smooth,  then  when  one  halves  A z one 
expects  the  error  to  be  approximately  quartered.  This  applies  to 
finite  difference  approximations  to  ellptic  partial  differential 
equations  for  sufficiently  small  Az  and  provided  the  solution  has 
continuous  fourth  derivatives  in  the  region  and  at  the  boundary. 

For  the  mathematical  models  decribed  above  for  the  potentials  $ 
and  Y,  the  solutions  have  derivatives  with  singularities  at  the 
end  of  the  electrode  and  at  the  ends  of  insulator  vanes.  Consequently, 
one  expect  that  the  error  does  not  decrease  as  the  square  of  the 
mesh  spacing  as  this  spacing  tends  to  zero.  The  error  probably 
decreases  as  the  first  power  of  the  mesh  spacing  as  It  tends  to  zero. 

We  now  describe  the  approximations  used  in  5-POINT  STAR  and 
their  modifications  which  allow  the  slit-region  problems  to  be  treated. 
We  use  some  of  the  variable  names  of  ELLPACK  and  our  routines,  such 
as  TEST  so  that  this  report  will  also  serve  as  partial  documentation 
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for  the  programs  involved. 

The  horizontal  sides  of  the  rectangle  are  AX  < X < BX  and  the 

vertical  sides  are  AY  < Y < BY.  Numerical  values  for  AX,  BX,  AY, 

and  BY  are  user  supplied  in  the  BOUNDARY,  segement  of  the  ELLPACK 

program.  In  our  case,  AX  = AY  = 0.,  BX  = X_  . and  BY  = 1 . 

max 

The  number  of  vertical  and  horizontal  mesh  lines  are  NGRIDX  and 
NGRIDY,  respectively.  Numerical  values  of  these  are  specified  in 
the  GRID  segment  of  the  ELLPACK  program.  Uniform  mesh  spacing  is 
used  by  5-POINT  STAR  and  the  spacing  are  HX  and  HY  defined  by 

HX  = (-AX+BX )/{NGRIDX-l ) , 

HY  = ( - A Y+B Y ) / ( NGR I DY- 1 ) . 

We  denote  the  mes.i  points  by  (Xj.Yj)  and  the  values  are  stored  in 
the  ELLPACK  arrays  GRIDX  and  GRIDY.  The  values  are 

Xj  = AX  + (I-l)HX  = GRIDX(I),  I = 1 ,... .NGRIDX, 

Yj  = AY  + (J-l)HY  = GRIDY(J) , J = 1 .... .NGRIDY. 

The  approximation  to  the  solution  of  the  partial  differential 
equation  is  denoted  by  U and  this  can  be  considered  as  a two-dimensional 
array  and,  depending  on  which  problem  is  solved, 

U(I,J)  = -XXj.Yj)  or  U(I,J)  « H-ttj.Yj). 

At  each  interior  mesh  point,  the  Laplace  equation  is  approximated 
with  its  usual  divided  central  difference  approximation  (this  is 
modified  at  certain  points  as  explained  below): 
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[U(I-1,J)  - 2U(I,J)  + U( 1+1 ,J)]/fu 

(6-la)  + CU(I,J-1)  - 2U ( I , J ) + U(I,J+1)]/HY2  =0,  I = 2,. . . .NGRIDX-l , 

J = 2,..., NGRIDX-l. 


The  vertical  mesh  line  number  of  the  endpoint  of  the  electrode  (conductor) 
is  denoted  by  ICONDR  and  its  value  is  set  in  a DATA  statement  in 
SUBROUTINE  TEST.  Since  AX  = 0,  the  X-coordinate  of  the  endpoint  of 
the  conductor  is 


Xg  = ( IC0NDR-1 )HX  = GRIDX( ICONDR). 

The  LOGICAL  variable  HAVANE  (HAve  VANE)  is  set  to  .TRUE,  or 
.FALSE,  in  a DATA  statement  in  TEST.  When  it  is  .TRUE.,  then  there 
is  an  insulator  strip  (vane)  on  the  X-axis,  and  in  this  case  the 
horizontal  mesh  line  numbers  must  also  be  set  by  a DATA  statement 
in  TEST.  The  endpoints  are  denoted  by  IVANEL  (L  for  "left") 
and  IVANER  (R  for  "right").  The  coordinates  of  the  endpoints  are 
then  determined  by  the  modified  5-POINT  STAR  to  be  at  (X  ^’^,0) 
and  (X^R’^,0)  where 

X(L,1)  = (IVANEL-1 )HX  = GRIDX( IVANEL), 

X^  = ( IVANER-1  )HX  = GRIDX(  IVANER). 

The  boundary  conditions  on  the  boundary  of  the  rectangle  are 
given  below  in  the  same  order  as  in  (3-3)  [see  also  (5-1)]: 

(6-lb)  U(I,NGRIDY)  = LOAD FA,  I = 1 ,. . . , ICONDR, 

(LOADFA  for  LOAD  FActor;  this  is  declared  REAL  and  set  in  a DATA  statment 
in  SUBROUTINE  TEST), 
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[U(I,NGRIDY-2)  - 4U( I .NGRIOY-1 ) + 3U(I ,NGRIDY)]/2HY 


(6-lc) 


= B(GRIDX(I)),  I * ICONDR+1 .... ,NGRIDX-1 , 


( 6-1 d) 

U(NGRIDX,J)  * 0, 

0 = 1 ,. . . jNGRIDY, 

(6-le) 

U(I.l)  - 0, 

I = 1,...,IVANEL-1  and 

I * IVANER+1 , . . . .NGRIDX-1 

[-3U(  1 ,1 ) + 411(1,2)  - U(  1 ,3) 3/2HY 

(6-1  e* ) 


= B(GRIDX(I)), 

[-3U(1 ,J)  + 4U(2,J)  - U(3,J)]/2HY 

(6-1  f ) 


* 0, 


I = IVANEL IVANER, 


J * 2,. . . .NGRIDY-l . 


The  total  number  of  algebraic  equations  in  (6-1)  is  equal  to  the  total 
number  of  mesh  points  NGRIDX*NGRIDY,  hence  there  is  one  equation  for 
each  mesh  point  and  thus  for  each  unknown  U(I,J). 

The  equations  (6-1)  give  the  set  of  algebraic  equations  generated 
by  5-POINT  STAR  and  which  are  used  in  the  case  that  there  are  no 
insulator  vanes  in  the  interior  of  the  rectangle.  We  now  consider 
ther  case  in  which  there  are  insulator  vanes  in  the  interior;  we  allow 
for  one  or  two  such  vanes  or  slits. 

The  horizontal  mesh  line  number  of  a slit  is  set  in  a DATA  statment 
in  TEST.  The  variable  names  are  JYSL1D  and  JYSL2D  (D  for  "down"). 
If  one  or  both  of  these  are  nonzero,  then  this  informs  the  system  that 
there  is  one  or  two  interior  slits  and  in  this  case  the  corresponding 
vertical  mesh  line  numbers  for  the  endpoints  must  be  set  in  a DATA 
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| statement  In  TEST.  The  endpoints  are  denoted  by 


IXSL1L,  IXSL1R  for  JYSL1D  and 
IXSL2L,  IXSL2R  for  JYSL2D, 


where  "L"  and  "R"  stand  for  "left"  and  "right".  The  boundary 
condition  at  the  slit  is  84/3Y  = B(X)  and  at  the  downward  side  of 
the  slit,  the  difference  equation  approximation  is 

[U(I,JY$LkD-2)  - 4U(I ,JYSLkD-l ) + 3U(I ,JSLkD)]/2HY 

(6-lg)  - B(GRIDX(I)),  I = IXSLkL IXSLkR 

for  k = 1 and/or  k = 2. 

The  value  of  34/3Y  is  continuous  across  a slit,  but  the  value 
of  4 is  discontinuous;  that  is 

!ime+0  + *U,Y(k^+E)3  / 0 for  X(L,k)  < X < X(R,k\  j 

I 

l 

Si  nee  the  value  of  U at  both  the  bottom  and  the  top  of  the  slit 
are  unknowns,  we  need  to  add  another  equation: 

[-3U(I,JYSLkD+l)  + 4U ( I , JYSLkD+2 ) - U(I ,JYSLkD+3)]/2HY 

( 6-1 h ) = B(GRIDX(I)),  I = IXSLkL IXSLkR 

for  k = 1 and/or  k * 2. 

The  mesh  points  (I,JYSLkD)  and  (I,JYSLkD+l)  ’ correspond  to  the 
same  point  j * YJYSLkD^  in  the  re9'ion*  Consequently,  in  order  to 
obtain  an  approximation  to  the  potential  in  a rectangular  region  with 


interior  slits,  we  double  the  number  of  horizontal  mesh  jioints 


along  each  interior  slit. 

This  is  done  by  a modification  of  5-POINT  st*p  uk  ,, 

. 3 ™1NT  STAR‘  When  the  module 

s ca  ed  to  be  executed,  it  calls  a subroutine  which  deters  whether 

"*  ther'  'nt'r,0r  SHts  - -««-  the  horizontal  mesh 

-e  coordinates.  W,  give  three  examples.  Suppose  that  on,  wants 

estimates  of  0 on  a mesh  with  horizontal  n*sh  spacing  HY  . M 

vertical  mesh  spacing  HX  « 1.  where  X . 20 

max  £U' 

Example  1:  no  Interior  slits. 

In  TEST  one  includes  the  statement 
MTA  JYSL1D,  JYSL2D  / 0,  0 / 

The  GRID,  segment  of  the  ELLPACX  program  Is 
GRIO.  UNIF0RM  X * 21  $ UNIFORM  Y » 11 
Then  2),  NGR1DY  “11,  HX-  20/20.,.  HY  * 1,10  ■ , 

T«mp,e  2:  one  Interior  s„t  with  endpoints  «..3,  .„d  (,0..3) 

In  TEST  one  Includes  the  statement 

DATA  JYSL10,  JYSL2D,  IXSL1L,  IXSL1R  / 4,  0,  5,  11  / 

The  GRID,  segment  of  the  ELLPACX  program  is 
GRID.  UNIFORM  X ■ 21  $ UNIFORM  Y - 12 

Then  NGRIDX  * 21,  NGRIDY  . 12,  and  the  ELLPACX  system  set 

HY  « 20/20  - 1,  HY  • l/n  „ 0909090  i*  i 

x .uy09090...  It  also  generates 

values  for  GRlox  ,„d  GRIOY  Wlth  spacing  .„d  ,/n> 
respectively.  The  modified  5-POINT  STAR  determines  that 
0 ts  nonzero  so  that  there  Is  an  Interior  slit  It 
^ ^ ‘^e  is.  pair  of  horizon,.,  mesh  lines  which 

correspond  to  the  same  value  of  Y along  the  s', t.  It 

changes  HY  from  l/n  to  l/lo  « l .nH  it 

/,u  •’  «na  it  resets  the 
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the  values  of  the  Y coordinates  of  the  mesh  lines: 

Yj  = (J-l)HY  = GRIDY(J),  J « 1,..., JYSLlD 

yjysud+i  ‘ ¥jyslid  + Hy*10'5  * GMDY(JYSUIW) 

Yj  = (J-2)HY  = GRIDY(J),  J « JYSL1D+2 NGRIDY. 

The  values  in  GRIDY  are  then: 

J = 1,  2,  3,  4,  5,  6,  7,  8,  9,10,11,12 

GRIDY(J)  = 0.  ,.l ,.2,.3,.3+10”6,.4,.5,.6,.7,.8,.9,l. 

Example  3:  two  Interior  slits  with  endpoints  (4, .3) ,(10, .3)  and 
(3, .7), (5,. 7): 

In  TEST  one  includes  the  statements 

DATA  JYSLlD,  IXSL1L,  IXSL1R  / 4,  5,  11  / 

DATA  JYSL2D,  IXSL2L,  IXSL2R  / 9,  4,  6 / 

The  GRID,  segment  of  the  ELLPACK  program  is 
GRID.  UNIFORM  X = 21  $ UNIFORM  Y = 13 
The  ELLPACK  system  sets  HX  = 20/20  =1.,  HY  * 1/12  s .08333... 
and  generates  values  of  GRIDX  and  GRIDY  with  these  spacing. 

The  modified  5-POINT  STAR  finds  that  both  JYSLlD  and 
JYSL2D  are  nonzero  so  that  there  are  two  interior  slits. 

It  then  knows  that  there  are  two  pairs  of  horizontal  mesh 
lines  which  correspond  to  the  two  values  of  Y along  the  slits. 
It  changes  the  value  HY  from  1/12  to  1/10  s .1  and  it 

resets  the  values  of  the  Y coordinates  of  the  mesh  lines: 
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Yj  * (J-l)HY  * GRIDY(J),  J * 1 JYSL10 

VJYSL1D+1  * YJYSL10  + HY*'°‘5  ’ GRIDY(JYSLmi) 

Yj  = (J-2)HY  = GRIDY(J) , J * JYSL1D+2 JYSL2D 

VJYSL20fl  ’ VJYSL2D  + HV*,0‘5  * GRIDY( JYSL2D+1 ) 

Yj  = (J-3)HY  = GRIDY(J) , J = JYSL2D+2,. . . ,NGRIDY 

The  values  in  GRIOY  are  then: 

0=1,  2,3,4,  5,  6,  7,  8,  9,  10,  11,12,13 

GRIDY(J)  = 0.,.1,.2,.3,.3+10‘6,.4,.5,.6,.7,.7+10'6,.8,.9,1. 

By  introducing  these  lines  of  "double"  mesh  points,  we  have  added 
additional  unknowns,  namely  the  values  of  U( I , JYSLkD)  and  U(I,JYSLkD+l), 
I = 1,,.., NGRIDY.  Equations  corresponding  to  some  of  these  are  given  in 

(C-lg)  and  (6-lh),  specifically  for  I = IXSLkL IXSLkR.  At  points 

off  of  the  slits,  the  potential  is  continuous,  so  we  add  the  equations 

U( I, JYSLkD)  - U(I ,JYSLkD+l ) =0,  I = 1 ,. . . ,IXSLkL-l  and 

(6-1 i) 

I = IXSLkR+1 , . . . .NGRIDX 

In  equation  (6-la)  we  give  the  approximation  for  the  Laplace  equation; 
but  now  that  does  not  apply  for  (1,0)  on  these  lines  of  double  mesh 
points.  In  (6-la)  the  range  of  tha  subscripts  changes  from 

I = 2,. . . .NGRIDX-l , J = 2 NGRIDX-1 

to 


-180- 


I « 2,...,NGRIDX-1,  J • 2,. , . ,JYSL2D-2 

(6-la 1 ) 

J « JYSL1D+2,. . . .JYSL2D-2,  J • JYSL2D+2,. . . ,NGRIDY-1 

The  approximation  to  the  Laplace  equation  for  a mesh  point  with 
JYSLkD  and  at  a distance  at  least  2HX  from  the  end  of  a slit 
Is  taken  as  [note  the  "+2"  In  the  second  pair  of  square  brackets] 

[U( 1-1, JYSLkD)  - 2U{ I, JYSLkD)  + U(I+1 ,JYSLkD)]/HX2 

(6-1 j)  + [U(I.JYSLkD-l)  - 2U{ I .JYSLkD)  + U(I ,JYSLkD+2)]/HY2  • 0, 

I - 2,. . . ,IXSLkL-2,  and  I - IXLSkR+2,. . . .NGRIDX. 

For  mesh  points  Immediately  to  the  right  and  to  the  left  of  the 
end  points  of  the  slits,  we  incorporate  the  average  value  of  U 
at  the  top  and  bottom  of  the  end  of  the  slit  and  use  [again,  note  the  "+2"] 

[U(I-1, JYSLkD)  - 2U( I, JYSLkD)  + U(I+1 ,JYSLkD)/2  + U(I+1 ,JYSLkD+l )/2]/HX2 

(6-1 k)  + [U(I.JYSLkD-l)  - 2U( I, JYSLkD)  + U(I ,JYSLkD+2)]/HV2  - 0, 

for  I - IXSLLL-1 


and 

[U(I-l,JYSLkD)/2  + U( 1-1 ,JYSLkD+l )/2  - 2U( I, JYSLkD)  + U(I+1 ,JYSLkD)]/HX2 
(6-U)  + [U(I.JYSLkD-l)  - sU(I, JYSLkD)  + U(I,JYSLkD+2)]/HY2  - 0 

for  I - IXSLkR+1 

These  last  two  equations  Involve  a linear  combination  of  six 
of  the  unknowns;  the  usual  approximation  to  Laplace's  equation  involves 
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five  as  In  (6-la).  The  ELLPACK  module  5-POINT  STAR  constructs  a 
set  of  difference  equations  vhich  has  at  most  five  of  the  unknowns 
In  them;  consequently,  the  ELLPACK  system  does  not  allocate  a 
sufficient  amount  of  storage  when  the  DISCRETIZATION,  segement  of 
the  ELLPACK  program  Is 

DISCRETIZATION.  5-POINT  STAR 

In  order  to  obtain  the  necessary  storage,  we  Include  a second 
discretization  module--th1s  is  not  subsequently  used,  its  purpose 
is  merely  to  obtain  extra  storage.  Hence  in  an  ELLPACK  program 
which  estimates  the  potential  of  a region  with  slits,  one  uses 

D I SC  RET I ZAT I ON ( 1 ) 5-POINT  STAR 

DISCRETIZATION^)  HOLR  9-POINT  (I0RDER-40) 

and  in  the  SEQUENCE,  segment,  DISCRETIZATION  is  replaced  with 
DISCRETIZATION(l). 

Caution:  The  locations  of  the  end  of  the  conductor,  and 

the  endpoints  of  the  slits  and  vane  are  given  in  terms  of  mesh 
line  numbers.  Consequently,  if  the  number  of  grid  lines  is  changed 
in  the  ELLPACK  GRID,  then  the  integer  values  ICONDR,  JYSL1D,  and 
so  on,  must  be  changed  Inorder  to  obtain  the  same  geometry. 

Caution:  Because  a three  term  approximation  is  used  for 

the  normal  derivative  at  an  insultor,  for  accuracy  one  must  have 
at  least  two  interior  (non-slit)  mesh  lines  between  any  pair 
of  insulators.  In  Example  3 above,  with  two  slits,  note  that 
there  are  two  interior  mesh  lines  between  the  lower  slit  and  Y * 0 
(J  * 2,3)  and  two  between  the  upper  slit  and  Y - 1 (J  *=  11,12) 
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and  these  are  the  minimal  number.  There  arc-  three  between  the 
pair  of  Interior  slits  (J  3 6,7,B)  and  this  is  one  more  than 
the  required  minimum  number. 

Caution:  The  derivatives  of  the  potential  are  singular 

at  the  end  point  of  the  conductor  and  at  the  end  of  the 
insulator  slits  and  vane.  One  should  run  a few  test  cases  to 
determine  the  mesh  spacing  which  yieldssufficient  accuracy. 

In  particular,  the  vertical  mesh  line  spacing  should  be  smaller 
than  the  length  of  the  electrode. 
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APPENDIX  D 


MAJOR  ACCOMPLISHMENTS  OF  THE  COMPLETE  PROGRAM 

This  report  covers  the  final  ye..:  of  an  eight-year  experimental 
and  theoretical  program  aimed  at  developing  practical  two-phase 
LMMHD  generators.  Thus,  the  significant  accomplishments  are  summa- 
rized below  by  subject  area--generator  experiments,  code  develop- 
ments, facility  developments,  local  measurements,  foams,  and 
analytical  studies.  The  goals  of  this  program  were  to  study  the 
basic  phenomena  and  to  demonstrate  that  practical,  high-efficiency 
LMMHD  generators  can  be  built,  thus  making  possible  the  development 
of  high-efficiency  LMMHD  energy-conversion  systems. 

Generator  Experiments 

The  four  generator  channels  tested  in  this  program  are  compared 
in  Fig.I.l.  The  basic  channel  dimensions  were  set  by  the  facility 
constraints--magnet  gap  length  and  pole-face  size,  available  flow 
rates,  etc. --and  the  designs  were  done  within  these  constraints, 
using  the  generator  design  code  extracted  from  the  system  code 
(Section  II.  1).  The  structural  design  for  the  first  three  channels 
used  the  same  copper  electrodes  and  compensating  bars  with  different 
side  walls  to  set  the  different  channel  contours.  For  LT-4,  an 
improved  design  with  a filament-wound  liner  was  developed  to  provide 
for  different  contours  with  a higher  pressure  rating  and  easier 
sealing. 

Diagnostic  instrumentation  has  been  an  integral  part  of  this 
program  since  its  initiation.  Gas  and  liquid  volume  flow  rates, 
temperatures,  pressures,  void  fractions  (by  means  of  Y-ray  attenua- 
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tion),  electrode  voltages,  load  currents,  and  magnetic  flux  densities 
have  been  measured  for  all  channels.  For  channels  LT-2,  LT-3,  and 
LT-4,  probes  were  built  into  one  insulating  wall  to  measure  voltage 
profiles  between  the  electrodes.  As  the  scope  of  the  program 
permitted,  the  measurement  capability  was  improved.  Recently  local 
measurements  were  initiated,  as  described  below. 

The  first  channel,  IT-1,  was  designed  for  a slip  ratio  of 

I 

unity,  because  the  actual  slip  ratio  or  bubble  parameter,  C^,  (based 
on  experimental  data)  was  unknown.  As  anticipated,  slip  ratios 
greater  than  unity  were  observed,  with  some  values  in  excess  of  two, 
and,  thus,  the  liquid  velocity  along  the  generator  decreased. 

I 

Experimental  values  of  were  determined. 

Next,  a constant-area  channel,  LT-2,  was  tested  to  obtain 

I 

further  basic  data  on  slip  ratios  and  values  of  C^.  In  general, 

l 

Cd  values  were  less  for  LT-2  than  for  LT-1.  The  efficiency  was 
higher  for  LT-2  than  for  LT-1  at  zero  or  low  void  fractions,  but 
lower  at  the  higher  void  fractions.  This  was  expected,  because 
LT-2  would  yield  liquid  velocities  that  were  almost  constant  for 
zero  or  low  void  fractions,  but  not  for  high  void  fractions. 

(Constant  liquid  velocity  is  desired  because  it  results  in  constant 
generated  voltage  and  no  or  minimal  circulating  current  loss.) 

Analysis  of  the  data  from  LT-1  and  LT-2  showed  that  much  of  the 
discrepancy  between  the  experiments  and  the  models  could  be  explained 
by  the  existence  of  a pure-liquid  shunt  layer  on  the  insulating 
walls,  where  the  shunt  layer  bypasses  some  of  the  current  around 
the  load.  To  eliminate  this  possible  shunt  layer,  gas  injection 
along  the  insulating  walls  was  proposed.  Gas  injection  was  tested 
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wUh  channels  LT-2  and  LT-3.  For  LT-3,  substantially  better  perfor- 
mance was  obtained  without  the  injectors.  Support  studies  showed 
that  performance  improved  as  the  amount  of  gas  injected  was  decreased, 
that  the  injected  gas  did  not  form  the  desired  pure-gas  wall  layer,  and  that 
the  partial  flow  obstruction  of  the  gas  injection  ports  degraded  performance. 
A later  analytical  study  has  shown  that  liquid  shunt  layer  losses  will  not 
be  significant  for  large  generators. 

The  second  diversing  channel,  LT-3,  was  designed  on  the  basis  of  data 
obtained  from  LT-1  and  LT-2.  Improved  generator  performance  was  obtained, 
with  efficiencies  of  up  to  0.50  at  high  void  fractions.  Channel  LT-3  was 
also  tested  at  higher  liquid  velocities  in  the  revised  facility  (described 
below),  and  the  reduction  in  slip  ratio  with  increasing  liquid  velocity  was 
demonstrated.  Substantial  increases  in  generator  efficiency  were  obtained 
with  higher  liquid  flow  rates. 

Experiments  to  demonstrate  that  good  liquid-to-gas  heat  transfer 
exists  in  the  generator  were  succssfully  completed  with  LT-3.  Good 
heat  transfer  is  essential,  becuse  it  is  the  almost-constant-tempera- 
ture  expansion  of  the  gas  (vapor)  in  the  generator  that  yields  the 
higher  system  efficiencies  for  LMMHD  power  cycles. 

The  design  of  LT-4  and  the  significant  test  results  obtained 
are  described  in  this  report,  ^ower  densities  of  up  to  32  MWe/m^ 
were  obtained  in  this  small  20-kWe  generator,  equal  to  or  better 
than  that  anticipated  for  LMMHD  generators  in  practical  power 
systems.  Efficiencies  higher  than  0.60  at  average  void  fractions 
well  above  0.50  were  measured.  The  rapid  reduction  in  slip  ratio 
with  increasing  liquid  velocity  was  confirmed  for  high  power  densities. 
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Testing  of  HT-1  was  also  completed  under  open-circuit  conditions 
over  the  temperature  range  of  approximately  490  to  approximately  740  K 
[Dunn  et  al  (1979)].  The  slip  ratio  was  shown  to  decrease  with  increasing 
temperature,  implying  higher  generator  and  system  efficiencies;  this 
anticipated  result  was  a prime  reason  for  performing  the  experiments. 

Code  developments 

The  channel  design  code  (GENDES)  was  developed  prior  to  this 
program  as  part  of  the  system  performance  code  [Petrick  et  al 
(1970)].  In  that  code  (Section  II. 1)  the  liquid  velocity  is  held 
constant  and  the  channel  contour  is  calculated,  so  that  it  is  of  no 
value  in  analyzing  experimental  data.  Three  new  codes  were  developed 
in  connection  with  the  experiments:  ANLSIS,  MHDPRE,  and  MHD9. 

ANLSIS  is  the  basic  code  used  to  process  all  of  the  raw  experi- 
mental data.  From  that  data,  it  calculates  generator  performance 
and  the  variation  of  the  parameters  (pressure,  liqu'd  velocity,  gas 
velocity,  etc.)  along  the  channel. 

MHDPRE  is  used  to  predict  generator  operation  for  given  inlet 
conditions.  The  fluid  equations  are  solved,  and  the  load  voltage 
and  input  pressure  are  varied  to  match  the  desired  load  resistances 
and  exit  pressure.  Rather  than  solving  Laplace's  equation,  an  end 
resistance  is  assumed. 

MHD9  is  the  most-detailed  code;  it  solves  both  the  fluid 
equations  and  Laplace's  equation  iteratively  until  convergence  is 
obtained.  This,  an  improved  electrical  model  is  used,  and  no  end 
resistance  is  assumed.  The  inlet  conditions  are  specified,  and  the 
load  voltage  varied  to  match  a specified  condition.  MHP9  is  an 
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expensive  code  to  run,  because  of  the  detail  and  because  convergence 
Is  not  always  obtained. 

Facility  Development'- 

The  Nak-nitrogen  ambient -temperature  facility,  developed  to 
test  two-phase  LMMHD  generators,  provides  known  flows  of  liquid 
(NaK)  and  gas  (nitrogen)  to  the  generator.  It  consists  of  a NaK 
loop  and  a nitrogen  loop  with  a common  mixer,  generator  test  section, 
and  separator /NaK  storage  tank.  The  NaK  is  recirculated,  and  the 
nitrogen  is  exhausted  to  the  atmosphere  front  the  separator  tank 
through  a secondary  separator.  The  nitrogen  exhaust  line  is  equipped 
with  a control  valve  and  controller  to  sense  and  fix  the  generator 
exit  pressure.  The  NaK  loop  includes  three  canned-rotor  pumps  in 
series  with  bypass  valves,  so  that  only  one  or  two  can  be  operated 
if  desired,  a turbine  flow  meter,  control  valves,  and  pressure 
gauges.  The  nitrogen  loop  comprises  a heater  and  control  system  to 
heat  the  gas  to  the  liquid  temperature,  a flowmeter,  control  valves, 
and  pressure  gauges.  A dc  electromagnet  provides  the  magnetic 
field. 

Throughout  this  program  the  facility  has  been  continuously 
upgraded,  to  improve  its  characteristics  and  controlabi 1 ity,  to 
provide  a wider  range  of  operationing  parameters,  and  to  provide 
improved  instrumentation.  In  1975-1976,  the  facility  was  completely 
rebuilt  and  most  major  components  replaced,  to  allow  operation 
at  higher  pressures  and  flow  rates.  Present  capacities  are  18  kg/s 
(330  gpm)  NaK  and  0.44  kg/s  (800  scfm)  nitrogen  maximum  at  maximum 
system  pressure  of  2.0  MPa  (275  psig)  and  maximum  fluid  temperature 
of  328  K (130*F) . 
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The  test  section  fits  between  conventional  flanges,  and  is 
inserted  into  the  (adjustable)  air  gap  of  the  electromagnet. 
Instrumentation  is  provided  to  measure  and  control  gas  and  liquid 
flow  rates,  pressure.,  and  temperatures,  and  to  measure  (dc)  generator 
voltages  and  currents,  pressures  along  the  generator  channel,  magnetic1 
field  strengths,  and  void  fraction  distributions  along  the  generator. 
Auxiliary  equipment  includes  1)  instrumentation  for  recording  experi- 
mental data;  2)  a dc  electrogmagnet  with  a maximum  flux  density  greater 
than  1.2  T at  0.09  m (3.5  in.)  air  gap  and  a working  volume  0.41  m 
(16  in.)  in  diameter  by  an  adjustable  air-gap  length;  3)  a Y-ray 
void-fraction  measuring  system,  and  4)  a hot-film  anemometry 
system. 


Local  Measurements 

Initial  measurements  have  been  made  of  local  flow  parameters  in 
a NaK-nitrogen  two-phase  LMMHD  generator,  using  hot-film  and  resisti- 
vity probes.  Although  the  results  are  preliminary,  there  is  evidence 
that  large  electromagnetic  force  densities  have  a significant  effect 
on  the  size  and  distribution  of  the  gas  bubbles  or  voids. 


Foams 

The  use  of  surface-active  agents  in  the  liquid  metal,  to 
produce  a more-stable  foam  flow  in  the  LMMHD  generator,  reduces  the 
slip  loss,  allows  the  generator  to  operate  at  higher  average  void 
fractions,  and  results  in  a more-efficient  conversion  system.  A 
program  to  evaluate  surface-active  agents  in  liquid  metals,  especialy 
NaK,  was  initiated  in  1976. 
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The  feasibility  of  generating  relatively  stable  bubbles,  hence, 
a foam,  in  liquid  metals  was  demonstrated.  These  phenomena  were  documented 
in  both  still  and  motion  pictures.  Surface  tension  measurements  and  foaming 
experiments  showed  that  viscosity  is  also  a factor  in  promoting  bubble 
formation  and  persistence.  Wetting  and  contact  angle  measurements  were  made 
for  stainless  steel  and  carbon  steel  immersed  in  eutectic  NaK. 

Barium  has  been  identified  as  an  attractive  additive  for  the 
generation  of  foams  with  NaK,  and  bismuth  is  a back-up  for  barium.  A 
guiding  hypothesis  for  understanding  the  phenomena  contributing  to 
foamability  in  liquid  metal  systems  has  been  developed,  and  basic 
data  on  surface  tension  obtained.  The  next  step  is  to  test  foams  in 
the  LMMHD  generator. 

Analytical  Studies 

In  the  early  phase  of  the  program,  the  analytical  work  was 
concentrated  on  developing  the  required  comptuer  codes  and  under- 
standing the  basic  loss  mechanisms.  Subsequently,  this  work  was 
extended  to  more-detailed  analytical  studies. 

The  annular  generator  geometry  has  certain  advantages,  especialy 
the  lack  of  insulating  walls  and  the  ease  of  fabrication.  The 
equations  were  developed  and  solutions  obtained  for  three  cases--con- 
stant  velocity  and  no  armature  reaction,  laminar  flow  with  no  armature 
reaction,  and  armature  reaction  with  constant  velocity.  Numerical 
examples  showed  that  1)  the  attainable  terminal  voltages  appear  to 
be  very  low,  2)  flow  reversal  and  large  viscous  loss  occur  at  or 
below  the  desired  power  densities,  and  3)  armature  reaction  effects 
are  important  and  compensation  techniques  appear  impractical.  Thus, 
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this  annular  geometry  does  not  appear  attractive  for  either  generator 
or  pump  operation. 

An  analytical  study  of  the  liquid  shunt  (wall)  layer  sizes  and 
losses  showed  that  these  losses  are  not  expected  to  be  significant 
for  large  generators,  i.e.,  less  than  1.0  X decrease  in  efficiency 
for  Hartmann  numbers  in  excess  of  1000. 

End  current  losses  in  LMMHD  generators  can  be  controlled  by  the 
proper  choice  of  1)  channel  aspect  ratio,  2)  magnetic  flux  density 
distribution  outside  of  the  electrodes,  and  3)  number  and  dimensions 
of  insulating  vanes.  A computer  code  has  been  developed  to  calculate 
the  end  current  loss  as  the  above  parameters  are  varied.  Initial 
studies  with  this  code  have  shown  that  generator  efficiencies  in 
excess  of  0.8  are  attainable. 


Conclusion 

The  results  obtained  have  provided  basic  data  on  two-phase 
flews  in  a magnetic  field  and  demonstrated  that  high  power  density 
and  good  efficiency  are  attainable  in  small  LMMHD  generators.  All 
indications  are  that  large,  high-efficiency  LMMHD  generators  can  be 
built. 
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